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Abstract 



Microquasars, X-ray binary systems that generate relativistic jets, were discovered 
in our Galaxy in the last decade of the XX th century. Their name indicates that 
they are manifestations of the same physics as quasars but on a completely different 
scale. 

Parallel to this discovery, the EGRET instrument on board of the Compton 
Gamma Ray Observatory detected 271 point like gamma-ray sources 170 of which 
were not clearly identified with known objects. This marked the beginning of 
gamma-ray source population studies in the Galaxy. 

We present in this thesis models for gamma-ray production in microquasars 
with the aim to propose them as possible parent populations for different groups of 
EGRET unidentified sources. These models are developed for a variety of scenarios 
taking into account several possible combinations, i.e. black holes or neutron stars 
as the compact object, low mass or high mass stellar companions, as well as leptonic 
or hadronic gamma-ray production processes. 

We also show that the presented models for gamma-rays emitting microquasars 
can be used to explain observations from well known sources that are detected in 
energy ranges other than EGRET's. Finally, we include an alternative gamma-ray 
producing situation that does not involve microquasars but a specific unidentified 
EGRET source possibly linked to a magnetized accreting pulsar. 
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This thesis deals with the possibility to link two of the most intriguing and recent 
discoveries of the astrophysics from the last decade of the XX th century: micro- 
quasars and THE GALACTIC POPULATION OF GAMMA-RAY SOURCES. 

The first complete map of the Galaxy in gamma-rays was obtained with the 
ESA COS-B satellite 1 , launched on August 9 th , 1975, and operating until April 
25 th , 1982. It was, however, only with the advent of the Compton Gamma Ray 
Observatory (CGRO) in the 1990s that population studies of gamma-ray sources 
properly started. EGRET instrument, on board of CGRO, detected 271 point like 
gamma-ray excesses with a significance of more than 5cr on the galactic plane and 
more than Aa out of it. On the other hand, X-ray binary systems able to generate 
relativistic radio jets were baptized as "microquasars" with the discovery of the 
sources 1E1740. 7-2942 and GRS 1758-258 in 1992, soon after the beginning of the 
"Compton Era". 

We will present in this thesis models for gamma-ray production in microquasars 
(MQs). These models will be developed through the following chapters, where the 
link with gamma-ray sources will be discussed for different scenarios. 



1.1 Microquasars 

When X-ray binary systems (XRB) are able to generate relativistic radio jets, we 
call them microquasars. The name was given by Mirabel et al. (1992) as a clear 
reference to the already quite studied phenomena associated with quasars. Quasars 
and microquasars share three common basic features: 

1) A spinning compact object (in the case of MQs it can be either a neutron star 
(NS) or a black hole (BH)). 

2) An accretion disk. 

1 COS-B carried a single large experiment, the Gamma-Ray Telescope. The energy range cov- 
ered was: 2 keV - 5 GeV. It was originally projected to last two years, but it operated successfully 
for 6 years and 8 months. 
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3) Collimated jets of relativistic particles. 

However, as it is indicated by their names, quasars and microquasars, are dif- 
ferent manifestations of the same physics at different scales. We can highlight some 
of their differences in the following table: 



Feature 


Quasars 


Microquasars 


Black hole mass 


Several million 
solar masses 


Few solar masses 


Accretion disk size 


~ 10 9 km 


~ 10 3 km 


Mean thermal T of 
the accretion disk 


Several thousand 
degrees 


Several million 
degrees 


Characteristic A 
of the disk radiation 


ultraviolet and 
optical 


X-rays 


Distance traveled 
by the jets 


Millions 
of light-years 


Few light-years 



All this information is summarized in a schematic way in Fig. 11.11 We can see 
from the table that the linear and time scales are proportional to the black hole 
mass (Sams et al. 1996). In fact, quasars have super-massive black holes located in 
the center of distant galaxies (active galactic nuclei, AGN), while MQs have stellar- 
mass black holes forming binary stellar systems and can be found in the Galaxy. In 
quasars, the accretion disk feeds itself from disrupted stars or from the interstellar 
medium of the host galaxy; in MQs, instead, the matter supply is from the com- 
panion star in the binary system. 

Taking into account the nature of the donor star, MQs are classified into high- 
mass microquasars (HMMQs) and low-mass microquasars (LMMQs). HMMQs have 
a young massive stellar companion and the mass transfer is done mainly through 
the stellar wind. On the other hand, LMMQs contain an old stellar companion that 
transfers mass by Roche lobe overflow. 

The presence of a relativistic jet is a fundamental characteristic in quasars as 
well as in MQs. The first evidence of a jet-like feature in astrophysics was discov- 
ered by Curtis (1918): he found the optical jet from the elliptical galaxy M87 in 
the Virgo cluster. The finding that jets can also be produced at smaller scales by 
binary stellar systems is much more recent. In the late 1950s several radio sources 
were matched with very dim optical objects that looked like stars, but had strange 
spectra with a strong ultraviolet excess. They were called quasars for "quasi-stellar 
radio source", and it was in the early 1960s that the first spectrum lines of these 
"stars" were identified showing a redshift of 0.158 (it was the case of 3C273). The 
discovery of such a large redshift established the cosmological nature of quasars. 
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Figure 1.1: Quasars and Microquasars are supposed to be different scales of the same 
physics (from Mirabel & Rodriguez 1998). 



Relativistic jets in galactic objects were first observed by Margon et al. (1979) 
in SS433. This binary remained the only known object of this kind for more than 
a decade. As it was already mentioned above, supermassive black hole disks emit 
strongly at optical and UV wavelengths. The more massive the black hole, the 
cooler the surrounding accretion disk. For a black hole accreting at the Eddington 
limit, the characteristic blackbody temperature at the last stable orbit in the sur- 
rounding accretion disk is given approximately by T ~ 2 x WM' 1 ^ (Rees 1984), 
with T in Kelvin and the mass of the black hole, M, in solar masses. This explains 
why whereas accretion disks in quasars have strong emission in the optical and 
ultraviolet, BH and NS binaries usually are identified the first time through their 
X-ray emission. It is therefore understandable that the jet phenomenon was first 
discovered far away in very distant galaxies, and that, there was an impasse in the 
discovery of stellar size sources with relativistic jets in our galactic surroundings 
until the recent developments in X-ray astronomy. It has to be noted that among 
galactic jet sources, SS 433 is unusual, because of its broad optical emission lines 
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Figure 1.2: Radio counterpart from microquasar 1E1740. 7-2942. Notice the clear two- 
sided jet feature (from Mirabel et al. 1992). 

and its brightness in the visible that made it detectable before such a technological 
development took place. 

As we remarked some lines above, MQs are firstly detected by their X-ray emis- 
sion as X-ray binaries. A subsequent multi-wavelength study of the system can 
establish the presence of a radio counterpart, which is the unmistakable signature 
of the non-thermal emission from a jet composed by relativistic particles. The con- 
fluence of observations made in hard X-rays and in radio wavelengths gave as a result 
the discovery, in 1992, of two such stellar sources of relativistic jets in the galactic 
center region: 1E1740.7-2942 and GRS 1758-258 (Mirabel et al. 1992, Rodriguez et 
al. 1992). The radio counterpart of 1E1740. 7-2942 is shown in Fig. 11.21 Based on 
the spectral shape, it was assumed that the steady radio emission in X-ray binaries 
comes from jets (Rodriguez et al. 1995; Fender et al. 1999). This assumption 
was corroborated recently with VLBA observations of GRS 1915+105 (Dhawan et 
al. 2000), Cyg X-l (Stirling et al. 2001) and other sources. The resu lting VLBA 
images of GRS 1915+105 are in agreement with the model of a conical expanding 
jet emitting through synchrotron processes (Hjellming & Johnston 1988; Falcke & 
Biermann 1999). It is the high brightness temperature, the rapid variability, and 
the linear polarization which suggest that the observed radio emission is due to 
synchrotron processes, from the relativistic leptons in the jet. 
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Another very important parallelism can be established between quasars and 
MQs. Apparent superluminal 2 transverse motions in the expansion of jets have 
been observed in quasars (Pearson & Zensus 1987; Zensus 1997) with values up 
to ten or more times the speed of light. Shortly after the quasar - microquasar 
analogy was proposed, the first superluminal source in the Galaxy was discovered: 
GRS 1915+105 3 (Mirabel & Rodriguez 1994). Thereafter, this characteristic was 
observed in MQs like GRO J1655-40 4 (Tingay et al. 1995; Hjellming & Rupen 1995) 
and XTE J1748-288 5 (Hjellming et al. 1998), among others. 

The ejection of two radiating components into opposite directions implies (if 
the viewing angle is smaller than 90°) that one of them will be approaching the 
observer and the other one, receding. The synchrotron emission from the ejecta 
will be affected by the Doppler effect. For example, in the case of GRS 1915+105, 
the approaching condensation shows an apparent velocity on the plane of the sky of 
1.25c and the receding one of 0.65c with the first one appearing brighter than the 
second (see Fig. II. 3|) . Using VLA data, it was inferred that the ejecta moves with 
a speed of 0.92c at an angle = 70° to the line of sight. We can then assume that 
the existence of apparent superluminal movements implies relativistic bulk motions 
in the jets. 

As well as there are relativistic bulk motions, recently, X-ray synchrotron emis- 
sion from the jet of XTE J1550-564 has been observed (Corbel et al. 2002, Kaaret 
et al. 2003), which implies the presence of extremely relativistic electrons with TeV 
energies and shock reacceleration in the outer jet (see Fig. 12. 8j) . 



1.2 Gamma- Ray Sources 

The study of gamma-ray emission from unidentified sources was the first motiva- 
tion of this thesis, searching for an understanding of the origin of some mysterious 
EGRET detections. 

EGRET, which stands for Energetic Gamma Ray Experiment Telescope, was 
one of the four instruments on board of the Compton Gamma-Ray Observatory. 
This satellite was launched on April 5 th 1991 aboard the space shuttle Atlantis; on 
June 4 th 2000 it was safely deorbited and re-entered the Earth's atmosphere. The 
main characteristics of the four instruments (that can be seen on board in Fig. II .4j) 
were: 

• BATSE was optimized to measure brightness variations in gamma-ray bursts 
and solar flares on timescales down to microseconds, over the energy range 30 

2 See Appendix |0 for a detailed explanation of apparent superluminal motions. 
3 Found with the satellite GRANAT (Castro-Tirado et al. 1994; Finogenov et al. 1994). 
4 Found with the satellite CGRO (Zhang et al. 1994). 

5 Found with the satellite X-Ray Timing Explorer (XTE) (Smith ct al. 1998). 
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Figure 1.3: Bright radio condensations ejected in opposite directions. The approaching 
component (left) shows an apparent velocity on the plane of the sky of 1.25c while the 
receding one (right) seems to move at 0.65c. Observations of GRS 1915+105 with VLA 
in 1994 (from Mirabel & Rodriguez 1994). 
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Figure 1.4: The Compton Gamma-Ray Observatory and its four instruments on board, 
BATSE, OSSE, COMPTEL, EGRET. The four scientific instruments on Compton were 
designed to operate simultaneously and cooperatively, to cover the gamma-ray energy range 
from 15 keVto 30 GeV. 

keV to 1.9 MeV. It was also continuously monitoring all transient sources 
and bright persistent sources in the gamma-ray sky. It consisted of 8 detectors. 

• OSSE was designed to undertake comprehensive spectral observations of as- 
trophysical sources in the 0.05 to 10 MeV range, with capability above 10 
MeV for solar gamma-ray and neutron observations. Each of the 4 detectors 
had a single axis pointing system which enabled a rapid OSSE response to 
target of opportunities, such as transient X-ray sources, explosive objects, and 
solar flares. 

• COMPTEL has performed the first sky survey in the energy range from 1 
to 30 MeV. Source mapping was provided over a field of view of about 1 
steradian. 

• EGRET was the highest energy instrument on Compton, and covered the 
broadest energy range, from 20 MeV to ~ 30 GeV. It had a wide field 
of view, good angular resolution and very low background. Because it was 
designed for high-energy studies, the detector was optimized to detect gamma 
rays when they interact by pair-production process, which forms an electron 
and a positron within the EGRET spark chamber. 

As we have already mentioned, we will focus our attention on the last instru- 
ment. The Third EGRET Catalog (3EG) (Hartman et al. 1999) of high-energy 
gamma-ray sources includes data from April 22 nd , 1991, to October 3 rd , 1995. The 
271 sources that conform the catalog were detected at E > 100 MeV. They include: 

• The single 1991 solar flare bright enough to be detected as a source. 

• The Large Magellanic Cloud. 

• 5 pulsars (their number has been recently extended up to 7). 
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Figure 1.5: This diagram indicates the positions on the sky of EGRET sources of celes- 
tial gamma-rays at energies above 100 MeV. The source positions are plotted in Galactic 
coordinates with color-coded classifications. The sources are represented as point-like ones. 

• One probable radio galaxy detection - Cen A - (there are two other candidates 
at the moment). 

• 66 high- confidence identifications of blazars (BL Lac objects, flat-spectrum 
radio quasars, or unidentified flat-spectrum radio sources). 

• 27 lower-confidence potential blazar identifications. 

• ~ 170 sources not yet identified firmly with known ob- 
jects. 

Fig 11.51 shows the positions in the sky of these sources. 

The detections included have a significance of more than 5<r at less than 10 
degrees from the galactic plane, where the background radiation is stronger, and 
more than 4a otherwise. Most of the sources can be fitted by a single power-law 
spectra, and the corresponding spectral photon index usually ranges between 1.5 
and 3.0. 

Our interest here will be focused on the unidentified sources. It has been shown 
in Romero et al. (1999) that the distribution of the unidentified EGRET sources 
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Figure 1.6: Distribution of EGRET sources with galactic latitude on the left panel and 
with galactic longitude on the right panel (from Romero et al. 1999). 



with galactic coordinates has a clear concentration of sources on the galactic plane 
plus a concentration in the general direction of the galactic center (see Fig. Il.fi)! . 
This indicates a significant contribution from galactic sources. A correlation ana- 
lysis between 3EG-low-latitude sources and bright and giant HII regions, the usual 
tracers for the galactic spiral structure, shows that there is a strong correlation at 
~ 7<7-level with the spiral arms of the Galaxy, the places where stars are formed 
(Romero 2001). This means that there is a significant number of Population I ob- 
jects 6 in the parent population of the low-latitude gamma-ray sources. 

The large error boxes in the location of the EGRET sources (typically ~ 1 degree 
in diameter), make correlation studies to find secure counterparts very uncertain. 
This technique has been successful mainly to identify sources at high latitudes (such 
as bright blazars), where the turmoil is significantly lower than in the galactic plane. 
Complementary tools like population studies using the known characteristics of the 
gamma-ray sources, like their spectra, variability and flux density had to be used 
also. The so-called log N — log 5* 7 studies can be particularly useful on this respect 
(e.g. Reimer & Thompson 2001, Bhattacharya et al. 2003, Bosch- Ramon at al. 
2005). Combining this with studies including variability and source distribution 

6 A term used to describe stars and other objects, such as star clusters, that tend to be found 
in, or near to, the plane of a spiral galaxy and follow roughly circular orbits around the center. 
They are younger than Population II objects, have a relatively high heavy element content, and 
have probably been formed continuously throughout the lifetime of the disk. Extreme Population 
I objects are found in spiral arms and consist of young objects, such as T Tauri stars, O stars, 
B stars, and stars newly arrived on the zero-age main sequence. Older Population I objects 
include stars like the Sun. All Population I stars are relatively rich in elements heavier than 
hydrogen and helium since they formed from clouds of gas and dust which contained the products 
of nucleosynthesis from previous generations of stars. 

7 A \ogN — log S plot displays the number of sources (N) with a flux (S) greater than a given 
value. 
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models that take into account the non-uniform detection sensitivity across the sky 
(Gehrels et al. 2000; Grenier 2 001, 2004; Grenier et al. 2004), it has been possible 
to establish a division of the sources into two well-defined groups: 

I) One at low latitudes, |6| < 5°, on the galactic plane and spiral arms. These 
are bright sources, relatively hard (photon index V ~ 2.18). This group should 
be formed by young sources (a few million years at most) with isotropic lu- 
minosities in the range of lfj 34-36 erg/ s. These sources contain a subgroup of 
clearly variable sources (Torres et al. 2001a, Nolan et al. 2003). 

II) The other group, at mid-high latitudes, includes 3 populations: 

• 45 ± 6 sources spatially associated with the Gould belt (star forming 
region at ~ 600 pc), with V ~ 2.25, stable and weaker. In this group 
young sources are also expected but with luminosities in the range 10 32 ~ 33 
erg/s. 

• 45 ± 9 sources with softer spectra (r ~ 2.5) and high variability forming 
a kind of halo around the galactic center with a scale height of ~ 2 kpc. 
The sources in this group might be formed in and ejected from globu- 
lar clusters (however, there is no significant correlation with individual 
clusters) or from the galactic plane. These sources should be old, age 
measured in Gyrs, and very luminous, in the range 10 35-37 erg/s. Signif- 
icant variability between different EGRET viewing periods is observed 
in many of these sources (Nolan et al. 2003). 

• Around 35 potentially extragalactic sources (isotropically distributed). 
The transition between both groups occurs at |6| ~ 5 degrees. 



1.3 Motivation of the Thesis: the Link 

A few words now about the link we will propose between MQs and gamma-ray 
emission. In the year 2000, Paredes et al. reported the discovery of a new micro- 
quasar, LS 5039 8 (see Fig. II. 7|) . This discovery involved a very interesting issue: 
the microquasar LS 5039 is inside the error box of the unidentified EGRET source 
3EG J1824-1514 (see Fig. II. 9J) . and was therefore proposed to be its counterpart 
(Paredes et al. 2000). 

Two similar cases followed. The discovery of radio jets in LS I +61 303 confirmed 
the microquasar nature of this binary system (Massi et al. 2001, 2004). There is 
an argument that seems to suggest that LS I +61 303 is also a gamma-ray emitting 
microquasar (Massi 2004a): the likely correlation of the short term variability of 

8 The presence of an asymmetric two-sided jet in LS 5039 was confirmed with observations from 
EVN and MERLIN (Paredes et al. 2002a). The longest jet arm was estimated to reach up to ~ 
1000 AU (see Pig EE). 
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Figure 1.7: VLB A observation where the relativistic radio jets of the microquasar LS 
5039 can be seen (from Paredes et al. 2000). 

3EG J0241+6103 with the orbital period of the binary. The other case has just been 
proposed by Combi et al. (2004) through the association between the microquasar 
candidate AX J1639.0-4642 and 3EG J1639-4702. This latter case is the most un- 
certain one, since there are other potential gamma-ray emitting objects inside the 
error box and the microquasar nature of AX J1639. 0-4642 is not demonstrated. 

The original proposal by Paredes et al. (2000) motivated the research project 
described in this thesis. We investigated the possibility that MQs could be the par- 
ent population of the gamma-ray sources distributed on the galactic plane and at 
high latitude forming a halo around the galactic center and whose spectral energy 
distribution peaks in the MeV-GeV range. 

Though our first aim was to propose a new parent population for different groups 
of EGRET sources, the model of microquasars emitting gamma-rays will also be 
used to explain observations from well-known sources at other wavelengths than 
EGRET's. The models described in this thesis can also be applied outside this 
energy range, as it will be shown for different cases. 



1.4 Outline of the Thesis 

In the next chapter we will review the main features of each component of a micro- 
quasar. The compact object, BH or NS, the accretion disk, the corona, the stellar 
companion and the jet: main characteristics and some models will be presented. 
We will focus on those aspects that are more relevant from the point of view of the 
possible high-energy emission. 
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Figure 1.8: Observations of LS 5939 made with EVN (left panel) and MERLIN (right 
panel) confirming the existence of a jet in that object. Axis units are in milliarcseconds 
(from Paredes et al. 2002a). 



After the second chapter, we will start the presentation of our proposal of MQs 
generating gamma-ray emission. Different scenarios and combinations of parame- 
ters will be studied. In Chapter 3 we will introduce a high-mass microquasar model, 
showing the fundamental role of a high mass companion. We will develop a leptonic 
model studying the interaction of the jet with the field of photons from the stellar 
companion (Inverse Compton, IC). The variability observed in some sources will 
also be considered. 

The toy model of Chapter 3 will become more complex in Chapter 4 where 
the interaction of the jet with the accretion disk and the corona photon fields will 
be taken into account in order to analyze the possible contributions to the total 
gamma-ray production. We will also discuss the application of the model as an 
attempt to interpret some repeated soft gamma-ray events observed in Cyg X-l. 

In Chapter 5 we present another kind of model to produce high-energy gamma- 
rays in HMMQs. It is based on hadronic interactions between multi-TeV protons of 
the jet and the ions present in the wind of the stellar companion. The gamma-ray 
emission results from the decay of neutral pions that are produced in the inelas- 
tic collisions between the protons from the jet and from the wind. This kind of 
scattering implies also the generation of neutrinos. We include estimations for this 
emission as well. 

Low- mass microquasars are finally considered in Chapter 6. We come back to 
a leptonic model and we propose LMMQs as the parent population of one of the 
groups of the unidentified mid-high latitudes EGRET sources. We focus on the 
sources that are forming a kind of halo around the galactic center. In this scenario 
it will be fundamental to consider not only the interactions between leptons of the 
jet and photons from external sources, but also the IC interactions with the syn- 
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Figure 1.9: Location map of the EGRET source 3EG J1824-1514- We can see the 

statistical probability that the gamma-ray source lies within the 50%, 68%, 95% or 99% 
contour. LS 5039, here represented by a filled circle inside an open cirle, has coordinates, 
I = 16. 88°, and b = —1.29°, and lies in the 95% contour. It's the only X-ray and radio 
source in the EGRET error box (from Ribo 2002). 



chrotron photons (Synchrotron-self-Compton, SSC). 

Four appendices complete the thesis. The first one is the presentation of an- 
other gamma-ray producing situation, also studied during the thesis working pe- 
riod. It does not involve MQs but a specific unidentified gamma-ray source, 3EG 
J0542+2610 that is proposed as the gamma-ray counterpart of the Be/X-ray tran- 
sient A0535+26. This is the case of gamma-rays from magnetized accreting pulsars, 
where the high-energy photons are produced in the interaction of a beam of rela- 
tivistic protons and a transient accretion disk. 

Relevant radiative processes, that are used in the models presented in the former 
chapters of this work, are described in Appendix |Bj A detailed explanation about 
superluminal motions is presented in Appendix O The last appendix lists the 
publications completed during the development of the thesis. 



Microquasars: Main Features 
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In the following sections we will review the main statistical and physical charac- 
teristics of MQs. MQs were discovered in recent times, and the number of objects 
recognized as such is still far from representing a statistically important figure. 
Questions like "Can they be assumed as rather common objects in our Galaxy?", 
are still lacking a clear answer. 

On the other hand, we have to deal with the observational data we already 
have in order to start to understand which are the main features of the different 
components of MQs and which are the active physical processes that are taking 
place in them. Observations made at different wavelengths are being interpreted 
and different models are being proposed. 

We need to ponder how significant are the chances to rise the number of con- 
firmed MQs and also to understand the physics involved, if we want to be able to 
propose them as a parent population of EGRET unidentified sources. 



2.1 Statistics and Projections 

Taking into account that MQs are until now considered as a particular subgroup 
of XRBs that generate jets, the first step is to wonder how many XRBs are hosted 
in the Galaxy. The most updated catalogs of XRBs are found in Liu et al. (2000) 
for High Mass X-ray Binaries (HMXBs) and Liu et al. (2001) for Low Mass X- 
ray Binaries (LMXBs). After an increase from ~ 30 known HMXBs in 1983 (van 
Paradijs 1983) to 69 in 1995 (van Paradijs 1995), the new catalog contains 130 of 
these sources. The number of registered LMXBs went from 33 in 1983 (van Paradijs 
1983, Bradt & McClintock 1983) to 119 in 1995 (van Paradijs 1995) to reach in the 
newest catalog, 150 confirmed of such XRBs. The two catalogs have then 280 XRBs 
in total (several new binaries have been discovered since then). 

Next step is to check how many of those 280 XRBs are also radio emitters 
(Radio Emitting X-ray Binaries, REXBs), since this is the wavelength at which 
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jets manifest themselves more clearly. Taking into account both catalogs, there 
are ~ 75 X-ray pulsars, which are suspected not to emit in radio. It is believed 
that the strong magnetic field of X-ray pulsars disrupts the accretion disk at some 
thousand kilometers from the neutron star, suppressing the possibility to generate 
a jet (Fender & Hendry 2000). In fact, synchrotron radio emission has never been 
detected from any of these sources. Otherwise, the two catalogs contain 43 radio 
emitting sources (Paredes 2004b). Once detected, the radio emission has to be re- 
solved in order to confirm that it is produced by a relativistic jet. Only when the 
presence of a jet is unquestionable we can call the XRB a microquasar. 

Nowadays there are 16 systems identified as MQs. They are listed in Ta- 
ble |2H](fro m Paredes 2004a and 2004b with some modifications) 1 . Different multi- 
wavelength cross-searchings have been carried out in X-ray and radio source catalogs 
in order to find more MQ candidates that could rise this number after being con- 
firmed as such with VLBA, ATCA or VLBI observations (e.g Paredes et al. 2002b). 
The chances to increase the number of MQs radically would enhance if we could 
know that there are even more XRBs in the Galaxy. Grimm et al. (2002) have 
made some estimate in this respect. Using data of XRBs from RXTE/ASM 2 they 
obtained a log N — log S diagram and the luminosity function. The total numbe 
r of XRBs brighter than 2 x 10 35 erg s" 1 is ~ 190 divided into ~ 55 high mass 
and ~ 135 low mass binaries. Extrapolating the luminosity functions towards low 
luminosities they estimated a total number of ~ 705 XRBs brighter than 10 34 erg 
s _1 , distributed as ~ 325 LMXBs and ~ 380 HMXBs. As a consequence of these 
expectations the existence of MQs in a larger number becomes more likely. In fact, 
some authors sustain that jets are a quite common feature of X-ray binary systems, 
suggesting that up to 70% of XRBs might produce them (Fender 2004). In such 
a context it would be a matter of time to detect them, providing that new instru- 
ments with more sensitivity and resolution in both X-rays and radio-band become 
available. We will come back over the arguments for such a high percentage in 
Section ESI 



2.2 Physical Components 

In the next four subsections we will make an overview of the main components of a 
microquasar, studying their physical properties, the classifications that arise from 
them and some of the models that attempt to describe the observational data. These 
components are: the companion star, the disk, the corona, and the jet. Fig. 12.11 
sketches a microquasar system showing the range of wavelengths corresponding to 
each emitting component. 



1 XTE J1118+480 is included in the list because even though its relativistic jet has not been 
resolved yet it is inferred from theoretical reasons. 

2 ASM, All-Sky Monitor, instrument on board of the satellite RXTE, Rossi X-Ray Timing 
Explorer, also known as XTE. 
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Notes: ( a ) NS: neutron star; BH: black hole. ^ p: persistent; t: transient. ( c ) jet inclination, 
(d) p rec: precession; XRJ: X-ray jet. *Reported by Corbel et al. 2002. ^Recently reported by Gallo et al. 2004. 
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Figure 2.1: Sketch of the physical components of a microquasar. We can see the binary 
system formed by a compact object (NS or BH) that accretes material from the stellar 
companion, forming an accretion disk and corona. The potential accretion energy is 
converted into high energy emission and mechanical energy of the jet, that emits from 
radio wavelengths to gamma-rays, as we will propose in the forthcoming chapters (from 
Fender & Maccarone 2004)- 



2.2.1 The Companion Star and the Compact Object 

According to the mass of the stellar companion, XRBs are divided in two groups: 
high mass and low mass X-ray binaries. HMXBs are systems where the donor star 
is a young massive star of ~ 8 — 20 M and orbital periods are several days. The 
nature of the star determines a further division of HMXBs in two subgroups: those 
in which the stellar companion is a Be star (Be/X-ray binary) and those in which it 
is an O or B supergiant (SG/X-ray binary). In general, these high mass companion 
stars transfer mass to the compact object through strong stellar winds, though in 
some cases it can also happen via Roche Lobe overflow. X-ray outbursts are ex- 
pected to happen in the Be/X-ray sources during the periastron passage. 

In the case of LMXBs, the companion is an older low-mass type of star, with 
M < 2 M ; it can be a white dwarf, a late-type main-sequence star, an A- type 
star or an F-G sub-giant (in which case, it may be the remnant of a star that 
originally was of intermediate mass ~ 1.5 to 4 M ). The orbital periods are in the 
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Figure 2.2: All-sky map where open circles representing 86 LMXBs and filled circles, 52 
HMXBs. The two spatial distributions in the Galaxy can be noted: HMXBs are signifi- 
cantly concentrated on the galactic plane whereas LMXBs are mainly forming a halo in 
the galactic bulge (from Grimm et al. 2002). 

range of 0.2 to 400 hours and the mass transfer occurs through Roche Lobe overflow. 

The already mentioned ASM catalog includes 340 sources being 217 galactic, 
112 extragalactic, and 10 unidentified. Studying the galactic sources, Grimm et 
al. (2002) found significant differences in the 3D spatial distribution of HMXBs 
and LMXBs in our Galaxy as it is shown in Fig. 12.21 (in good agreement with 
theoretical expectations and earlier results - van Paradijs & White 1995; White 
& van Paradijs 1996; Koyama et al. 1990; Nagase 1989). Whereas HXMBs are 
more concentrated towards the galactic plane with a vertical scale height of 150 
pc, and clear indications of a distribution following the spiral structure, LMXBs 
have a strong tendency to concentrate towards the galactic bulge and their vertical 
distribution has a scale height of 410 pc. This important difference in the angular 
distribution of high and low mass XRBs is also illustrated by Fig. 12.31 where the 
angular distribution is represented against the galactic latitude and longitude. 



The compact object that is present in an X-ray binary system can be either a 
neutron star or a black hole. While stars are on the main sequence, the thermal 
pressure of hot gas maintains an equilibrium with gravity. The source of energy 
that provides the pressure comes from nuclear reactions generated in their cores. 
The life of the star evolves out from the main sequence through the so-called giant 



26 



2 MICROQUASARS: MAIN FEATURES 




-40 -20 20 40 -100 100 

gal. latitude b" [deg.] gal. longitude l 11 [deg.] 



Figure 2.3: Angular distribution of galactic HMXBs (solid lines) and LMXBs (thick 
green lines) against the galactic latitude (left panel) and longitude (right panel). These 
two graphs illustrate the well-known fact that HMXBs are strongly concentrated towards 
the galactic plane. An important difference in the longitude distributions of HMXBs 
and LMXBs can be noticed, with the second ones significantly concentrated towards the 
galactic center/bulge and the former distributed in clumps approximately coinciding with 
the location of tangential points of the spiral arms ( whose position is marked by arrows 
in the right panel). The LMXBs number is divided by 3 on the right panel (from Grimm 
et al. 2002). 



branch towards the final phases when the outer layers of the star are ejected. During 
all that path the nuclear processing occurs using all the available nuclear energy re- 
sources of the star. In the most massive stars, M > 10M Q , it is quite probable that 
the nuclear burning continues till the formation of iron. At the end of this phase 
when the core of the star runs out of nuclear fuel, gravity overtakes and the star col- 
lapses until another form of pressure appears to establish a new kind of equilibrium. 
The new possible configurations in this new stag e after the star collapses are are 
known as white dwarfs, neutron stars and black holes. These are called "dead stars" . 

Statistics made with observed cases seems to indicate that low mass stars, 1 < 
M < 4M Q , evolve into white dwarfs, whereas more massive stars, 4 < M < 1OM , 
towards NSs. Many of these latter cases are formed in supernova explosions. Super- 
nova explosions associated with even more massive stars, M > 10M Q , are assumed 
to end up in the formation of a stellar-mass BHs. 

As we stated, the compact objects that are of our interest are NSs and BHs. In 
NSs, the pressure that provides the force opposed to gravity is the quantum mechan- 
ical pressure associated with the fact that neutrons are fermions and therefore only 
one particle is allowed to occupy one quantum state at a time. NSs are then held 
up against gravitational collapse by neutron degeneracy pressure and have masses 
in the range of ~ 1.5M to 3M & with a diameter that goes from 10 to 20 km. The 
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corresponding escape velocity is about half the speed of light and hence the path of 
light is strongly curved by gravity. Theories support the concept that NSs should 
have a solid crust and a liquid core with properties of superfluidity, having almost 
no viscosity. 

NSs are the last stable-star form. In the case of more compact objects the attrac- 
tive force of gravity can not be balanced by any known physical force and the last 
collapse to a physical singularity is unavoidable. Even electromagnetic radiation 
cannot escape from this singularity since the escape velocity exceeds the speed of 
light, hence the origin of the name black holes. The escape velocity from the surface 
of a star of mass M and radius r can be obtained through classical calculations: 
v = (2GM/r) 1 ^ 2 . If we set v as the speed of light (c) the radius of such a star would 
be r = 2GM/c 2 . This is the so-called Schwarzschild radius, rs c h, of a black hole 
of mass M. For non-rotating BHs there is a last stable circular orbit around this 
point mass at 3rs c h- At radii less than this value the particles start spiraling into 
the singularity, located at r = 0. 

Setting the values of the constant in the expression of the Schwarzschild radius 
one obtains a very practical expression: rg c h = 3(M/M ) km. This straightfor- 
wardly implies that for the Sun, rg c h = 3 km, which is negligible compared to the 
solar radius, R Q ~ 7 x 10 5 km. Instead, in the case of a NS with M ~ 3M , we 
find r Sc h ~ 10 km which, being of the order of the neutron star radius, means that 
general relativity becomes important in determining the stability of this case. 

To distinguish whether a binary system harbors a neutron star or a black hole 
as the compact object is not a trivial task since the gravitational field near a NS is 
nearly as strong as near a BH. Differences in the X-rays emitted by inf ailing matter 
in the two cases are expected to be subtle. Though regular pulsations are a sure 
proof that an X-ray emitter is a neutron star 3 , they are not always present. Many 
sources like Sco X-l, are suspected to hold a NS despite the fact that they show no 
sign of regular pulsations. 

On the other hand there are severe differences between NSs and BHs that should 
imply certain differences in the radiation spectra involving one or the other. NSs 
have real surfaces whereas BHs have horizons, which represent the complete absence 
of a physical surface. The characteristic magnetic field present in NSs would also 
be expected to mark a distinction. Nonetheless, up to now, no one has succeeded 
in predicting the spectral differences only from theoretical considerations. 

There is however a decisive discriminant that can at least ensure the absence of 
a NS from some XRB systems: if the present compact object has a mass greater 
than 3M , then one could assume the presence of a BH. Such a mass estimation can 
be obtained by measuring the orbital period and the parameters of the normal star 
in the binary system. In this way the X-ray source in the binary system Cygnus 

3 See e.g. Begelman & Rees (1998) or Longair (1997), for a detailed explanation on the physics 
of pulsars. 
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X-l was inferred to have a mass of at least 6M , becoming in this way the first 
example of a black hole candidate. 

By definition microquasars have relativistic jets. Bright radio detections asso- 
ciated with the production of a relativistic jet are by now a common property of 
accreting black holes whereas they are not so commonly associated with accreting 
neutron stars. Nonetheless the most relativistic flow observed so far in the Galaxy 
has been very recently discovered by Fender et al. (2004a) in the source Circinus X- 
1, whose compact object is a neutron star (Tennant et al. 1986). It was established 
that the apparent superluminal velocity is /3 app > 15. Knowing that for an intrinsic 
velocity (3 at an angle 9 to the line of sight, /3 app has a maximum value of /3 app = Yf3 
when (3 = cos 8 (see Appendix EJ) , then Tf3 > 15, and it follows that the minimum 
velocity solution is T > 15, with (3 > 0.998. In order to have a comparison, let 
us mention MQs with outflowing components are estimated to have bulk Lorentz 
factors in the range of 2 < T < 5; in the case of active galactic nuclei (AGN) typical 
Lorentz factors are T ~ 10. 

An important comment is that the lower limit to T may correspond to a pertur- 
bation (shock) propagating along the flow instead of the underlying jet. We should 
then consider the relation Tbuik = Tshock/ \/2~, that is valid in the approximation of a 
lD-relativistic shock propagating along a relativistic fluid, in the limit /3 s hock ~~ > !■ 
This would imply Tbuik > 10. 

It seems that Sco X-l is another neutron star binary that may be generating 
ultrarelativistic jets (Fomalont et al. 2001). The discovery of ultrarelativistic flows 
in XRBs with a NS instead of a BH as the compact object, leads to the conclusion 
that the generation of such jets is not related to any property that would be ex- 
clusively of BHs, like the event horizon or the ergosphere. It rather highlights the 
role of the accretion flow (which is a common feature of both NSs and BHs) as the 
main element that might lead to the formation of the jet. 



2.2.2 The Accretion Disk 

We will devote this section to the physical description of the accretion disk that is 
present in MQs as it can be inferred from X-ray observations. The accretion disk 
component is usually modeled by the so-called "Standard Disk Model" in which 
the accreting gas forms a geometrically thin and optically thick disk, producing a 
quasi-blackbody spectrum due to thermal emission. The effective temperature of 
the accreting gas is in the range 10 5 — 10 7 K, depending on the compact object mass 
and the accretion rate, T e g oc M _1//4 M 1//4 . 

Accretion processes around compact objects imply the presence of rotating gas 
flows. The situation is therefore modeled by hydrodynamic equations of viscous 
differentially-rotating flows. The "standard disk model" is the name of the most 
famous solution to this set of equations. The model was developed by Shakura & 
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Sunyaev (1973) and Novikov & Thorne (1973). 

If a particle is in a circular orbit around a central gravitating body, it will stay in 
that orbit. If then energy and angular momentum are extracted from the particle, 
it will spiral slowly inwards. In the context of the standard disk model, the viscosity 
will cause the gas to lose energy and will have the effect of transporting angular 
momentum outwards, allowing the accreting gas to spiral in towards the central 
mass. The amount of energy that can be extracted by such a process is equal to the 
binding energy of the innermost stable orbit. Around a neutron star, ~ 10% of the 
rest mass can be extracted and up to 40% for orbits around a black hole (Pringle 
1981). Viscosity also acts as a source of heat; some or all of this heat is radiated, 
leading to the observed spectrum. The accretion process can then be an efficient 
converter of rest mass to radiation. 

It is known already from the 1920s (Jeffreys 1924) that the action of viscosity 
on an initial ring is to spread it out. Most of the mass moves inwards losing energy 
and angular momentum, but a tail of matter moves out to larger radii in order to 
conserve the angular momentum. 

Let us develop this idea in a more quantitative way 4 . A thin disk means that 
the matter lies very close to the plane, which in cylindrical coordinates (R, 0, z) 
means z = 0. We assume that the gas moves with angular velocity Q in circles 
around the accreting compact object which has mass M and radius i?* 5 . When the 
angular velocity has the Keplerian value then, 



Since the gas is being accreted, its radial velocity, v^(R,t), near the compact ob- 
ject, is negative. Another important magnitude is the disk surface density, S(i2, t), 
which is the mass per unit surface area of the disk, obtained by integrating the gas 
density, p, in the z-direction. 

An annulus of the disk material lying between R and R + AR has a total mass 
2ttRARI1 and total angular momentum 2ttRARHR 2 Q. Taking into account that 
the rate of change of both quantities is given by the net flow from the neighboring 
annuli and making the limit AR — > 0, we get the mass conservation equation: 

4 Skipped steps in any of the calculations included in this section can be followed in Shapiro & 
Teukolsky (1983) and Frank et al. (2002). 

5 For the case where the compact object is a BH, i?* should be replaced by i?i n in all the 
following calculations, with R m the innermost stable orbit. 




(2.1) 



which means that the circular velocity is, 



= RQ K (R). 



(2.2) 
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To obtain the conservation equation for angular momentum the procedure is the 
same except that we have to include the transport due to the net effects of the 
viscous torques, G(R,t): 



where 



G(R,t) = 2nRui:R 2 n', (2.5) 



with fi' = dfi/dr, and v, the kinematic viscosity. 

Assuming dVt/dt = 0, which holds for orbits in a fixed gravitational potential, 
performing some algebraic combinations with all these formulas, we obtain: 



^ _ 3 d f „, „ d 
~dt ~ RdR 



{R 1/2 ^R 1/2 ]}, (2-6) 



that is the basic equation governing the time evolution of surface density in a Ke- 
plerian disk. In general it is a nonlinear diffusion equation of E, because v may be 
a function of local conditions in the disk, like S, R and t. Therefore, in order to 
solve ([2.6)1 . we need some prescription for v. Unfortunately, the nature of viscosity 
is still quite unclear, and as a consequence, its mathematical form is still uncertain 
and not unique. We gain some insight into the evolution of disks by choosing forms 
for v that makes (|2.6|) solvable, which is just a mathematical help to show how the 
disk's dynamics works, rather than a physically based prediction of the theory. For 
instance, assuming v = constant, ()2.6|) can be solved by separation of variables as: 

S(x,r) = r- 1 x- 1 /4 e xp|-ii±^}/ 1/4 (2x/r), (2.7) 



where / 1 / 4 (2x/r) is a modified Bessel function. Taking as the initial matter distri- 
bution a ring of mass m at R — Ro, x = R/Rq and r = YlvtR^ 2 are dimensionless 
radius and time variables. Fig 12.41 shows £( function of x for different 

values of r. We can see in the figure that viscosity has the effect of spreading the 
original ring in radius. 



The asymptotic behavior of 7 1 /4(2x/r), 
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Figure 2.4: A matter ring of mass m placed in a Keplerian orbit at R = Rq spreads out 
under the action of viscous torques. Surface density, X, vs x = R/Rq, the dimension- 
less radius, for different values of the dimensionless time, r = YlvtR^ 2 , where v is the 
viscosity. 
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2x 



> 1 



OC 



2x\ 1/4 2x 



< 1 



(2.8) 



implies that 



3u ( 1 2x 2 . , 

V R ~ ~5- F + ^ > for 2x > r, 

it [4:X T T 



(2.9) 



and 



3v f 1 2x 



A, 2i r 



< for 2x < r. 



(2.10) 



Thus, we confirm that the outer parts of the matter distribution (2x ^> r) move out- 
wards, taking away the angular momentum of the inner parts, which move inwards 
towards the accreting compact object. 
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In the steady state 6 the disk structure is determined by solving simultaneously 
four conservation equations: conservation of mass, angular momentum, energy, and 
vertical momentum (where an expression for the pressure is needed). In addition a 
prescription for v must be specified as well as a law describing the vertical energy 
transport from the center to the surface (where an expression for the opacity is 
needed). The set of equations is the following: 

1. REST- MASS CONSERVATION 



M = 2-kRLvh = constant. 



(2.11) 



2. ANGULAR MOMENTUM CONSERVATION 



3tt 



1/2 



(2.12) 



3. ENERGY CONSERVATION 

The kinematic viscosity, u, generates dissipation in the disk at a rate D(R) 
per unit area per unit time, such that: 



D(R) 



3GMM 
8ttR 3 



1- — 



R* 
R 



1/2 



(2.13) 



An interesting issue is that expression ()2.13|) is independent of the major un- 
certainty of the disk's theory, the viscosity. This is at the expense of the 
assumption that the chosen viscosity prescription can be adjusted to provide 
the steady mass flux M. 

Using (|2.13|) we obtain the total disk luminosity: 



(2.14) 



-disk = / D(R) 2nR dR 
1 GMM 

2 R* 

6 We impose the steady state condition by setting d/dt = in 1)2.3)1 and 1)2. 4[l . obtaining in 
this way 1)2.11)1 and 1)2.12(1 . The steady assumption will be valid whenever the external conditions 
change on timescales longer than the time t V i SC ~ =^-, in which changes in the radial structure of 
the disk are registered, i.e. the timescale of the spreading because of viscosity effects. 
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This is only one half of the total available accretion energy, since the total 
potential drop from infinity to is GM/R*. The point is that the matter 
just outside the boundary layer still retains as kinetic energy one half of the 
potential energy it has lost in spiraling in. Hence, the rest of the accretion 
luminosity is emitted in the boundary layer. This has a practical implication: 
the study of the details of emission from the inner edge of an accretion disk 
can be just as important as studying the emission from the disk itself. 

4. VERTICAL MOMENTUM CONSERVATION 

Since there is no net motion of the gas in the vertical direction, momentum 
conservation along the ^-direction reduces to a hydrostatic equilibrium condi- 
tion. Equating the component of the gravitational force of the compact object 
along the ^-direction to the vertical pressure gradient in the disk, 



ldP GM z , v . 



Replacing the differentials by finite differences, setting AP m P and z = H, 
yields: 

R 3/2 

fl =W (2 ' 16) 



where c s is the sound speed, (c^ = P/p), P, pressure, and p, density of the 
disk. 

5. PRESSURE 

The total pressure of the disk material is the sum of gas and radiation pres- 
sures: 

p = e m + p t (217) 



where k is the Boltzmann constant, p is the mean molecular weight, m p the 
proton mass, a the Stefan-Boltzmann constant, and T c the central tempera- 
ture: T(R,z) = T(R,0). 

6. RADIATIVE TRANSPORT 

The temperature T c can be given by an energy equation relating the energy 
flux in the vertical direction to the rate of generation of thermal energy by 
viscous dissipation. The heat generated internally by viscosity dissipation is 
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transported vertically through the disk before being radiated at the surface. 
The vertical energy transport mechanism can be either radiative or connec- 
tive. We will here assume the radiative one (see Frank et al. 2002, Section 
5.8, for a justification of this choice). 



We then equate the heat input per unit area, D(R), with the heat loss per 
unit area given by the radiative transport, 



4a 
3t" 



T 



D(R) 
3GMM 



~R 



1/2" 



(2.18) 



where r is the optical depth that is defined through 

T = k kP H = k r S, (2.19) 



with kr = «r(p, T c ), the total Rosseland mean opacity. The opacity in the 
disk is related with the photon absorption (mainly due to "free-free" transi- 
tions) and the photon scattering (Thomson scattering in this case). 

The optically thick condition of the standard disk model implies that r 1. 
If r becomes < 1 the radiation escape directly; the radiative transport is then 
not required and equation ()2.18|) is no longer valid. 

7. VISCOSITY 

The standard model uses for the viscosity the so-called a-prescription of 
Shakura & Sunyaev (1973). Suppose that the dominant process for redis- 
tributing angular momentum is a turbulent viscosity. The effective kinematic 
viscosity of a turbulent process is given by v ~ Iv, where I the size and v 
the velocity of the largest eddies in the flow. In an accretion disk, we can 
assume that the scale of the eddies is less than the disk thickness, H, and the 
turbulence is subsonic 8 . Consequently, 

u = ac s H, (2.20) 



with a < 1. It is important to realize that ()2.20|) is a parametrization that 
helps to isolate in a the present lack of knowledge on the viscosity. We have 

7 Since we are working under the geometrically thin disk condition (H <C R), the disk medium is 
then essentially "plane-parallel" , so that the temperature gradient is effectively in the z-direction, 
and we need only to consider the energy transport perpendicular to the disk. 

8 In the case that the turbulence would be supersonic, the turbulent motions would probably 
be thermalized by shocks. 
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at least obtained a constraint over the new parameter, a < 1. Though just a 
parametrization, the a-prescription has been very useful and has promoted a 
semi-empirical approach to the viscosity issue, which looks for estimating the 
magnitude of a by comparison of theory and observation. 

An alternative to the turbulent viscosity process and in fact one of the more 
plausible, is the magnetic viscosity (Lynden-Bell 1969, Shakura & Sunyaev 
1973), where the transfer of angular momentum occurs through magnetic 
stresses due to a field anchored in the disk. This kind of viscosity can also be 
contained within the a-prescription, with a ~ v\/c^, where v& is the Alfven 
speed in the disk. 

A more radical approach to the idea of using magnetic fields to transfer away 
angular momentum has been proposed by Blandford (1976) and by Lovelace 
(1976). The energy and angular momentum of the disk can be carried away 
in the form of a magnetized relativistic wind (See also Blandford & Znajek 
1977; Blandford & Payne 1982). 

Solving the above listed set of equations we obtain the magnitudes p, E, H, vr, 
P, T c , t, and v as functions of M, M, R and a. It can be shown that for fixed 
values of M and M, the disk can be divided into three distinct regions, depending 



• An outer region, at large R, in which gas pressure dominates over radiation 
pressure and in which the opacity is mainly due to free-free absorption. 

• A middle region in which gas pressure still dominates over radiation pressure 
but the opacity is mainly due to electron scattering. 

• An inner region, at very small R, in which radiation pressure dominates over 
gas pressure and again the opacity is dominated by electron scattering. 



Let us now concentrate on the emitted disk spectrum. A very important conse- 
quence of the assumption that the disk is optically thick is that each element of the 
disk radiates roughly as a blackbody with a temperature T(R). This temperature 
is given by equating the dissipation rate D(R) per unit area to the blackbody flux: 



on R: 



<jT\R) = D(R). 




Using (EHSj) for D(R) we obtain, 
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Figure 2.5: The integrated spectrum of a steady accretion disk that radiates a local black- 
body spectrum at each point. The units are arbitrary, but the frequencies corresponding to 
T ou \, = T(R ont ) and are marked (from Pringle 1981). 



87ri? 3 



1/2' 



1/4 



(2.22) 



This temperature is analogous to the effective temperature of a star. We approxi- 
mate the emitted spectrum from each element of area of the disk as: 



2hi> 3 

I v = B V [T{R)] = c2(eWfcT(fl) _ 1) (erg s" 1 cm" 2 Hz" 1 sr" 1 ). (2.23) 



For an observer at a distance D whose line of sight makes an angle 9 to the normal 
to the disk plane (8 is the binary inclination), the flux at frequency v from the disk 
is 

27TCOS0 f^ ut 



S„ = D2 I I v R dR, (2.24) 



where R out is the outer radius of the disk. With the blackbody assumption, we get 



4nhu 3 cos9 f Rout R dR 
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Notice that the expression we obtain for S v is independent of the disk viscosity. This 
is a consequence of the assumptions of steady and optically thick disk. This implies 
that if there is any important difference between the spectrum obtained with (J2.25|) 
and the observed ones, it will not be due to our ignorance of the viscosity but to 
the assumptions made. 

The integrated spectrum given by (j2.25|) is shown in Fig. 12.51 and its functional 
form is the following: 



The integrated spectrum is therefore a stretched-out blackbody. The flat part 
corresponding to S v oc v 1 ^ 3 characterizes the disk spectrum: unless T out = T(R out ) 
is significantly smaller than T* this part of the curve may be quite short and the 
spectrum is not very different from a blackbody. 

Finally, some remarks with respect to how the nature of the compact accreting 
object may affect the disk emission. If it is a star with a strong magnetic field, we 
already mentioned in previous section that the disk flow can be disrupted at some 
distance from the compact object itself. In this case most of the accretion lumi- 
nosity comes from the matter when it actually strikes the compact object surface 
and not from the disk. If the accretion disk is around a black hole, the motion 
near the inner regions deviates from the simple Keplerian law. It is assumed that 
the disk extends down to the innermost stable circular orbit R[ n . Within the ra- 
dius -R; n the gas spirals into the hole without radiating (see Stoeger 1980 for a 
justification). The energy released by the accreting matter corresponds then to the 
binding energy at the radius R m . For a Schwarzschild black hole this corresponds 
to about the 6% of the rest mass energy and, for a maximal Kerr hole, to about 42%. 

We end this section mentioning that the major success of the standard disk 
model is that the surface disk temperature it predicts approximately corresponds 
to the observed UV emission in AGNs and soft X-ray emission in galactic black 
holes. 



A strong motivation for much of the work on nonthermal and optically thin emis- 
sion from accretion disks was provided by the observations of the X-ray binary 
source Cyg X-l. There is reasonable evidence that the X-ray source is powered 
by disk accretion onto a black hole, but the observed X-ray spectrum is too hard 
to be produced by the standard, optically thick disk model that we just described 
in the previous section. With that model a maximum disk surface temperature is 




S v oc v 2 

S v OC I/V3 

the spectrum drops exponentially 



2.2.3 The Corona 
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T* ~ 10 5 K for AGNs (M ~ 10 7 - 1O 9 M ) and % ~ 10 7 K for stellar-mass black 
holes (M ~ 1OM ). A hard X-ray component would, however, correspond to a 
temperature ~ 10 9 K. The standard disk model may yield such high temperatures 
only if the accretion rate approaches the critical value, M crit , that implies a high 
radiative efficiency of accretion and the situation then corresponds to the critical 
Eddington luminosity. Near-critical disks c an strongly deviate from the blackbody 
state provided a is large (Shakura & Sunyaev 1973). An overheating then occurs in 
the inner region of the disk where the inflow time-scale is shorter than the time-scale 
for relaxation to thermodynamic equilibrium. 

The region supposed to be responsible for emission that gives rise to the ob- 
served hard X-rays is named "Corona" (see Fig. 12. lj) . Since the observations of Cyg 
X-l, different models appeared, based on different assumptions and physics, but 
having the same aim to interpret and fit those observed spectra that could not be 
fully explained with the standard disk model. We will focus on two of them, that 
are still being developed and tested: the "Disk - Corona" model and the "Advection 
- Dominated Accretion Flow" model or ADAF. 

The ADAF model originates as an alternative solution to the same hydrody- 
namic equations of viscous differentially-rotating flows that were solved in the stan- 
dard disk model. Instead the Disk - Corona model implies to add a new component 
to the standard model. 



The ADAF model is a solution of the hydrodynamic equations of viscous differen- 
tially-rotating flows (listed in Section 12.2. 2|) that works with low, sub-Eddington, 
accretion rates (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994, 1995a, 1995b; 
Abramowicz et al. 1995) 9 . In this solution, the accreting gas has a very low den- 
sity which implies on one hand that it is optically thin and on the other that it is 
unable to cool efficiently within the accretion time. The viscous energy is therefore 
stored in the gas as thermal energy instead of being radiated, and is advected onto 
the central compact object. Since most of the viscously generated energy in the 
ADAF model is stored in the gas as internal energy, rather than being radiated, 
the gas temperature is quite high. This causes the gas to swell, then H ~ R and 
therefore, geometrically, ADAF s are quasi-spherical. In this model the gas adopts 
a two-temperature configuration that will be explained later on. 

There is a compact way to express the just stated definition: 

<f dv = q + ~ 9", (2.26) 

9 ADAF solutions are also possible at very high accretion rates (see e.g. the review by Reynolds 
k Nowak 2003). 
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where q represents the advective transport of energy, q + is the energy generated 
by viscosity per unit volume, and q~ is the radiative cooling per unit volume. 

Depending on the relative magnitudes of the terms in this equation, three 
regimes of accretion may be identified: 

• q + = q~ ^> g adv : this corresponds to a cooling-dominated flow where all the 
energy released by viscous stress is radiated; the amount of energy advected 
is negligible. The thin disk solution corresponds to this regime. 

• g adv ^ q + ^> q~\ this corresponds to an ADAF where almost all the viscous 
energy is stored in the gas and is deposited into the black hole. The amount 
of cooling is negligible compared with the heating. For a given M, an ADAF 
is much less luminous than a cooling-dominated flow. 

• — g adv = q~ ^> q + : this corresponds to a flow where energy generation is negli- 
gible, but the entropy of the inflowing gas is converted to radiation. Examples 
are Bondi accretion, Kelvin-Helmholtz contraction during the formation of a 
star, and cooling flows in galaxy clusters. 

It is important to point out that models based on the two-temperature ADAF 
solution make certain important assumptions. The validity of these assumptions is 
not yet completely proved and is currently under testing. 

1. Equipartition Magnetic Fields: It is assumed that magnetic fields contribute 
a constant fraction (1 — (5) of the total pressure: 

Pm = ^ = (l-/3)pc s 2 , (2.27) 

where p m is the magnetic pressure due to an isotropically tangled magnetic 
field. The assumption of a constant j3 is fairly innocuous since, in general, 
we expect equipartition magnetic fields in most astrophysical plasmas. In 
particular, Balbus & Hawley (1991) have shown that differentially rotating 
disks with weak magnetic fields develop a strong linear MHD instability which 
exponentially increases the field strength to near equipartition values. ADAF 
models assume (5 = 0.5, which means equipartition between gas and tangled 
magnetic pressure. 

2. Thermal Coupling Between Ions and Electrons: ADAF models assume that 
ions and electrons interact only through Coulomb collisions and that there is 
no non-thermal coupling between the two species. In this case the plasma is 
two-temperature, with the ions much hotter than the electrons. 

3. Preferential Heating of Ions: The two-temperature ADAF model assumes 
that most of the turbulent viscous energy goes into the ions (Shapiro et al. 
1976; Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1995b), and that only a 
small fraction 5 <C 1 goes to the electrons. The parameter 5 is generally set to 
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~ 1CT 3 ~ m e /m p , but none of the results depend critically on the actual value 
of 5, as long as it is less than a few percent. There have been several theoretical 
investigations which consider the question of particle heating in ADAFs like 
Bisnovatyi-Kogan & Lovelace (1997), Blackman (1999), Gruzinov (1998), and 
Quataert (1998). 

4. a Viscosity: The viscosity parameter a of Shakura & Sunyaev (1973) is used 
to describe angular momentum transport; a is assumed to be constant, inde- 
pendent of radius. Some authors have proposed that it may vary as a function 
of (H/R). Since ADAFs have H ~ R, no radial dependence is expected, and 
a constant a appears to be a particularly good assumption (Narayan 1996a). 

The optically thin ADAF has an important constraint: it exists as a solution 
only for m 10 less than a critical value m crit (Ichimaru 1977; Rees et al. 1982; 
Narayan & Yi 1995b; Abramowicz et al. 1995), such that for m < m crit , we have 
q + > q~ and a consistent ADAF solution, whereas for m > fn CT \ t no ADAF is 
possible. At low densities Coulomb coupling between protons and electrons is very 
weak and the amount of viscous energy that is transferred to the electrons is very 
small. Coulomb coupling therefore restricts the amount of energy that can be lost 
by radiation. With increasing m, Coulomb coupling becomes more efficient, and 
at a critical density the coupling is so efficient that a large fract ion of the viscous 
energy is transferred to the electrons and is radiated. Above this accretion rate, 
the flow ceases to be an ADAF and becomes a standard cooling-dominated thin 
disk. The critical accretion rate can be estimated by determining the m at which 
the viscous heating, q + , equals the rate of energy transfer from the ions to the elec- 
trons, q ie . Observations suggest that the two-temperature ADAF solution exists up 
to m crit ~ 0.05 - 0.1. This suggests that a ~ 0.2 - 0.3 in ADAFs. 

It is interesting to remember, at this point, that in the case of the standard 
model there is also a critical value, m crit , that implies a deviation from the thin disk 
feature to an emission model like for instance the ADAF model. The confluence of 
the two possible relations between m and m cr i t gives: 

J Thin disk to ADAF m < m crit 
y ADAF to thin disk m > m crit . 

The spectrum from an ADAF around a black hole ranges from radio frequencies 
~ 10 9 Hz to gamma-ray frequencies > 10 23 Hz, and can be divided into two parts 
based on the emitting particles: 

1. The radio to hard X-ray radiation is produced by electrons via synchrotron, 
Bremsstrahlung and IC processes (Mahadevan 1997). Comptonization process 
is over the locally produced synchrotron and Bremsstrahlung radiation and 
dominates the spectrum. 

10 m = — : , MEdd = — —4, with i] c ff, the efficiency of converting matter to radiation. 

M Edd Vcsc 2 
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Figure 2.6: Schematic spectrum of an ADAF around a black hole. S, C, and B refer 
to electron emission by synchrotron radiation, IC scattering, and Bremsstrahlung, respec- 
tively. The solid line corresponds to a low rh, the dashed line to an intermediate rh, and 
the dotted line to a high rh ~ rhcrit- The gamma-ray spectrum is due to the decay of 
neutral pions created in proton-proton collisions (from Narayan, Mahadevan & Quataert 
1998). 



2. The gamma-ray radiation results from the decay of neutral pions created in 
proton-proton collisions (Mahadevan et al. 1997) 11 . 

Fig. 12.61 sketches the spectrum obtained with the ADAF model, showing the 
places where the emission is generated and the radiative processes that causes it for 
each frequency. 

The low energy end of the spectrum, is due to synchrotron cooling by semi- 
relativistic thermal electrons (Mahadevan et al. 1996). The synchrotron emission is 
highly self-absorbed and is very sensitive to the electron temperature (Mahadevan 
1997). The emission at the highest (peak) frequency comes from near the black hole, 
while that at lower frequencies comes from further out. The peak frequency varies 
with the mass of the black hole and the accretion rate, roughly as ^p eak oc m~ 1//2 m 1//2 

11 This remains speculative since gamma-ray emission has never been observed from an accretion 
disk. As we will see most models for gamma-ray production in MQs invoke processes occurring in 
the relativistic jets. 
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(Mahadevan 1997). 

The soft synchrotron photons inverse Compton scatter off the hot electrons in 
the ADAF and produce harder radiation extending up to about the electron tem- 
perature ~ 100 keV (/i^max ~ kT e ). The relative importance of this process depends 
on the mass accretion rate. At high m, because of the increased optical depth, the 
IC component dominates the spectrum. 

As m decreases, Comptonization becomes less efficient and the IC component 
of the spectrum becomes softer and less important. At low m, the X-ray spectrum 
is dominated by Bremsstrahlung emission, which again cuts off at the electron tem- 
perature (his ^ ax »i kT c ). 

Gamma-ray emission from an ADAF is via the decay of neutral pions produced 
in proton-proton collisions. The results depend sensitively on the energy spectrum 
of the protons. If the protons have a thermal distribution, the gamma-ray spectrum 
is sharply peaked at ~ 70 MeV, and the luminosity is not very high. If the protons 
have a power-law distribution, the gamma-ray spectrum is a power-law extending 
to very high energies, and the luminosity is much higher. The photon index of the 
spectrum is equal to the power-law index of the proton distribution function. 

Another outstanding feature is that ADAFs are much less luminous than thin 
disks at low values of m. This is because most of the energy in an ADAF is advected, 
rather than radiated, leading to a low radiative efficiency. In fact, the luminosity of 
an ADAF scales roughly as ~ Th 2 . In a thin disk, on the other hand, the luminosity 
scales as m (Narayan, Mahadevan & Quataert 1998). 

In determining the radiation processes in ADAFs, the electrons are assumed to 
be thermal, while the protons could be thermal or non-thermal. It is clear that 
the spectrum will depend significantly on the energy distribution of the particles. 
Mahadevan & Quataert (1997) considered two possible thermalization processes in 
ADAFs: (1) Coulomb collisions and (2) synchrotron self-absorption. In the case 
of the protons they found that, for all accretion rates of interest, neither Coulomb 
collisions nor synchrotron self-absorption lead to any significant thermalization. 
The proton distribution function is therefore determined principally by the char- 
acteristics of the viscous heating mechanism, and could therefore be thermal or 
non-thermal. 

ADAF models have been applied to a number of accreting black hole systems. 
They give a satisfying description of the spectral characteristics of several quiescent 
black hole binaries (Narayan et al. 1996; Hameury et al. 1997) and low luminosity 
galactic nuclei (Narayan et al. 1998; Manmoto et al. 1997) which are known to ex- 
perience low efficiency accretion. ADAF models have also been applied successfully 
to more luminous systems which have higher radiative efficiencies (Esin et al. 1997, 
1998). 
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As we already mentioned, the Disk - Corona model is not just another solution 
of the set of equations that the standard disk or the ADAF fulfill; an extra compo- 
nent is added to the already defined thin disk, "the corona". Likely a low density 
corona is heated by reconnecting magnetic loops emerging from the disk (Galeev 
et al. 1979). This implies that the corona is coupled to the disk by the magnetic 
field (for alternative models, where the corona accretes from above the disk, see, 
e.g., Esin et al. 1998). According to the model for the corona formation by Galeev 
et al. (1979), a seed magnetic field is exponentially amplified in the disk due to 
a combination of the differential Keplerian rotation and the turbulent convective 
motions. The amplification time-scale at a radius r is given by £q ~ r/3v c where 
v c is a convective velocity. They showed that inside luminous disks the field is not 
able to dissipate at the rate of amplification. Then bu oyant magnetic loops ele- 
vate to the corona where the Alfvenic velocity is high and the magnetic field may 
dissipate quickly. The magnetic reconnection-heated corona process has therefore 
two phases (Liu et al. 2002a): a mass evaporation happens at the bottom of the 
magnetic flux tube and consequently builds the corona up to certain density; this is 
followed by Compton scattering cooling that radiates away the previous magnetic 
heating. Liu et al. (2002a) showed the dependency of this forming and heat- 
ing mechanism of the corona with the accretion rate. They concluded that for low 
luminosities [L < 0.1L Edd ) most of the accretion energy is transferred to the corona. 

An alternative mechanism is proposed by Liu et al. (2002b) as a frictionally 
heated corona. The existence of the corona is assumed in their model. Heat re- 
leased by friction in the corona flows down into lower and cooler parts. If there the 
density is low, cool matter is heated up and joins the coronal gas. Liu et al. (2002b) 
found that a larger value of the viscosity implies an increase of the evaporation ef- 
ficiency. The evaporated gas loses angular momentum because of the friction and 
drifts towards the central object. This is compensated by a steady mass evaporation 
flow from the underlying disk. Energy from the gravitational potential is released 
by friction in the form of heat in the corona. 

Before going deeper into corona models it is convenient to mention some rele- 
vant spectral features that were the primary motivation for their formulation. We 
will also introduce some basic but useful notions of plasma physics in this context. 

Broad-band X/7-ray spectra of galactic black holes can be explained in terms of 
successive Compton scatterings of soft photons - Comptonization- in a hot electron 
cloud. Besides, both galactic black holes and Seyfert galaxies (a specific case of 
AGN) show a hardening of the spectra at ~ 10 keV, which is attributed to Comp- 
ton reflection (combined effect of photo-electric absorption and Compton down- 
scattering) of hard radiation from a cold material (White, Lightman & Zdziarski 
1988; George & Fabian 1991). This implies the presence of that close cold mate- 
rial. Hard radiation, reprocessed in the cold matter, can form a significant fraction 
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of the soft seed photons for Comptonization. The energy balance of the cold and 
hot phases determines their temperatures and the shape of the emerging spectrum 
(Poutanen & Svensson 1996). 

The situation becomes more complex when a notable fraction of the total lumi- 
nosity escapes at energies above ~ 500 keV. Then hard photons can produce e + e~- 
pairs which will be added to the background plasma. Electrons (and positrons) 
Comptonize soft photons up to gamma-rays and produce even more pairs. Thus, 
the radiation field, in this case, has an influence on the optical depth of the plasma, 
which in turn produces this radiation. This makes the problem non-linear. 

Another complication appears when the energy distribution of particles starts to 
deviate from a Maxwellian. In the so-called non-thermal plasma models, relativistic 
electrons are injected to the soft radiation field. 

Let us first assume a thermal pair plasma. We consider properties of an electron- 
positron plasma cloud in energy and pair equilibria. There are four parameters that 
describe the properties of hot thermal plasmas: 

1. /h = LhO"r/(m e c 3 r c ), the hard compactness which is the dimensionless cloud 
heating rate. L^, is the heating rate of the hot cloud, o^, the Thomson 
scattering cross-section and r c , the cloud size. 

2. The soft photon compactness, / s , which represents the cold disk luminosity 
that enters the hot cloud (corona). L s is the luminosity of the seed soft 
photons that cool the plasma. 

3. t p , the proton (Thomson) optical depth of the cloud (i.e., the optical depth 
due to the background electrons). 

4. The characteristic temperature of the soft photons, T bb . 

Pietrini & Krolik (1995) proposed a very simple analytical formula that relates 
the observed X-ray spectral energy index to the amplification factor, lh/l s ' 



When l h /l s increases, the source becomes more "photon starved" and the observed 
spectrum becomes harder. 

Let us now consider various geometrical arrangements of the hot plasma cloud 
and the source of soft photons. If soft seed photons for Comptonization are pro- 
duced by reprocessing hard X/7-ray radiation, then the geometry will define the 
amplitude of feedback effect and the spectral slope (Liang 1979). What would be 
the most probable geometry taking into account all those concepts and observation 




1/4 
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features?. The main possibilities are: 

" Sandwich" : The simplest solution is to assume that a hot corona covers most 
of the cold disk (a sandwich, or a slab-corona model). The radiative transfer in such 
a geometry was considered by Haardt & Maraschi (1993) who showed that in the 
extreme case, when all the energy is dissipated in the corona, the emitted spectra 
resemble those observed in Seyfert galaxies. Dissipation of energy in the cold disk 
(with subsequent additional production of soft photons) would produce too steep 
spectra in disagreement with observations. Even harder spectra observed in Galac- 
tic BHs cannot be reconciled with the slab-corona model (Gierlihski et al. 1997), 
and alternative models with more photon starved conditions and smaller feedback 
of soft photons are sought. 

"Magnetic flares": A patchy corona (Galeev et al. 1979; Haardt et al. 1994), 
where the cold disk is not covered completely by hot material, has certainly a 
smaller feedback, and the resulting spectra are harder. A patchy corona can be 
described by a number of active regions above the cold accretion disk. Both patchy 
and slab-corona models predict an anisotropy break (i.e. a break in the power-law 
spectrum due to the anisotropy of the seed photons) that should appear at the 
energy corresponding to the second scattering order. 

" Cloudlets" : Another possible solution of the photon starvation problem is to 
assume that the cold disk within the hot corona is disrupted into cold dense opti- 
cally thick clouds (Lightman 1974; Kuncic et al. 1997) that are able to reprocess 
hard X/7-ray radiation and produce soft seed photons for Comptonization. If the 
height-to-radius ratio of the hot cloud is small, we can approximate this geometry 
by a plane-parallel slab. We assume further that the cold material is concentrated 
in the central plane of the hot slab and has a covering factor f c . Compton reflection 
comes from these cold clouds (cloudlets) as well as from the outer cooler disk. The 
seed soft photon radiation is much more isotropic and the emerging high energy 
spectrum does not have an anisotropy break. The covering factor defines the am- 
plitude of the feedback effect. The total soft seed luminosity (with corresponding 
compactness, / s ) is the sum of the reprocessed luminosity and t he luminosity in- 
trinsically dissipated in the cold disk (with corresponding compactness, /g ntr )- 

"Sombrero" : In this model, the cold disk penetrates only a short way into the 
central coronal region (see, e.g., Bisnovatyi-Kogan & Blinnikov 1977, and Poutanen 
et al. 1997). We can assume that the X/7-ray source can be approximated by a 
homogeneous spherical cloud of radius r c situated around a black hole (probably, 
a torus geometry for a hot cloud would be more physically realistic, but then it 
would be more difficult to compute the radiative transfer). The inner radius, r- m , 
of the cold geometrically thin, infinite disk is within the corona (r in /r c < 1). This 
geometry is also similar to the geometry of the ADAF model already discussed. 
Spectra from the sombrero models are almost identical to the spectra expected 
from the cloudlets model, with the only difference that the amount of Compton 
reflection would be a bit larger for the same configuration of the outer cooler disk. 
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From the observational point of view, these models are almost indistinguishable. 



Let us add a new ingredient, by assuming that the plasma can be a "hybrid 
thermal/non-thermal pair plasma" . There are reasons to believe that in a physically 
realistic situation, the electron distribution can notably deviate from a Maxwellian. 
A significant fraction of the total energy input can be injected to the system in form 
of relativistic electrons (pairs). In this hybrid model, the injection of relativistic 
electrons is allowed in addition to the direct heating of thermal electrons and then 
the particle energy distribution is approximately a Maxwellian plus a power law tail. 

The most important input parameters of this version of the model are: 

1. The thermal compactness, Z t h, which characterizes the heating rate of electrons 
(pairs). 

2. The analogous non-thermal compactness, Z nt h, which characterizes the rate of 
injection of relativistic electrons. 

3. The soft photon compactness, Z s . 

4. Tinj, the power-law index of the non-thermal electron injection spectrum. 

5. t p , the proton (Thomson) optical depth. 

6. T bb 

By Zh = Zth + Znth, we now denote the total hard compactness. 

Following this model, the spectrum of escaping radiation then consists of the 
incident blackbody, the soft excess due to Comptonization by a thermal population 
of electrons and a power-law like tail due to Comptonization by a non-thermal elec- 
tron (pair) population (see Fig. 12. 7J) . 

For large Zh/Z s , most of the spectrum is produced by Comptonization of a ther- 
mal population of electrons (pairs), while the tail at energies above m e c 2 is produced 
by non-thermal electrons. For low Zh/Z s , the resulting spectrum is produced by a 
single Compton scattering off non-thermal electrons. 



2.2.4 The Jet 

The main observational aspect that determines the identification of jets as such 
is their radio emission (see Fig. 11.21 and Fig. I1.3J1 : high brightness temperatures, 
"non-thermal" spectra and high degree of linear polarization measurements indicate 
an origin as synchrotron emission from relativistic electrons. However, CHANDRA 
images of moving X-ray jets from the black hole transient XTE J1550-564 (Corbel 
et al. 2002 - see Fig l2.8|) made clear that the non-thermal electromagnetic radiation 
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Figure 2.7: Spectra from the hybrid pair plasmas. Solid curves show dependence on 
hard compactness lh- Other parameters: l\Jl s = 10, Z n thAh = 0.1, t v = 1, T^j = 2.5, 
7bb = 0.2 keV. The resulting electron temperature and optical depth are (kT e ,TT )=(126 
keV, 1.0002), (123 keV, 1.02), and (82 keV, 147) for l h = 1, 10, 100, respectively (l h 
increases from the bottom to the top of the figure). For a higher compactness, the spectrum 
has a sharper cut-off at energies above 1 MeV due to larger optical depth for photon-photon 
pair production. These spectra are similar to the spectra of Galactic BHs in their hard 
state. Dashed curves show dependence on lh/l s - Here we fixed 1^ = 10, l n th/hi = 0.5. The 
resulting electron temperature and optical depth (kT e ,TT) are (104 keV, 1.07), (34 keV, 
1.02), and (5 keV, 1.01) for lh/l s = 10,1,0.1, respectively. Increase in 1^ results in a 
more pronounced blackbody part of the emerging spectrum. The blackbody is modified by 
Comptonization on thermal electrons (from Poutanen 1998). 



from X-ray binary jets may extend to at least the X-ray band 12 . The spectrum of 
these moving jets ranges from radio to X-ray wavelengths and can be fitted by a 
single power law with spectral index a = —0.660 ± 0.005, consistent with optically 
thin synchrotron e mission. Corbel et al. (2002) deduced that under minimum 
energy conditions, the magnetic field in the jets is of 0.3 mG, which implies that 
the leptons emitting in soft X-rays have been accelerated to TeV energies through 
the interchange of bulk kinetic energy to particles via shocks. 

Detailed investigations of the jets from X-ray binaries, both in the radio band 
and at shorter wavelengths, have revealed certain features that seem to indicate a 



12 This observation corresponds to a transient relativistic jet. No direct evidence of X-ray non- 
thermal electromagnetic radiation has yet been found in the case of steady jets. 



2 MICROQUASARS: MAIN FEATURES 



strong coupling between the accretion and the outflows occurring close to relativis- 
tic objects. Let us review some of the main characteristics; interpretations will be 
developed in the next section. 

Observed bright events reveal an optically thin spectrum above some frequency, 
from which the underlying distribution of the electron synchrotron emitting popu- 
lation can be derived. If the electron distribution is a power law of the form 

N(E)dE oc E~ p dE, (2.29) 

then determinations of the spectral index in the optically thin part of the syn- 
chrotron spectrum, a, 

S u oc z/\ (2.30) 

can be used through the relation p — 1 — 2a to find the index of the particles. 
Observed optically thin spectral indices —0.4 > a > —0.8, indicate 1.8 < p < 2.6. 
This is the same range derived for the majority of AGN jets and also for synchrotron 
emission observed in other astrophysical scenarios (e.g. supernova remnants) and is 
consistent with an origin for the electron distribution in diffusive shock acceleration 
(e.g. Longair 1997; Gallant 2002). 

Observations of apparent superluminal motion in galactic jets is a characteristic 
that we already mentioned and developed in Section 11.11 (and with more details 
in Appendix . Apparent superluminal motions, /?t > 1, require intrinsic velo- 
cities (3 > 0.7, indicating the presence of relativistic bulk motions. The associated 
relativistic Doppler shift is given by 

D = r _1 (l =f j3 cos0) _1 , (2.31) 

where is the angle to the line of sight and V is the bulk Lorentz factor of the flow; 
the =F signs refer to approaching and receding components respectively. 

In addition to the proper motions and Doppler-shifting of frequencies, there is 
a boosting effect (or de-boosting depending if we refer to the approaching or re- 
ceding component) due to the relativistic Doppler factor expressed in (|2.31|) . An 
object moving at angle <ft to the line of sight with velocity (3 (and resultant Lorentz 
factor T) will have an observed surface brightness D k ~ a times brighter (or fainter 
in the receding component case) than in the proper system, where 2 < k < 3 and 
k = 2 corresponds to the average of multiple events in e.g. a continuous jet, and 
k = 3 to an emission dominated by a singular evolving event, and a as defined above. 
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Figure 2.8: X-ray emission from large-scale jets produced during an outburst of the 
black hole transient XTE 1550-564- The leptons radiating in soft X-rays should have TeV 
energies (from Corbel et al. 2002). 
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Until now it has been assumed, quite reasonably in the absence of other informa- 
tion, that the jet inclination is perpendicular to the plane of the binary. However, 
at least two jet sources, GRO J1655-40 and V4641 Sgr, appear to show significant 
misalignments (Maccarone 2002). The clearest example of a precessing jet is SS 
433. The ~ 162.5-day precession of these jets has been assumed to reflect the pre- 
cession period of the accretion disk (see e.g. Ogilivie & Dubus 2001). Hjellming 
& Rupen (1995) suggested a precession period for GRO J1655-40 which was very 
close to the subsequently determined orbital period; similarly there seems to be 
marginal evidence for precession in the jets of GRS 1915+105 (Fender et al. 1999), 
LS I +61 303 (Massi et al 2004), and Cyg X-l (Stirling et al. 2001; Pooley et al. 
1999; Romero et al. 2002). 

Because in X-ray binary jets only the synchrotron emission from the leptonic 
(electrons and/or positrons) component has been detected, there is little direct in- 
formation on the baryonic jet content (or lack of it). This is also valid in the case of 
AGN. The exception is SS 433, whose jets are associated with a variety of emission 
lines in optical, infrared and X-rays (e.g. Margon 1984), from where the presence 
of ions can be inferred. 

Finally, a comment about the distinction between jets from neutron stars and 
those from black holes. Fender & Kuulkers (2001) have found that defining a 
quantity "radio loudness" as the peak radio flux of transients (in mJy), divided 
by their peak X-ray flux (in Crab), black hole transients are more radio loud than 
neutron stars. Black holes seem to be about one or two orders of magnitude more 
radio loud than neutron star systems. 



A very active research field related with jets is the investigation of their ori- 
gin, how they may be generated. Blandford & Payne (1982) provided a seminal 
model in which magnetic fields anchored in an accretion flow may produce "ra- 
dio" jets by magneto centrifugal effects. Fig. 12.91 shows four different configurations 
of accretion with magnetic fields which may induce the formation of jets. In all 
these theoretical models a magnetized accretion flow is the basis for jet formation 
by magnetohydrodynamic (MHD) processes. The likely key role of magnetic fields 
in the coupled accretion - outflow system have been reviewed by Meier et al. (2001). 

There are two ways to study these MHD jet formation processes: 

• One is to straightforwardly set up and solve the full set of eight partial differ- 
ential MHD equations: continuity equation (1), Navier-Stokes equation (3), 
induction equation (3) and energy equation (1), making as few assumptions 
as possible. Such demanding simulations run for several hours or days on 
large supercomputers, calculating the evolution of a rotating disk or stellar 
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Figure 2.9: Four ways to make jets with magnetic fields. A: dipole field of a rotating 
neutron star. B: A collapsing object drawing and winding up an initially uniform field. 
C: Poloidal magnetic field from a magnetized accretion disk. D: Frame- dragging near a 
rotating black hole resulting in strong coiling of the magnetic field lines. Types C and D 
(possibly also A) may be relevant for X-ray binaries; type A for isolated pulsars; types C 
and D for AGN (from Meier et al. 2001). 



core of the magnetized plasma and showing the effects of the field on the flow 
and the flow on the field. Time evolution of multidimensional and supersonic 
processes, such as jets and shock waves, can be followed in detail. 

• The other approach is the semi-analytic one. The main idea is to make several 
simplifying assumptions in order to reduce the set of eight partial differential 
equations that can then be solved on a modest computer. The most widely 
used assumptions are: i) Plasma flow and magnetic field structure do not 
change with time - steady state assumption. Such solutions clearly focus only 
on a long-term structure that forms after transient effects due to initial con- 
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ditions have died away, ii) Simplification of the geometry - self- similarity, 
i.e. the structure scales with the distance from the compact object, so that it 
looks similar at every radius. 

Historically speaking, the semi-analytic was the first approach to be considered; 
even now, in the era of modern computers, semi-analytic models still complement 
supercomputer numerical simulations. The simplifications of the semi-analytic ap- 
proach may lead to less realistic solutions. For example, it is more difficult to 
generate highly relativistic narrow jets with such approach (Meier et al. 2001). 
However, because the results are expressed in terms of simple parameters (central 
object mass, magnetic field strength, etc.), they are often applicable to many dif- 
ferent situations. On the other hand, whereas semi-analytic studies often ignore 
the effects of an inner accretion disk boundary or of a central black hole horizon, 
and they usually do not treat the time-dependent nature of the flow. In the case 
of numerical simulations, boundary conditions are taken into account, but care- 
ful attention should be put when choosing them, because they can cause spurious 
effects on the nu merical results. Though having a careful perfect control of the 
initial and/or boundary conditions, numerical simulations lack of a monitoring of 
the physics developed step by step, whose control is the main philosophy of the 
semi-analytic methods. 

Besides the differences between these two approaches, that are definitely to be 
used in a complementary way, they have a common essence with respect to the 
basic steps existing in the MHD jet generation approach, 

1. Generation from a seed magnetic field of a specific structure, always containing 
a poloidal component, B p , and a toroidal component, B^, which state the 
existence of rotational symmetries. 

2. Magnetic braking of the existent angular momentum (of the disk, of the com- 
pact object or of the young star). 

3. Transfer of the extracted angular momentum and rotational energy through 
torsional Alfven waves (TAWs). 

4. Transfer of the magnetic energy from the TAWs (Poynting flux) to some of 
the present charged particles, pushing them upwards and outwards. 

5. Collimation of those particles extracted through pinch effect. 

A rather different but still magnetically working model is proposed in Tagger & 
Pellat (1999). Radiative acceleration (O'Dell 1981) as an alternative to magnetic 
acceleration, is unlikely to be able to push jets to the highest observed bulk veloci- 
ties (Phinney 1982) but may still be operating in the case of SS 433 (Shapiro et al. 
1986). 

It seems quite unlikely that the jet formation region in X-ray binaries will ever 
be directly imaged. Junor et al. (1999) report the direct imaging of jet formation 
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around the (low-luminosity) AGN M87 at a distance of ~ 100 Schwarzschild radii 
from the black hole. Comparing M87 to X-ray binaries in our own galaxy, the ratio 
of distances is so much smaller than the ratio of black hole masses and therefore 
than the one between Schwarzschild radii, that such imaging will not be possible. 
For example, a structure of size 100 Schwarzschild radii around a 1OM black hole 
at a distance of 5 kpc would have an angular size of ~ 10 -9 arcsec (Fender 2004). 
Therefore the key for studying jet formation in X-ray binaries will remain in careful 
multiwavelength studies at the highest time resolution, such as those performed 
with such a success on GRS 1915+105 (e.g. Mirabel et al. 1998; Klein-Wolt et al. 
2002). 

Another family of mechanism for launching jets invoke one-shot acceleration in 
large differential drops established close to the BH (e.g. Blandford & Znajek 1977). 
However, it has been recently argued that the mechanism is not efficient when the 
MHD equations are solved in a regime of strong gravity (Punsly 2001). There is an 
open discussion about the validity of these models. 



2.3 The Dynamics: Different Spectral States - A 
Cycle? 

XRBs are observed in different spectral states. The two most distinctive states 
(and somehow opposed in their characteristics) are the so-called "Low/Hard" and 
"High/Soft" states. These and the other existing states are sketched in Fig. 12.101 as 
a function of the total mass accretion rate, rh. It has been proposed by different au- 
thors (Narayan 1996b; Esin et al. 1997, 1998; Fender et al. 2004b with new variants) 
that different spectral states observed in XRBs can be understood as a sequence of 
thin disk plus coronal region models with varying rh and r in (inner thin disk radius). 

A brief description of each state follows: 

• Quiescent state: This is the lowest luminosity state and has rh < 10~ 2 . Due 
to the low accretion rate, Comptonization is weak, and the X-ray flux is 
much lower than the optical flux. The radiative efficiency is very low and the 
systems are extremely dim. 

• Low/Hard state: For rh above 10~ 2 and up to ~ 10 _1 the geometry of the 
accretion flow is similar to that of the quiescent state, but the luminosity and 
radiative efficiency are larger (and increase rapidly with rh). Comptonization 
becomes increasingly important, giving rise to a very hard power law spectrum 
which peaks around 100 keV. Even though in this state the coronal emission 
predominates, an excess at energies < 1 keV is interpreted as radiation from 
the retracted thin disk (Balucihska-Church et al. 1995) with a characteristic 
temperature, kT^, of order 0.1 - 0.3 keV (usually quite poorly determined, 
due to strong interstellar absorption in that spectral range). 
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• Intermediate state: At still higher accretion rates, m approaches m cri t ~ 0.1, 
the coronal region progressively shrinks in size, the inner radius decreases, and 
the X-ray spectrum changes continuously from hard to soft. This occurs at 
roughly constant bolometric luminosity. In this state, the thin disk becomes 
radiatively comparable to, or even brighter than, the coronal region emission. 

• High/Soft state: At still higher accretion rates, fn > m crit , the coronal region 
cannot exist as an independent entity at any radius, and the thin disk comes 
all the way down to the last stable orbit; there is, however, a weak coronal 
region (with m < m crit ). A typical high state spectrum resembles a standard 
thin disk spectrum with a characteristic temperature kT\^ rs 0.5 — 1 keV, 
and a power-law tail due to the weak coronal emission. Although signatures 
of Compton reflection are also observed in this state (e.g. Tanaka 1991), its 
amplitude is much more difficult to determine since it depends on the assumed 
distribution of ionization and detailed modeling of the continuum, which is 
rather curved in the spectral region around the iron edge. 

This would then be a sequence through which the coronal region gains impor- 
tance and then looses it in favor of the thin disk and so on. While the corona 
is in its main activity period, the cold disk slowly regenerates and its edge moves 
inwards up to the last stable orbit, gradually becoming the main site of emission. 
The emission is then mainly soft with a little hard contribution. When the matter 
in the innermost part of the disk is engulfed by the BH there is a quick change 
of state usually accompanied by a mayor radio outburst. There is some evidence 
that seems to indicate a disk/jet coupling, i.e a relation between the inflow and the 
outflow of matter (we will come back to this point with more detail further on in 
this section). The disappearance of the inner part of the thin disk is then followed 
by the formation of the corona and the generation of the jet 13 . 

It is important to remark that in addition to the explanation it gives of the dif- 
ferent possible spectral states, the combination of a cold disk with a Comptonizing 
corona provides also a satisfactory explanation to other observed spectral charac- 
teristics, i.e. a fluorescent iron line at ~ 6.4 keV, an iron edge at ~ 7 keV, and 
the reflection "hump" at 10-30 keV. During the hard state, some of the photons 
from the coronal region are intercepted by the disk and Compton reflected by the 
colder matter back to the observer, producing this way the mentioned features. In 
the High/Soft state both, the iron line and the iron edge, appear to be smeared 
due to probably gravitational redshift and Doppler effect, implying that the cold 
disk extends very close to the central black hole (see Reynolds & Nowak 2003 for a 
comprehensive discussion) . 

* 

Two important issues can be mentioned at this point. On one hand, the presence 
of jets or the lack of them, as well as their characteristics, seem to be also related 

13 Similar sequences have been proposed in the context of AGNs, e.g. Donea & Biermann (1996). 
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Figure 2.10: The configuration of the accretion flow in different spectral states shown 
schematically as a function of the total mass accretion rate in. The coronal region is 
indicated by dots and the thin disk by horizontal bars. (Adapted from Esin et al. 1997) 



with the sequence of spectral states we have described above. On the other hand, 
the coupling between the thin accretion disk and the jet generation is apparently 
supported by a strong correlation between radio and X-ray flux. 

Hannikainen et al. (1998) were the first to notice a broad correlation between 
the radio and X-ray fluxes for the black hole binary GX 339-4 on its Low/Hard 
state. A similar correlation between radio and X-ray fluxes was found for Cyg 
X-l (Brocksopp et al. 1999), and Fender (2001) suggested that the magnitude of 
the radio/X-ray flux ratio was similar for all Low/Hard state black holes. In the 
past years the understanding of this coupling between radio and X-ray emission has 
advanced significantly. Corbel et al. (2000, 2003) and Corbel & Fender (2002) in 
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a detailed long-term study of GX 339-4, have found that the radio emission in the 
Low/Hard state scales as 



-^radio OC Lr^, (2.32) 



where b ~ 0.7 for X-rays up to at least 20 keV. This relation holds over more than 
three orders of magnitude in soft X-ray flux. In fact, by compiling data for ten 
Low/Hard state sources, it was found that in the luminosity range 10~ 5 L Edd < 
Lx — 10~ 2 LEdd ; all systems are consistent with the same coupling with a very small 
scatter (less than one order of magnitude in radio flux). Jets in the Low/Hard state 
are steady and self-absorbed (radio spectral index, a > 0). Bulk Lorentz factors of 
these jets were deduced by Gallo et al. (2003) to be T < 2. 

The radio/X-ray coupling that appears in the Low/Hard state under the form 
()2.32|) induced some authors to suggest that both, radio and hard X-rays, have the 
same origin: synchrotron emission from the compact jet. This jet model was pro- 
posed by Markoff et al. (2001) 14 , where it was applied to the MQ XTE Jl 118+480. 
Since then it was improved and applied to GX 339-4 in Markoff et al. (2003) and 
further generalized in Falcke et al. (2004). The radio emission is well known to be 
generated through optically thick synchrotron from the compact jet. It has been 
pointed out by Corbel & Fender (2002; though apparently firstly noted by Motch 
et al. 1985) that the infrared/optical bands show a spectral break from optically 
thick to optically thin synchrotron emission which seems to indicate that the X-ray 
spectrum is an extension of the second. 

At the beginning, this approach to the hard X-ray emission was received as an 
alternative explanation to the by then established idea that this kind of X-rays 
arise via thermal Comptonization of softer photons in a hot corona (such models 
are described in Section 12.2.31 a different interpretation of the here discussed cor- 
relation, but in the context of the coronal models, can be found in Zdziarski et al. 
2003). Nowadays the need of both mechanisms of hard X-ray production starts to 
be suggested (Markoff et al. 2003, Falcke et al. 2004, Bosch-Ramon et al. 2005 and 
Fender et al. 2004b) since they seem to be able to account for different character- 
istics that are present in the spectra of MQs. 

A third alternative to mimic this high-frequency part of the observed MQs spec- 
tra was proposed by Georganopoulos et al. (2002). In that scenario the hard X-rays 
would be due to EC scattering of photons from the companion star or the accre- 
tions disk by the leptons of the jet (more on this proposal will be said along the 
forthcoming chapters). 

The next step in the sequence of spectral states is the transition from the 
Low/Hard to the High/Soft state is the so-called "intermediate" or recently also 



1 The formalism of this proposal is based on previous works developed by Falcke & Biermann 
(1995, 1999) for AGNs. 
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called "very high" state, VHS/IS, which seems to be associated with the formation 
of discrete, powerful jets. These jets have important differences with the middling 
relativistic steady optically thick ones present in the Low/Hard state: they are op- 
tically thin, transient and have T > 2 (Fender et al. 2004b). 

The VHS/IS is itself divided into hard- VHS/IS and soft-VHS/IS. Fender et al. 
(2004b) showed through radio observations how the compact jet persists into the 
hard- VHS/IS before the outcome of the outbursts that occur in the transition to 
the soft-VHS/IS while X-rays are softening. These outbursts or very relativistic 
transient jets are optically thin at radio as well as X-ray wavelengths. After these 
flares happen the fact that the radio emissions disappear is taken as a proof that the 
jet does not exist anymore in the soft-VHS/IS and nor will it exist in the following 
High/Soft state. 

The components associated with radio flares have higher bulk Lorentz factors 
than the previously ejected steady jet. These components might be due to internal 
shocks that could arise when the fastest discrete ejections catch up the slower steady 
jet. Fender et al. (2004b) proposed that the very relativistic jets that are produced 
in the middle of the VHS/IS are likely to be formed by the coronal matter that is 
ejected, which, in this way, will also disappear in the transition to the High/Soft 
state. 

Before 1998, radio monitoring of GX 339-4 in the Low/Hard state had already 
established the existence of a weak, mildly variable radio counterpart (Hannikainen 
et al. 1998). When in 1998 the source spent a year in the High/Soft X-ray state, no 
radio counterpart was detected despite multiple observations (Fender et al. 1999). 
The source subsequently returned to the Low/Hard X-ray state and the weak radio 
counterpart reappeared (Corbel et al. 2000). This was the strongest evidence that 
in "soft" disk-dominated states the radio jet was either non-existent or more than 
an order of magnitude weaker. This conclusion was supported, among other obser- 
vations of different XRBs, by the detailed studies of GRS 1915+105 reported by 
Klein- Wolt et al. (2002), in which steady "soft" X-ray states are never associated 
with bright radio emission. 

This cycle of spectral states and the jet characteristics in each step are repre- 
sented in a schematic way in Fig 12. Ill 

A second issue of the relation between accretion disks and generation of jets 
has been observed many times in the collimated ejections of GRS 1915+105, which 
provides one of the best studied cases supporting the proposed disk/jet connection. 
In Fig. 12.121 simultaneous observations realized on September 8th 1997 (Mirabel et 
al. 1998) are presented at radio, infrared and X-ray wavelengths. The data show 
the development of a radio outburst, with a peak flux density of about 50 mJy, as 
a result of a bipolar ejection of plasma clouds. Previous to the radio outburst there 
was a former outburst in the infrared. The simplest interpretation is that both 
flaring events were due to synchrotron radiation generated by the same relativistic 
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Figure 2.11: The central box panel represents an X-ray hardness-intensity diagram. HS 
= high/soft state, VHS/IS = very high /intermediate state, LS = low/hard state. The 
lower panel indicates the variation of the bulk Lorentz factor of the outflow with hardness. 
The dashed loop and dotted track indicate the paths that GRS 1915+105 and some other 
transients take in repeatedly hardening and then crossing zone iii - the "jet line" - from 
left to right, producing further optically thin radio outbursts. Sketches around the outside 
illustrate the relative contributions of jet (blue), corona (yellow) and accretion disk (red) 
at these different stages of the spectral states cycle (from Fender et al. 2004b). 



electrons of the ejected plasma. The adiabatic expansion of the jets causes losses of 
energy of these electrons and as a result the spectral maximum of their synchrotron 
radiation is progressively shifted from the infrared to the radio wavelengths. At 
the same time of the emergence of plasma clouds that produce the infrared and 
radio flares, a sharp decay and hardening of the X-ray emission occur. This can 
be noticed at 8.08 - 8.23 h UT in the figure. The X-ray dip is interpreted as the 
disappearance of the inner regions of the accretion disk (Belloni et al. 1997). Part 
of the matter content in the disk is then ejected into the jets, while the rest is finally 
captured by the central black hole. The recovery of the X-ray emission level at 8.23 
h UT is interpreted as the progressive refilling of the inner accretion disk with a new 
supply of matter that will gradually reach the last stable orbit around the black 
hole. Similar effect has been observed in time scales of years in the extragalactic 
source 3C 120 (Marscher et al. 2002). 
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Figure 2.12: Multi-wavelength behavior of the microquasar GRS 1915+105 as observed 
in September 8th 1997 (Mirabel et al. 1998). The radio data at 3.6 cm (blue squares) 
were obtained with the VLA interferometer; the infrared observations at 2.2 micron (red 
squares) are from the UKIRT; the continuous line is the X-ray emission as observed by 
RXTE in the 250 keV range. These observations show the disk/jet coupling, i.e a relation 
between the inflow seen in X-rays and the outflow of relativistic plasma clouds observed 
as synchrotron emission at infrared wavelengths first and later at radio wavelengths (from 
-original version- Mirabel et al. 1998, modified in Paredes 2004b). 

* 



At this point the thesis sustained by some authors that most likely the 70% 
(Fender 2004) of the galactic X-ray binary systems are able to generate jets looks 
more than reasonable considering that: 

• All known BH systems which are either in the Low/Hard X-ray state, or in the 
intermediate, are associated with the formation of a jet (though of different 
kind). Then the majority of known binary BHs are, or have been in the past, 
associated with a jet. 

• The six brightest low magnetic field neutron star systems, the "Z sources", 
are all associated with jets. The lower luminosity, low magnetic field systems, 
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which are dubbed as "Atoll" sources, may be associated with radio emission 
(although as in the case of black holes there may be bright soft states without 
jets), implying that the lack of radio detections of the majority is a sensitivity 
issue (for detailed description on these two kind of sources, see Fender 2004). 

• The high magnetic field neutron stars, including all but two of the accreting 
X-ray pulsars, are not associated with radio emission. 



This basically means that all the systems except the high magnetic field X-ray 
pulsars, and a small number of black hole and neutron star systems which are in 
persistent "soft" states, are expected to have or to have had jets. 



Chapter 




High-Mass Microquasars as Gamma-Ray 
Sources 



ENGLISH3.1 UNIDENTIFIED EGRET SOURCES: THE GALACTIC PLANE GROUP g3 



3.1 Unidentified EGRET Sources: the Galactic 
Plane Group 

As it is mentioned in Section IT21 about 170 out of 271 point like gamma-ray sources 
included in the third EGRET catalog have not yet been clearly identified with ob- 
jects detected at lower energies. Many of these sources form a group that is con- 
centrated towards the galactic plane and well correlated with the spiral arms of 
the Galaxy; something that indicates a significant contribution from Population I 
objects. Several types of objects have been proposed as possible counterparts of 
these gamma-ray sources, including early-type stars (both isolated and forming bi- 
nary systems), accreting neutron stars, radio-quiet pulsars, interacting supernova 
remnants, and black hole candidates (see Romero 2001, 2004 and references therein). 

In this chapter we will present a model for galactic variable gamma-ray sources 
based on the idea of precessing high- mass 1 microblazars 2 . In our model, the jet 
precession will be induced by the gravitational torque of the high-mass stellar com- 
panion, which is in a non-coplanar orbit, on the accretion disk of the compact object 
(BH or NS). We propose that gamma-rays can be produced by external Compton 
(EC) scattering of UV stellar photons of the massive companions by relativistic 
leptons (electrons and/or positrons) far from the base of the jet. We shall discuss 
cases where the particles are in the form of a "blob" ejected from the central engine 
and where the particles are continuously injected in the form of a persistent jet. 
In the case of high-mass binaries, this mechanism for generation of gamma-rays is 
more efficient than Compton scattering upon disk photons or self-synchrotron jet 
emission (Georganopoulos et al. 2002) and, in combination with the gravitational 
effects of the stellar companion, it can be used to explain gamma-ray sources vari- 
able on time scales from weeks to years. 



1 Let us recall the reader from Section 12.2.11 that HMXBs are distributed along the galactic 
plane. 

2 Aharonian & Atoyan (1998) and Mirabel & Rodriguez (1999) proposed that microquasars 
with jets forming a small angle with the line of sight -by analogy with the unified model for 
AGNs- could appear as microblazars, sources with highly variable and enhanced nonthermal flux 
due to Doppler boosting. 
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3.2 Inverse Compton: the Jet - Photon Field In- 
teraction 

Let us consider a binary system where accretion onto the compact object results 
in the production of twin, relativistic e + e~-pair jets propagating in opposite direc- 
tions. A sudden variation in the injection rate can result in the formation of a blob 
(for which we adopt spherical geometry) which propagates down the jet with a bulk 
Lorentz factor T. We shall follow the general treatment derived by Georganopoulos 
et al. (2001) for EC scattering in extragalactic blazars, adapting it to microblazars. 

In the blob frame the relativistic leptons are considered to have an isotropic 
power-law density distribution, 

n'(7') = ^7' ^(71,72,7'), (3.1) 



where 7 ' is the Lorentz factor of the leptons, k is a constant, and 

^(71,72,7 = 



1 for 71 < 7 ' < 72 
otherwise. 



Using the Lorentz invariant ra/7 2 and the relation 7 = D^' between the Lorentz 
factor of the electrons/positrons in the lab frame and the blob frame respectively, 
we get that the electron density in the lab frame is 



n(j) = ^-D 2+ ^~ p P(liD, 72A 7), (3.2) 



where D = [T (1 q= /3cos0)] 1 is the usual Doppler factor: j3 and are the bulk 
velocity in units of c and the viewing angle, respectively. 

The energy distribution of the effective number of electrons is 

N cS { 1 )=n{ 1 )V eS) (3.3) 
where V e g is the effective observed volume of the blob defined by, 



hoo 

Kff = / d 3 x I dt <J(t-x-n/c). (3.4) 

'blob 



To relate V e g with V, the volume of the blob in its frame, we have to rewrite ()3.4j) 
in terms of the spatial coordinates in the blob frame, keeping the time t in the lab 
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frame. Knowing that dx = dx'/T, dy' = dy, dz' = dz and x' = T(x — c/3t), we 
perform the integral over t and we obtain 3 : 

r(i-/?cos0) y blob 

Using at this point expressions (|3.2J) and (|3.5J) into N e s, (I3.3|l . we get, 

kV 

iV cfr (7) = —D^rvPilxD^D^). (3.6) 

47T 

The blob is injected into an isotropic monoenergetic 4 photon field and will then 
produce the inverse Compton upscattering of some of these photons. The lab frame 
rate of IC interactions per final photon energy is, 



dN p 3cttc 
dtde 4e 7 2 



f(x), (3.7) 



where <jt is the Thomson cross section, e and eo are the energies of the scattered 
photons and seed photons in units of m e c 2 . For the case of Thomson scattering 
and assuming isotropic scattering in the electron frame (Blumenthal & Gould 1970), 

/(x) = ^(l-x)P(l/4 7 2 ,l,x), x=-4-. (3.8) 
3 47^0 

The specific luminosity is then obtained by integrating the scattering rate ()3.7j) 
over the particle energy distribution ()3.6|) . and multiplying by the observed photon 
energy em e c 2 and the number density n p = U / eom e c 2 (where U is the energy density 
of the photon field). The limits of integration are 7 min = (e/4e ) 1//2 and 7 max = 72-D. 
The final expression is: 

dL D3+p kVa T cU e ^ 



dedVt 87reo eo 

£ \ "(l+P)/2 
+ I T- 



4e ; (l+p)(3 + p) 



3 To solve the time integral we use: 



5(g(x)) = —— — —S(x — Xi), Xi such that g(xi) — and g'fxA^O. 
, \9{Xi)\ 



In the coming sections and chapters we will apply the model to specific photon fields. Some 
of them will have blackbody distributions and we will work with the peak of such distribution, eo, 
adapting the situation in a first approximation to a monoenergetic photon field. 
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When p > — 1 and 7 m i n <C 7max, the specific luminosity reduces to: 



dL kVa T cU2P~ 1 



dedVt 



7re (l +p)(3 + p) \e 



-(p-l)/2 



(3.10) 



The luminosity per sr is given by L = eL t and the scattered photons have: 



e m in,T = 4e 7?£ )2 
e ma x,T = 4e 7l-D 2 



(3.11) 



In the Klein-Nishina regime 5 the expression for f(x) in ()3.7j) becomes 6 



/(*) 



o i , , 1 o 2 . (4eo7^) S 
2a; In x + x + 1 — 2x + 



2(l + 4e 7^)J 



P(l/4 7 2 ,l,x), 



(3.12) 



x 



4e 7 2 (l-e/ 7 ) 



and hence the procedure to calculate L e requires numerical integrations instead 
of working with an analytical expression like in the Thomson case (|3.9J) . In this 
context: 



Cmin,KN 



Cmax,KN 



;i + 4e o7lJ D) 

46q7 2 2 ^ 2 
(1 + 4eo72£) 



(3.13) 



In the case of a single blob, the interaction time with the photon field, and 
hence the duration of the gamma-ray flare, will be strongly limited by the appar- 
ently superluminal speed of the feature. For instance, in the case of an 07 star 
and a blob with (3 = 0.98 (r = 5), the flux will decrease to 1/e of its original value 

E E 

5 When — — — 2lL. > i ; which is equivalent to 760 > 1. 

{m e c 2 y 

6 Obtained by Jones (1968). 
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in ~ 1 s. Gamma-ray production through EC in blobs in microblazars are then 
transient phenomena difficult to detect on Earth. This is not the case, however, of 
microblazars with persistent jets. 

The amplification due to the Doppler factor in (j3.10j) is given by D 3+p . This 
factor changes to D 2+p when we consider a continuous relativistic flow instead of 
a blob ejection. The change comes from the fact that in the jet-like case the ob- 
served flux is obtained by integrating over the area (and thereafter the luminosity as 
well) of the feature instead of integrating over the volume as we did in the blob case 7 . 



3.3 External Photon Field: the High-Mass Com- 
panion Star 

3.3.1 Gamma- Ray Emission 

Let us now consider some concrete examples of binary systems, where the jet is 
injected into the photon field due to the high-mass companion star. 

We have calculated the spectral energy distribution for a specific model with 
p = 2, bulk Lorentz factor T = 5, viewing angle <fi = 10 deg, photon field of an 
07 star (photon energies ~ 10 eV), and a high-energy cut-off of 72 = 10 3 >> 71. 
The results are shown in Fig. 13.11 where we also show the energy distribution of the 
stellar photons. The important point is that luminosities of ~ xo 36-37 erg s" 1 can 
be obtained in the observer frame at EGRET's energy range, i.e. 100 MeV - 20 GeV. 

Another similar example can be seen in Fig. 13.21 where the important difference 
with the case of Fig. 13.11 is that the leptons of the jet are assumed to have a high 
energy cut-off in the TeV range 8 . This situation corresponds to an electron Lorentz 
factor of 7 = 10 7 5 . The stellar companion is also a high-mass type corresponding 
to an 09 I star. 

The spectral energy distribution of the IC emission has a peak at MeV energies, 
due to the effect of losses on the particle spectrum 9 . At energies above 10 10 eV 
the Klein-Nishina effect becomes important and the spectrum softens significantly. 
The Klein-Nishina spectral energy distribution resulting from a power law electron 
energy distribution is not a power law, and it is not possible to assign a unique spec- 
tral index to it. The photon spectral index at the MeV-GeV band is T ~ 2. At TeV 
energies, where the spectral index is T > 3, the luminosity is L(E > 1 TeV) ~ 10 32 

7 See Kembhavi & Narlikar (1999) and references therein for the calculations that justify this 
change. 

8 We saw in Section T2.2.4I that the existence of TeV electrons in microquasar jets (at least at 
large scale) have been recently inferred from the detection of synchrotron X-ray emission in XTE 
J1550-564 (Corbel et al. 2002). 

9 This physical effect will be described in Chapter 21 
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Figure 3.1: Spectral high-energy distribution of the scattered photons for a microquasar 
injecting a power-law spectrum of electrons in the photon field of the high-mass stellar 
companion (an 07 star in this example). A cut-off at Lorentz factors 7 ~ 10 3 has been 
assumed. The spectral energy distribution of the star is also shown (left top corner). 

erg s" 1 . Hence, one of these microquasars located at a few kpc might be a detectable 
TeV gamma-ray source. For instance some of the unidentified MeV sources detected 
in the galactic plane by COMPTEL (Zhang et al. 2002, 2004) might correspond 
to sources of this type, and they are potential targe ts for high-energy Cherenkov 
imaging telescopes like HESS or MAGIC. 



3.3.2 Precession vs Variability 

The companion star in a high-mass microblazar system not only provides a pho- 
ton field for inverse Compton interactions, but also a gravitational field that can 
exert a torque onto the accretion disk around the compact object. The effect of 
this torque, in a non-coplanar system, is to induce a Newtonian precession of the 
disk. If the jets are coupled to the disk, as it is usually thought, then the precession 
will be transmitted to them. This situation, which should not be confused with 
the geodetic precession (a purely General Relativity effect), has been extensively 
studied in the case of SS433 (Katz 1980) and extragalactic sources like 3C 273 and 
A0235+16 (e.g. Romero et al. 2000, 2003b). 

A sketch of the situation is presented in Fig. 13.31 The disk Keplerian angular 
velocity is, 




(3.14) 
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Figure 3.2: Inverse- Compton spectral energy distribution for a microquasar with a mas- 
sive stellar companion. Leptons in the jet are assumed to have a power-law energy distri- 
bution with an index p = 2 and a high-energy cut-off at multi-TeV energies. Notice the 
softening of the spectrum at high-energies due to the Klein-Nishina effect. 



where M is the mass of the compact object and r<± is the radius of the precessing 
part of the disk. The orbital period is T m and the orbital radius is given by Kepler's 
law, with m the mass of the star: 

rl=^(m + M)Tl . (3.15) 
The angular velocity of the tidally induced precession can then be approximated 

by, 

3Gmcos9 

where 9 is the half-opening angle of the precession cone (Katz 1980, Romero et al. 
2000). We can now introduce a time-parametrization of the jet's viewing angle as 
(see Abraham & Romero 1999): 



(j){t) = arcsin [x 2 + y 2 ] 1/2 (3.17) 

x = (cos # sin 0o + sin # cos 0o sin cOpt) cos 770 — 

— sin 6 cos uipt sin rj 
y = (cos 9 sin O + sin 9 cos 0o sin u p t) sin 770 + 

+ sin 9 cos uipt cos i] Q , 
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Figure 3.3: Precessing jet model. 



where all angles are defined in Fig. 13.31 

In Fig. 13.41 we show the time evolution of the boosting amplification factor of 
the gamma-ray emission for the case of a continuous jet for two different sets of 
geometrical parameters (viewing angle of 10 deg and precession half-opening angle 
of 1 and 10 deg). Time units are normalized to the precession period T. 

We see that the flux density can change by a factor of ~ 6 x 10 3 in a single 
period. A very weak, otherwise undetected gamma-ray microblazar, can increase 
its flux due to the precession and then enter within the sensitivity of an instrument 
like EGRET, producing a variable unidentified gamma-ray source. The duty cycle 
(i.e. the fraction of time in which the source is highly variable) in this example is 
~ 0.2T, so the source could appear in several EGRET viewing periods. 

Just to give a feeling of the magnitudes involved, we mention that for a black 
hole of 4 M , an 071a stellar companion, an orbital period of 10 days, a half-opening 
angle of 10° for the precession, and an accretion disk of ~ 5 10 11 cm (« 5 10~ 2 a.u.), 
we get a precession period T ~ 100 days. The source could then be detectable ~ 70 
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Figure 3.4: Variation of the amplification factor for continuous jet emission as a func- 
tion of time in the precessing microblazar model (for two different opening angles). Time 
units are normalized to the precessing period. 

days per year. This is consistent with some EGRET observations of highly variable 
sources (e.g., the case discussed by Punsly et al. 2000). 

3.4 Discussion 

Taking into account that there exist more than ~ 130 high-mass X-ray binaries de- 
tected so far 10 and that this number should be a small fraction of the total number 





See Section IO 
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of these objects in the Galaxy, it is not unreasonable to expect the existence of a 
few tens of microblazars at mid and low galactic latitudes that could be responsible 
for the variable galactic gamma-ray sources detected by EGRET. In particular, we 
already mentioned that Paredes et al. (2000) have proposed that the microquasar 
LS 5039, a massive X-ray binary with persistent non-thermal radio emission, could 
be physically associated with the gamma-ray source 3EG J1824-1514. In fact they 
suggested that the observed gamma-ray flux is the result of EC scattering of UV 
photons from the high-mass stellar companion. Our calculations support that sug- 
gestion (see also the more recent treatment by Bosch- Ramon & Paredes 2004). 

The recently discovered X-ray transient V4641 Sgr, which seems to harbor a 
~ 9 M -black hole and displays extreme superluminal velocities, could also be a 
microblazar (Orosz et al 2001). Since the stellar companion is a late B-type star, 
external Compton gamma-ray production is not expected to be very efficient in 
this case, but its high-energy emission could fall within GLAST sensitivity. An- 
other interesting candidate, from the theoretical point of view, is the high-mass 
X-ray binary LS I +61 303, which presents a one-sided jet at milliarcsecond scales 
and evidence for a precessing accretion disk (Massi et al. 2001, 2004). The intrinsic 
jet velocity, however, seems not to be very high: ~ 0.4c, but the source location 
is consistent with its identification with a highly variable gamma-ray source 3EG 
J0241+6103 (Tavani et al. 1998, Torres et al. 2001a). A recent discussion of this 
object as a potential gamma-ray source can be found in Massi (2004a, b) 

Perhaps the best way to identified precessing gamma-ray microblazars is through 
the detection of the electron-positron annihilation feature in their spectra. This an- 
nihilation signature should appear as a broad, blueshifted (by a factor D) line in the 
spectrum at a few MeV, exactly within the energy range of IBIS imager of the IN- 
TEGRAL satellite. Due to the precession of the jet, the Doppler factor will change 
periodically with time, and hence the position of the annihilation peak should os- 
cillate in energy in the lab frame around a mean value (for a detailed discussion of 
the phenomenon see Abraham et al. 2001). Chandra X-ray observations of non- 
thermal radio sources within the EGRET location error boxes (a complete list of 
these sources is given by Torres et al. 2001a) could help to find candidates to new 
microquasars (through the detection of X-ray disk emission) , and then INTEGRAL 
exposures could be used in an attempt to find the annihilation line. Such a detec- 
tion would be a remarkable discovery, since it would establish, at a same time, the 
matter content of the jets in microquasars, and would help to clarify the nature of 
some variable unidentified high-energy gamma-ray sources. Unfortunately, INTE- 
GRAL sensitivity seems not to fulfill pre-launch expectations and the detection of 
the annihilation line might require prohibitively long integration times in most cases. 

To finish this chapter, a comment that concerns the fact that not all high-mass 
microblazars should be confined to the galactic plane. Recent direct measurements 
of proper velocities in high-mass microquasars (Ribo et al. 2002) show that these 
objects can present very high velocities and then could be ejected from the galactic 
plane. This is also supported by the discovery of V4641 Sgr at b ~ —4.8 deg 
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(estimated distance: 7.4 < d < 12.3 kpc). In the case of microquasars with high- 
mass companions, which are young objects, we could expect to find them up to 
distances ~ 100 pc (Ribo et al. 2002) or even more (e.g. V4641 Sgr) from the 
galactic plane. 
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4.1 External Photon Fields: the Disk and the 
Corona 

We briefly recall the reader from Chapter |21 that the soft X-ray blackbody compo- 
nent in the SED of MQs is usually understood as emission from a cold, optically 
thick accretion disk, whereas the power law component is thought to be originated 
in an optically thin hot corona by thermal Comptonization of photons from the disk 
-disk/corona model- or locally produced by synchrotron and Bremsstrahlung radi- 
ation -ADAF model. The hot corona would fill the inner few tens of gravitational 
radii around the black hole. The accretion disk penetrates only marginally in the 
coronal region. We already saw that in the hard state the thermal X-ray emission 
is dominated by the corona, with typical luminosities of a few times 10 37 erg s _1 . In 
the soft states, the disk approaches to the black hole and then most of the energy 
dissipation occurs through it. 

All photon fields that interact with the jet are shown in Fig. 14.11 The photon 
field of the corona is assumed to be isotropic. External inverse Compton interac- 
tions with these photons can be treated in the Klein-Nishina regime since 760 ^> 1, 
with eo the peak of the spectral energy distribution of the coronal photons in units 
of the electron rest energy (eom e c 2 ~ 100 keV). 

In the case of the photons from the accretion disk and the companion star we 
can make calculations in the Thomson regime. Taking into account that the pho- 
tons from the disk come from behind the jet, an additional factor (1 — cos 4>Y P+1 ^ 2 
that reduces the effects of beaming has to be introduced (Dermer et al. 1992). 
The stellar photons, instead, can be treated as an isotropic field, at least as a first 
approximation as we already saw when we studied the interaction of the jet with 
this photon field in Chapter El 

The emerging spectrum of the specific luminosity can be approximated by a 
power law of index (p — l)/2 in the Thomson regime, whose complete expression 
is given by f)3. lOj) . In the Klein-Nishina regime, where numerical integrations are 
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Figure 4.1: External photon fields to which the jet is exposed. 



necessary, the results significantly depart from a power law, resulting in a softer 
spectrum (see Fig. 



The Compton losses in the different regions will modify the injected electron 
spectrum, introducing a break in the luminosity power law at the energy at which 
the cooling time equals the escape time. This will occur at (e.g. Longair 1997): 



7b 



3m e c 2 



Aa T T 2 Ut c 



(4.1) 



Here t e sc is the average time spent by the particles in the field region (typically 
^esc ~ l/c, with I the linear size), and the rest of the symbols have their usual mean- 
ings. The spectrum will steepen from an index p to p + 1 for energies higher than 
7b- After the interaction of the jet with the photons from the disk and the corona, 
the modified spectrum will be injected in the stellar photon field region, suffering 
further losses and modifications. 



4.2 Application to Cygnus X-l 

4.2.1 Main Properties of the System and Observations 

Cygnus X-l is the most extensively studied black hole candidate in the Galaxy. It 
is a very bright X-ray binary with a compact object of ~ 10.1 M & and a companion 
09.7 lab star of ~ 17.8 M Q (Herrero et al. 1995), at an estimated distance of ~ 2 
kpc (e.g. Gierlihski et al. 1999 and references therein). As in other sources of this 
type, the X-ray emission switches between soft and hard states, being most of the 
time in the latter. The spectrum in both states can be approximately represented 
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Figure 4.2: The radio image of the continuous jet of Cygnus X-l at 8.4 GHz obtained 
by Stirling et al. (2001) when the source was in the Low/Hard state. 



as the sum of a blackbody plus a power law with exponential cut-off (e.g. Poutanen 
et al. 1997). During the soft state the blackbody component is dominant and the 
power law is steep, with a photon spectral index F ~ 2.8 (e.g. Frontera et al. 2001). 
During the hard state more energy is in the power law component, which is then 
harder, with photon index ~ 1.6 (e.g. Gierlihski et al. 1997). 

Cygnus X-l has a persistent, mildly variable, compact continuum counterpart 
of flat spectrum (e.g. Pooley et al. 1999). During many years, evidence for non- 
thermal radio jets in Cygnus X-l was lacking, despite the efforts of the observers 
(e.g. Marti et al. 1996). Finally, the jet was detected by Stirling et al. (2001) 
at milliarcsecond resolution using VLBA observations - see Fig 14.21 The jet-like 
feature extends up to ~ 15 mas with an opening angle of less than 2 degrees. The 
spectrum seems to be flat, and no counterjet is observed. The total radio emission 
at 8.4 GHz is ~ 11 mJy, with variations of ~ 2 mJy over timescales of 2 days (Stir- 
ling et al. 2001). The average angle with the line of sight, if the jet is perpendicular 
to the disk, seems to be ~ 30° (Fender 2001). 

Recently, the Interplanetary Network detected a transient soft-gamma ray event 
from the general direction of Cygnus X-l (Golenetskii et al. 2002). Analysis of pre- 
vious data indicates that at least other two events were observed during 1995. These 
latter events were also detected by BATSE instrument on the Compton Gamma- Ray 
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Observatory, suggesting that they were originated in Cygnus X-l (Schmidt 2002). 
The luminosities above 15 keV of the outbursts were in the range 1 — 2 x 10 38 erg 
s -1 , much higher than the typical thermal luminosity in the hard state. 



4.2.2 Spectral Energy Distributions 

In this section we suggest that these soft-gamma ray flaring events can be inter- 
preted in terms of non-thermal microblazar activity. We study the effects of the 
interaction of the relativistic jet with the ambient photon fields from the accretion 
disk, the corona, and the companion star, and we calculate the expected non- 
thermal contribution to the keV-MeV spectrum. In this context the recurrent char- 
acter of the events can be explained through variable Doppler boosting originated 
in the precession of the jet. 

A sketch of the injection of a relativistic leptonic jet at a few Schwarzschild 
radii from the central black hole and its subsequent propagation through the ambi- 
ent photon fields is shown in Fig. 14.31 The individual electrons have Lorentz factors 
7 in the lab frame and the flow is assumed to have a bulk Lorentz factor T. In 
accordance to the disk/jet symbiosis model (e.g. Falcke & Biermann 1999, Markoff 
et al. 2001) and the estimated accretion rate of Cygnus X-l, M ~ 1O~ 8 M yr _1 
(e.g. Poutanen et al. 1997), we adopt a mean particle density n = 10 14 cm -3 for the 
jet at 10rg c h from the black hole. In our calculations we have adopted two different 
values for the original electron energy index 1 : a hard index p = 1.5 and a steeper 
index p = 2.3. The jet was assumed forming an angle of 30 degrees with the line of 
sight, and a bulk Lorentz factor T = 5 was adopted. The initial part of the jet was 
modeled cylindrical structure. 

The jet will traverse first the coronal region. The parameters that characterize 
this region and its photon field in the hard state change with time. We shall assume 
here the typical values given in Table |4~T1 Following Poutanen et al. (1997), the 
coronal region was considered as a spherical region of ~ 500 km in radius. 

The results obtained for the two assumed electron distributions are shown in 
Fig. 14.41 and Fig. 14.51 There, we show all the components due to external Comp- 
ton scattering of the different photon fields by the injected electrons in the lab 
frame. This non-thermal contribution should be added to the thermal compo- 
nents shown in Fig. 14.11 in order to recover the total emission. We can see that 
the non-thermal emission is dominated by the up-scattering of the stellar photons. 
Moreover, gamma-rays produced within the coronal region will be mostly absorbed 

'"Let us remember that the electron energy distribution is assumed to be a power law (see 
Section EH|) . 



n( 7 ) = A^2 + p 7 - Pp(7ijDj72Dj7) _ 
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Figure 4.3: Sketch of the general situation discussed in this chapter. A relativistic jet is 
injected close to the black hole in Cygnus X-l. This jet must traverse photon fields created 
by the cold accretion disk, the hot corona, and the stellar companion. Inverse Compton 
up-scattering of some of these photons is unavoidable. Here is the half-opening angle of 
the cone defined by a possibly precessing jet. 

in the local field through pair creation. Using the simple formulae by Herterich 
(1974) along with the adopted parameters for the coronal region, we estimate an 
optical depth r ~ 1 for photons of 5 MeV. The probability for a 10 MeV photon to 
escape from the corona is only 0.1. 

For a viewing angle of 30 degrees, we still have the thermal emission dominating 
by a factor of 5. But if we introduce the gravitational effects of the companion star, 
then precession of the disk should occur. In Fig. 14.61 we show the modification of 
the beaming amplification factor of the external inverse Compton emission for a 
precessing angle of 16.5 degrees. The time axis is normalized in units of the period 
T. We see that there is a variation of about 1 order of magnitude in the emission 
measured in the observer's frame because of the precession. This means that when 
the jet is closer to the line of sight, the non-thermal luminosity can reach values of 
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Table 4.1: Observational characteristics of Cyg X-l in its soft and hard states (from 
Poutanen et al. 1997). 

Parameter Hard state Soft state 

L° bs a 1 x 1(F erg/s 4 x 1(F erg/s 

L obs b 4 x 1Q 37 erg/ / g l x 1Q 37 erg / g 

a c 0.6 1.6 

T s d 0.13 keV 0.4 keV 

E c e 150 keV > 200 keV 

C f 0.4 0.55 



a 0bserved soft luminosity (total luminosity below ~ 1 keV in the hard state and 

below ~ 3 keV in the soft state) 

b Observed hard luminosity 

c Energy spectral index 

d Temperature of the soft component 

e Cut-off energy of hard component 

f Covering factor of the cold matter 

~ 10 38 erg s _1 , as observed in the recurrent outbursts detected by the Interplane- 
tary Network. The transit through the peak of flux magnification can be very fast 
in the observer frame, leading to quick and transient states when the total flux is 
dominated by the non-thermal contribution. 

The angular velocity of the tidally induced precession can be approximated by 
(13.16(1 . For a half-opening angle of the precession cone, 9 = 16.5 deg as we have 
assumed, the observed orbital period of 5.6 days, and a precessing period of ~ 140 
days (Brocksopp et al. 1999, see below), we get a disk size of ~ 3.9 x 10 11 cm, quite 
reasonable for a wind-accreting system like Cygnus X-l. 



4.3 Discussion 

The above outlined model incorporates the different known components of Cygnus 
X-l, accretion disk, corona, stellar companion, and relativistic jet, in an integrated 
picture where transient non-thermal outbursts are a natural and expected result. 
The amplitude of these outbursts can be similar to what has been recently observed 
in some intriguing flaring episodes in this source. 

Our model is different from the model proposed by Georganopoulos et al. (2002) 
not only because we incorporate the effects of precession, but also because we do 
not attempt to explain the bulk of X-ray emission as non-thermal all the time. 
This emission is normally dominated by thermal Comptonization in the hot corona 
around the black hole, except during the microblazar phase, and in this case we in- 
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Figure 4.4: Results of the model for an injection electron spectrum with index p = 1.5 in 
a cylindrical jet forming a viewing angle of 30 degrees. The bulk Lorentz factor is V = 5 
and the electron power law extends from 71 = 2 to 72 = 10 3 . Three different components 
are shown, resulting from the up-scattering of photons from the star (top panel), the 
disk and the corona (bottom panel, solid and dashed lines respectively). Notice that the 
contribution from the coronal photons is not a power law because of the Klein-Nishina 
effect. 



corporate the effects of the interaction of the jet with the corona in our calculations. 
We emphasize that, as it is shown in Fig. 14.41 and Fig. I4.5| during the transient mi- 
croblazar phase the X- and soft gamma-ray spectrum will be softer than in the 
normal hard state, when the coronal emission dominates. This is an unavoidable 
consequence of the steepening produced by Compton losses in the injected electron 
spectrum and can be used to test our proposal, not only through new observations 





of Cygnus X-l, but also of other potentially precessing MQs as LS 5039 and LS I 
+61 303 (Paredes et al. 2002a; Massi et al. 2004). 

Brocksopp et al. (1999) have found multiwavelength evidence for the presence 
of a 142.0 ± 7.1 days period in Cygnus X-l. The optical and X-ray period seem to 
originate in the precession of the accretion disk (Brocksopp et al. 1999), whereas 
the modulation at radio wavelengths is probably produced by the associated preces- 
sion of the jet (see Pooley et al. 1999). The morphology of the extended radio jet, 
with a clear bend, is also consistent with a precession of the inner beam (Stirling et 
al. 2001). The periodic signal in the radio lightcurve, however, is not expected to 
be as strong as at high energies since the magnification factor for the synchrotron 
emission goes as D^ 3+p ^ 2 , whereas as we already saw for the external Compton 
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Figure 4.6: Variation of the amplification factor for the external inverse Compton 
emission from the jet (electron power-law indices of p = 1.5 -dots- and p = 2.3 -solid 
line-) as a function of time in the precessing microblazar model for Cygnus X-l (precessing 
opening angle of 16.5 degrees). Time units are normalized to the precessing period. 

component, it goes as D 2+p . 

The time lag between the two high-energy flares observed in 1995 is ~ 75 days, 
about a half of the value reported by Brocksopp et al. (1999), but since Cygnus 
X-l is a wind-accreting system variations in the period along a span of several years 
are possible. Certainly, more observations on longer time spans are necessary to 
constrain the dynamical models. In the model presented here the duty cycle of 
the microblazar phase is rather small, ~ 10%. Future X-ray observations of non- 
thermal flares can be used for a better determination of the geometric parameters. 

We saw that in the case of Cyg X-l an interesting resulting characteristic is 
that most of the gamma-rays produced within the coronal region will be absorbed 
by pair production. Sooner or later these pairs will annihilate producing a broad, 
blueshifted feature in the MeV spectrum that might be detected by the INTEGRAL 
satellite. It is expected that this satellite will probe Cygnus X-l spectrum and its 
temporal evolution at this energy range, helping to test and constrain the model 
here proposed 2 . 



2 It is expected though not certain, taking into account the remarks we made in Section 13.41 
about the after-launched resulting INTEGRAL sensitivity. 
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5.1 Components of the Model 

The model for gamma-ray production in MQs proposed in Chapters H3 and 0] was 
based on the process of IC upscattering of seed photons from the high-mass stel- 
lar companion, the disk and the corona. In that context the resulting gamma-ray 
emission is dominated by the EC of the stellar photon field. Besides, pair creation 
absorption mechanism in the disk and coronal X-ray field might quench the contri- 
butions from these two regions to the total gamma-ray emission. 

In this chapter we present a new mechanism for the generation of high-energy 
gamma-rays in MQs that is based on hadronic interactions occurring outside the 
coronal region. The gamma-ray emission arises from the decay of neutral pions 
created in the inelastic collisions between relativistic protons in the jet and the ions 
of the stellar wind. The requisites for the model are a windy high-mass stellar com- 
panion and the presence of multi-TeV protons in the jet 1 . The presence of hadrons 
in MQ jets like those of SS 433 has been inferred from iron X-ray line observations 
(e.g. Kotani et al. 1994, 1996; Migliari et al. 2002), although direct and clear 
evidence exists only for this source so far. In what follows we describe the model 
and present the results of our calculations. 



5.1.1 Hadronic Jet 

The general situation discussed in this chapter is shown in Fig. 15.11 A binary sys- 
tem is formed by a black hole and a high-mass early-type star. A relativistic e — p 
jet is ejected perpendicularly to the accretion disk plane. For simplicity, we shall 
assume that this is also the orbital plane, but this condition can be relaxed to al- 
low, for instance, a precessional motion or more general situations as discussed by 
Maccarone (2002), Butt et al. (2003), and Romero & Orellana (2004). 



interactions of hadronic beams with moving clouds in the context of accreting pulsars have 
been previously discussed in the literature by Aharonian & Atoyan (1996). For an early discussion 
in a general context see Bednarek et al (1990). 
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Figure 5.1: Sketch of the general situation in the hadronic windy MQ: relativistic e—p jet 
is injected close to the black hole in a MQ with a high-mass stellar companion. The stellar 
wind penetrates on the jet from the sides. The resulting interaction produces gamma-ray 
emission. In the figure perpendicularity is assumed between the jet and the orbital plane, 
but this particular assumption can be relaxed in a more general situation. 
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The jet axis, z, is assumed to be normal to the orbital radius a. We shall allow 
the jet to expand laterally, in such a way that its radius is given by 
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R(z) = iz\ 



(5.1) 



with e < 1 and zq < z < z mgx . For e = 1 we have a conical beam. The jet starts 
to expand at a height zq ~ a few hundred km above the black hole, outside the 
coronal region. The particle spectrum of the relativistic e — p flow is assumed to be 
a power law, 



N e ,p(E eiP ) — K eiP E ep p , 



(5.2) 



valid for E' e ™ n < E' e p < E' e ™ x , in the jet frame (denoted by the prime). The 
corresponding particle flux will be, 



J'e, P {K, P ) = -^K, P (K, P )- 



(5.3) 



Since the jet expands, the proton flux can be written as: 



(5.4) 



where n > (a value n = 2 corresponds to the conservation of the number of 
particles, see Ghisellini et al. 1985). Using relativistic invariants, it can be proven 
that the proton flux, in the observer (or lab) frame, becomes (e.g. Purmohammad 
& Samimi 2001) 



Jp(E p , 9) 



An \ z 



cK { z»yn ^~ ap+1 [E p - (3 J El - mycosO 



sin 2 # + r 2 cosfl 



PEp 



1/2 



(5.5) 



where T is the bulk Lorentz factor, 6 is the angle subtended by the emerging photon 
direction and the jet axis, and (3 is the corresponding velocity in units of c. The ex- 
ponential dependence of the cross section on the transverse momentum (p t ) of the in- 
cident protons beams the gamma-ray emission into an angle < cp t /m p T ~ 0.17/r 
along the proton direction, hence justifying the assumption that both directions are 
similar. Note that only photons emitted with angles similar to that of the inclina- 
tion angle of the jet will reach a distant observer, and thus 9 can be approximated 
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by the jet inclination angle. 

In order to determine the matter content of the jet we will adopt the jet-disk 
coupling hypothesis proposed by Falcke & Biermann (1995) and applied with success 
to AGNs (see also Falcke & Biermann 1996), i.e. the total jet power scales with the 
accreting rate as: 

Qi = gjMdiskC 2 , (5.6) 



with q s = 10 1 — 10 3 . The number density n' Q of particles flowing in the jet at 
Ro = R(zo) is then given by, 

cnR 2 n > = (5.7) 
m p c z 

where m p is the proton rest mass. From here we can obtain n' : 



< = ^T- (5-8) 



Additionally, 



n'n 



/ V N>(E' p , zo) dE' p . (5.9) 



Then, if E' p ma - X » E' p mm , which is always the case, we have, 



K = n' (a p -l)(E' p ^r p -\ (5.10) 



which gives the constant in the power-law spectrum at z . 



5.1.2 The Wind from the High-Mass Companion Star 

Early-type stars, like OB stars, lose a significant fraction of their masses through 
very strong supersonic winds. Typical mass loss rates and terminal wind veloci- 
ties for O stars are of the order of 10~ 5 M yr -1 and 2500 km s _1 , respectively 
(Lamers & Cassinelli 1999). At the base of the wind, the density can easily reach 
10~ 12 gcm~ 3 . Such strong winds provide a field of matter dense enough as to pro- 
duce significant hadronic gamma-rays when they penetrate on the jet from the sides. 
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The structure of the matter field will be determined essentially by the stellar 
mass loss rate and the continuity equation: M* = 4nr 2 p(r)v(r), where p is the 
density of the wind and v is its velocity. Hence, 
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Aixr 2 v(r) 

The radial dependence of the wind velocity is given by (Lamers & Cassinelli 1999): 

v(r) = v 00 {l- r -^ , (5.12) 



where is the terminal wind velocity, r* is the stellar radius, and the parameter 
( ~ 1 for massive stars. Hence, using the fact that r 2 = z 2 + a 2 and assuming a 
gas dominated by protons, we get the particle density of the medium along the jet 
axis: 



n(z) = - 1 - - (5.13) 

KJ 4nm pVoo {z 2 + a 2 ) \ ^z 2 + a 2 ,' 



The stability of a relativistic jet under the effects of an external wind has been 
recently investigated by Hardee & Hughes (2003) through both theoretical analysis 
and numerical simulations. Their results indicate that jets surrounded by outflowing 
winds are in general more dynamically stable than those surrounded by a stationary 
medium. Note that protons pertaining to the wind can diffuse into the jet medium. 
The wind penetration into the jet outflow depends on the parameter 



zu~vR(z)/D, (5.14) 

where v is velocity of wind, R(z) is the radius of the jet at a height z above the 
compact object, and D is the diffusion coefficient, w measures the ratio between 
the convective and the diffusive timescale of the particles. In the Bohm limit, with 
typical magnetic fields B ~ 1 — 10 G, w < 1, and the wind matter penetrates the 
jet by diffusion. 
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5.2 Gamma- Ray Emission in the Hadronic Model 

The p — p interaction results in the production of pions whose decay chain leads to 
gamma-ray and neutrino generation: 



p + p — > p + p + got + b(ir + + tt ), a and b G N 

7T° ► 2 7 

7T+ — > + fi + — > v p + e + + v e + i>n 

tt~ — ► + \T — ► Vft + e" + u e + 

The differential gamma-ray emissivity from 7r°-decays is: 



g 7 (£ 7 ) = ^o pp {E p )—^ J p (£ 7 , % A . (5.15) 



Here, the parameter t^a takes into account the contribution from different nuclei 
in the wind and in the jet (for standard composition of cosmic rays and inter- 
stellar medium tja = 1.4 — 1.5, Dermer 1986b). J P {E 1 ) is the proton flux dis- 
tribution evaluated at E — E T The cross section a pp (E p ) for inelastic p — p 
interactions at energy E p « 1(LE 7 can be represented above E p « 10 GeV by 
<x pp (£ p ) ^ 30x [0.95 + 0.06 log(£ p /GeV)] mb. Finally, Z^ is the so-called 
spectrum-weighted moment of the inclusive cross-section. Its value for different 
spectral indices a p is given, for instance, in Table Al of Drury et al. (1994). Notice 
that g 7 is expressed in ph s -1 erg -1 when we adopt CGS units. 

The spectral gamma-ray intensity (photons per unit of time per unit of energy- 
band) is: 

J 7 (£ 7 ,#) = I n(P)q 7 (r')d 3 P, (5.16) 
Jv 



where V is the interaction volume. 

Since we are interested here in a general model and not in the study of a par- 
ticular source, the spectral energy distribution 



L* (E 71 6) = E%(E 7 , 9) (5.17) 



is a more convenient quantity than the flux. Using eqs. (|5.5p . ()5.13|) . (|5.15|) and 
f|5.16|) . we get: 
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(5.18) 



This expression gives approximately the 7r°-decay gamma-ray luminosity for a windy 
MQ at energies £7 > 1 GeV, in a given direction # with respect to the jet axis. 

Depending on the characteristics of the primary star and the geometry of the 
system, high-energy gamma rays can be absorbed in the anisotropic stellar photon 
field through pair production, and then inverse Compton emission from these pairs 
can initiate an e ± -pair cascade. This effect has been studied in detail by Bednarek 
(1997) and Sierpowska & Bednarek (2005). The main effect of these cascades is 
a degradation of TeV gamma-rays into a form of softer MeV-GeV emission. Very 
close systems (a ~ 10 11 cm), with O stars and perpendicular jets, can be optically 
thick for gamma-rays between ~ 10~ 2 — 10 TeV. If the jet is inclined towards the 
star, the effect can be stronger (Bednarek 1997). For systems with larger separa- 
tions, the opacity rapidly falls below unity. 



5.2.1 Specific Spectral Energy Distributions 

In order to make some numerical estimates, we shall adopt the specific MQ model 
presented in Table 15.11 The values chosen for the different parameters are typical 
for MQs with O stellar companions, like Cygnus X-l. We shall consider a conical 
jet (e = 1) with conservation of the number of protons (n = 2) and a high-energy 
cut-off for the population of relativistic protons of _£^ max = 100 TeV. The minimum 
distance from the jet to the primary star is a = 70 R Q ~ 5 x 10 12 cm. 

In Fig. l5.2l we show the spectral high-energy distribution for models with proton 
index a p = 2.2 and a p = 2.8, and for different values of the jet/disk coupling pa- 
rameter gj, with the results obtained using a numerical integration routine. We have 
added an exponential-like cut-off at ~ 0.1£p max . Since the forward momentum 
of the protons in the jet is so great, the gamma-rays will be highly beamed, within 
an angle G ~ arctg (R/z). Hence, only hadronic microblazars would be detected 
(see, nonetheless, Romero & Orellana 2004). We assumed a jet inclination with 
respect to the line of sight of 10 degrees. 
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Table 5.1: Basic parameters of the model 



Parameter Symbol Value 

Type of jet e 1 

Black hole mass M bh 10 M ( 



Injection point 


Z() 


50 R* 


Initial radius 


R 


5 R g 


Radius of the companion star 


r* 


35 R Q 


Mass loss rate 


M* 


10~ 5 M yr~ 


Terminal wind velocity 




2500 km s- 1 


Black hole accretion rate 


Miisk 


10~ 8 M yr" 


Wind velocity index 


c 


1 


Jet's expansion index 


n 


2 


Jet's Lorentz factor 


r 


5 


Minimum proton energy 


771/min 


10 GeV 


Maximum proton energy 


771/ max 

e p 


100 TeV 


Orbital radius 


a 


2 n 



*R g = GM hh /c 2 . 



5.3 Production of Neutrinos and Secondaries 

Neutrinos are also generated by pion decay chains (see Section f5.2|) . The signal- 
to- noise (S/N) ratio for the detection of such a //-signal can be obtained analyzing 
the event rate of atmospheric //-background and comparing it with the event rate 
from the source (see e.g. Anchordoqui et al. 2003 for details). The S/N ratio in 
a km-scale detector (like IceCube) in the 1-10 TeV band (including the effects of 
neutrino oscillations) is ~ 3 for one year of operation, assuming an inclination angle 
of 30 degrees, a p = 2.2 and that the neutrino spectrum roughly satisfies (e.g. Dar 
& Laor 1997): 

dF dF 

-0.7^. (5.19) 



dE v dE^ 



This S/N is high enough as to justify speculations on the possibility of detecting 
a hadronic microblazar first from its neutrino signal (a serendipitous discovery in 
a detector like IceCube), and only later from its gamma-ray emission (through a 
pointed observation). 

The here presented hadronic model contemplates the production of gamma-rays 
in MQs as the final stage that follows the decay of neutral pions previously gener- 
ated in p — p interactions. However, neutral pions can also be produced by p — 7 
interactions (also called "photo-meson" processes), which as an alternative model 
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Figure 5.2: Spectral high-energy distribution for windy MQs with proton index a p = 2.2 
and a p = 2.8, for different jet/disk coupling constants (qj). The jet inclination with 
respect to the line of sight is assumed to be 10 degrees. An angle of 30 degrees reduces the 
luminosity for about two orders of magnitude. 



also gives rise to the generation of TeV photons and neutrinos 2 (Levinson and Wax- 
man 2001, Distefano et al. 2002). The site where the emission is generated in this 
case is within the coronal region in the innermost part of MQs' jets. This implies a 
fundamental difference of this model with the one used in this chapter. Pions are 
then produced in collisions of hadrons from the jet with the strong external X-ray 
photons (e.g. from the hot corona) and with the synchrotron photons produced by 
leptons inside the jet. In this process protons are required to be accelerated up to 
very high energies (~ 10 16 eV) in the inner jet, in order to produce a significant 
flux of multi-TeV neutrinos. This requirement might represent an important dif- 
ficulty for the photo-meson model considering that the possibility for particles as 
protons to be accelerated up to such energies in MQs is still under strong discussion. 
Besides, neutrinos are supposed to extract more energy from the original protons 



2 This kind of model is similar to the AGN proton initiated cascade model of Mannheim & Bier- 
mann (1992), Mannheim et al. (1992) and Mannheim (1993). See also Biermann & Strittmatter 
(1987) for the presentation of pioneering ideas in this sense. 
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in the p — p process than in the p — 7 one: £™ ax = j^E^ and E™ x = ^j-E™ ax 
respectively ( Alvarez- Muniz & Halzen 2002). 

It is not possible to separate in the neutrino signal the contributions from the 
photo-meson and p — p channels, but simultaneous X-ray observations can help to 
determine the characteristics of the relevant photon fields to which the inner jet is 
exposed, making then possible estimates of each contribution in particular cases. 
That one channel might dominate over the other will depend on the relative tar- 
get density of the photons and protons in the source region where the protons are 
accelerated (Halzen & Hooper 2002). In any case the detection of TeV neutrinos 
from MQs could be used as an important diagnostic of jets' content. This kind of 
detection would immediately imply the presence of hadrons in the jets. Therefore 
it would help to discriminate e — p jets from pair dominated ones (see Bednarek et 
al. 2004 for additional discussion). 

Another interesting aspect of the hadronic microblazar described in this chapter 
is that e ± are injected outside the coronal region through ^ decays (see Sectioning- 
These leptons do not experience the severe IC losses that affect to primary elec- 
trons and pairs when crossing the disk and coronal photon fields as we mentioned in 
Chapter These secondaries will mainly cool through synchrotron radiation (at X- 
rays, in the case of TeV particles) and IC interactions with the stellar seed photons 
(that would result into an additional source of MeV-GeV gamma-rays). The spec- 
trum of secondary pairs roughly mimics the shape of the proton spectrum. Hence, 
synchrotron emission from these particles will present indices a syn ~ (a p — l)/2, 
which for values of a p ~ 2 are similar to what is observed in MQs' jets at radio 
wavelengths. The losses of the primaries in the inner source lead to a soft par- 
ticle spectrum that is then injected in the region where the particles produce IC 
gamma-rays through interactions with stellar UV photons. Pure leptonic models 
for the gamma-ray flux (e.g. Kaufman Bernado et al. 2002, Romero et al. 2002) 
then require particle re-acceleration in order to explain the flat radio-spectrum of 
sources like LS 5039 far from the core. In the hadronic model, the leptons might 
be injected in situ with the right spectrum through hadronic decays, avoiding the 
problem. 

At TeV gamma-ray energies hadronic microblazars which are optically thin to 
pair production can be detected as unidentified, point-like sources with relatively 
hard spectra. This kind of sources could display variability. In the near future, 
new ground-based Cherenkov telescopes like HESS and MAGIC 3 might detect the 
signatures of such sources on the galactic plane. Hadronic microblazars might 
be part of this population, as well as of the parent population of low latitude 
unidentified EGRET sources. 



3 The VERITAS project is also expected to contribute to the detection of this kind of sources. 



Low-Mass Microquasars and Halo Sources 
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6.1 Variable Gamma- Ray Sources off the Galac- 
tic Plane 

A brief description of the mid-high latitude EGRET sources was included in Chap- 
ter H In what follows, we will describe in some more detail the properties that 
induce to the classification of these sources in three subgroups. The 93 unidenti- 
fied sources detected away from the Galactic plane are displayed in Fig 16. II Their 
concentration at 3° < |fe| < 30° and in the inner half steradian indicates (at a 7a 
confidence level) that 70 to 100 % of them have a Galactic origin, depending on the 
choice of Galactic scale height. Their temporal and spatial characteristics reveal 
an heterogeneous sample. Fig. 16.21 shows their variability index distribution, which 
closely follows that of the variable AGN sources (Nolan et al. 2003). The average 
5 index is 0.79 ± 0.08, 0.66 ± 0.06, and 0.42 ± 0.06 for the halo, AGN, and Belt 
sources, respectively. The halo sources are clearly variable. Taking into account 
the detection biases (Grenier 2000), the spatial distribution of the variable sources 
implies an origin in a thick galactic disk with a scale height of 1.3 ± 0.6 kpc. At 
typical distances of 5 to 10 kpc, the luminosities of the halo sources range from 2 to 
30 x 10 34 erg s _1 sr _1 above 100 MeV. They exhibit large luminosity L 1 /Lx ratios 
of a few hundred. 

The compact objects likely to power sources located high above the plane in- 
clude millisecond pulsars and microquasars with a low-mass star companion, both 
having migrated away from the galactic plane or escaped from globular clusters. In 
the first case we refer to MQs that could have been ejected by kicks imparted in the 
natal supernova explosions to galactic centric orbits (Mirabel & Rodrigues 2003). 
In the second case they might have been born in globular clusters from where they 
might have been injected in large galactic orbits (Mirabel et al. 2001). Until now 
no gamma-ray source is found to be positionally coincident with a globular cluster. 

Pulsed gamma-rays have been detected from the ms pulsar PSR J0218+4232, 
in phase with the radio and X-ray peaks (Kuiper et al. 2002). This object shares 
many traits with the halo sources: a distance of 5.7 kpc and an altitude of 1.6 kpc, 
a luminosity of 1.6 x 10 34 erg s _1 sr _1 and a spectral index of 2.6 above 100 MeV. 
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Figure 6.1: All-sky plot, in galactic coordinates, of the unidentified EGRET sources at 
latitudes \b\ > 3°. The steadier sources associated with the Gould Belt and the variable 
sources are marked as circles and stars, respectively. 

PSR J0218+4232 does not belong to a globular cluster either. Yet, no long-term 
variability is expected from theory (Zhang & Cheng 2003). 

In the next section, we explore whether low-mass microquasars can produce the 
halo sources despite their intrinsic faintness and softness compared with the young 
high-mass systems. The stellar luminosity in LMMQs is reduced by 4 to 6 orders 
of magnitude with respect to the case of a high-mass stellar companion. There is 
also a reduction of 1 order of magnitude in the emission from the accretion disk. 
The thermal emission in both components peaks a decade or two lower in energy. 
Compare Fig. 14.11 with Fig. 16.61 to check these differences luminosities and peak 
energies. 

Only leptonic models can be considered for halo sources since hadronic ones 
require the presence of an early type of companion star in order to ensure the exis- 
tence of strong stellar winds (Romero et al. 2003a). 



6.2 Emission from Low-Mass Microquasars 

XTE J1118+480 is a LMMQ at a high galactic latitude above the plane, b = 62°, 
and a distance of 1.8 ± 0.6 kpc. It lies at 1.6 ± 0.5 kpc above the galactic plane 
(McClintock et al. 2001a). The central object accretes matter through a disk from 
a low-mass star via Roche lobe overflow. McClintock et al. (2001a) and Wagner et 
al. (2001) constrained the companion spectral type to be between K5 V and Ml V. 
A large mass function, f(M) 6 MO, strongly suggests that the compact object is 
a black hole in a fairly compact binary system with a short orbital period of 4.1 hr 
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Figure 6.2: Distributions of the 5 variability indices from Nolan et al. (2003) for the 
firm AGN, Gould Belt, and halo sources at latitudes \b\ > 3°. 



(Wagner et al. 2001). The disk may be precessing under the stellar tidal influence 
(Torres et al. 2002). 

The rapid and correlated UV-optical- X-ray variability in the low-hard state is 
interpreted as a signature of the strong coupling between a hot corona and a jet 
emitting synchrotron radiation up to, at least, the UV band (Hynes et al. 2003; 
Chaty et al. 2003; Malzac et al. 2004). The outflow has remained steady through 
the outburst evolution (Chaty et al. 2003). 

The coronal emission extends up to ~ 150 keV (McClintock et al. 2001b). The 
optical to hard-X-ray data can be modeled by the Comptonization in a hot corona 
or in the inner accretion flow of the soft photons emitted by the outer cold disk 
with an inner radius of ~ 55i?sch and a temperature of ~ 24 eV (Esin et al. 2001; 
McClintock et al. 2001b; Malzac et al. 2004). 



XTE Jl 118+480 therefore serves as a good example for a low- mass microquasar 
and we adopt its characteristics as input to our model (see table 16. lj) . SED of 
the thermal stellar, disk and coronal components, taken from the aforementioned 
publications, are displayed in Fig. 16.61 The case of an F star companion is also 
considered in Fig. 16.71 We calculate the EC emission from interactions of the jet 
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Table 6.1: Parameter set corresponding to the MQ XTE J1118+480 and used in the 
model 



Black hole mass 

Mass accretion rate 

K-M star bolometric luminosity 

F star bolometric luminosity 

Star temperature 

Star orbital radius 

Jet /accretion power ratio 

Corona luminosity 

Corona outer radius 

Corona photon index (dNx / dE oc E 

Corona cut-off energy 

Disk luminosity 

Disk temperature 

Initial jet radius 

Jet bulk Lorentz factor 

Jet viewing angle 

Jet electron index (dN e /dE oc E~ p ) 
Maximum electron energy 
Minimum electron energy 



M bh = 6.5 MO 

M = 3 x 1(T 8 MO yr- 1 

L KM = 4 x 10 32 erg/s 

Lp — 1.5 x 10 34 erg/s 

IcTkm = 1 eV and kTp = 1.8 eV 



D, = 


1.7 x 10 11 cm 




9jet = 


Pjet/Mc 2 = 10 


-3 - 


-^cor — 


= 7.8 x 10 34 erg 


•/b 


Rcor ~ 


= 10 8 cm 




^cor 


: 1.8 




E 


= 150 keV 




-^disk = 


= 8.6 x 10 35 erj 


y/s 




= 24 eV 




-Rjet = 


= 1.9 x 10 7 cm 




Fjet = 


3 to 10 





<f) = r to 30° 

p = 2 to 3 

-Eemax = 5 GeV to 5 TeV 
E emin = 1 to 5 MeV 



with the three radiation fields following the model developed in Chapters El and 01 
We therefore assume a population of e + — e~ pairs in a persistent, cylindrical jet, 
with a power-law distribution in number density per unit energy, N(E e ) = kE~ p , 
between £™ in and £'™ ax . The jet is assumed to be parallel to the disk axis, at an 
angle to the line of sight. It moves with a bulk Lorentz factor T and carries a 
total power Pj et = q^ et Mc 2 (adopting the jet-disk coupling proposed by Falcke & 
Biermann 1995). 

The coronal region is assumed to fill a sphere inscribed in the inner disk radius. 
Let us recall that the Compton losses in the different regions can modify the in- 
jected electron spectrum, introducing a break in the power-law from an index p to 
p+1 at the energy at which the cooling time equals the escape time. In the scenario 
studied here, this occurs with the disk radiation field. We have seen from previous 
chapters that another important ingredient can be the absorption from two-photon 
pair creation in the ambient radiation. It turns out to be quite effective near the 
disk (r > 1) whereas the coronal and stellar fields are optically thin (r < 1) to 
the gamma-rays. These conclusions are supported by Fig l6.3l and Fig 16.41 where we 
show the computed opacity, r, of the stellar photon field on one hand and of the 
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disk and corona fields on the other. 

The ambient radiation from the companion star, E+, has an optical depth for a 
photon with energy E 1 given by (Jauch & Rohrlich 1976): 

r(£ 7 ) = (r) J N{E+) a(E+, EjdE*, (6.1) 

where (r) is the mean dimension of the pair production region and <j(E+, E^) is the 
photon-photon pair creation cross section given by: 



7rr: 



a e+e -(E„ E,) = ^(l-e) 



2£(£ 2 -2) + (3-£ 4 )ln 



1+j 
1-^ 



(6.2) 



where r e is the classical radius of the electron and £ 



1 - 



(jn e c 2 ) 
E+E~ 



,2N2-|l/2 



The seed 



photon field from the star is assumed to be a blackbody type and its differential 
number density, dn{E) = N(E)dE, is (Leung et al. 1993): 



N(E+)dE+ 



EldE* 



7i 2 {hc) 3 {e E */ kT - 1)' 



(6.3) 



The coronal region was considered as a spherical region of ~ 55i?s c h in radius. 
Considering that this area is filled by seed photons from the disk and from the hot 
corona, we calculate the total opacity taking into account its dependency on r: 



/»oo /»oo 

r(£7 7 , r)= / [N disk (E x , r') + N COTOIlll (E x , r')] x 

Jr JO 

xa e - e+ (E x ,E^)dE x dr'. (6.4) 



The disk photon distribution is that of a blackbody, 



firB(E x )\ fl|,„ 



^ r >1wjf' (6 - 5) 

with R disk = 55i? S ch, and 



2E- 



B ^ = 7, , 2(M &^ ~ (6-6) 



(he) 2 (e 



The corona photon distribution is a power law such that: 
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Figure 6.3: The optical depth for a photon with energy E~ in the ambient radiation from 
the companion star. The stellar field is clearly optically thin to the gamma-rays since 
t < 1 for the whole range of E~ . 
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(6.7) 



As it has been discussed, in the observer frame, the IC contributions from in- 
teractions with the stellar and coronal photons to the SED are amplified by D 2+p 
whereas the disk one by D 2+p {\ — cos0)^ 1+p ^ 2 . As it can be seen in Fig. 16.51 the 
amplification factor corresponding to the star and corona peaks along the jet axis 
in = 0°, whereas one for the disk peaks in ~ 15° to 20°, for T ~ 3 and p = 2 to 3. 

In Section POl we clarified the differences between the application of our leptonic 
model to the HMMQ Cygnus X-l and the work presented by Georganopoulos et 
al. (2002) on the same object. Now we have to point out that with respect to 
the LMMQ XTE J1118+480, the scenario we are discussing also differs from that 
of Georganopoulos et al. (2002) who have imposed a much lower energy cut-off 
to the electrons so that the hard X-ray emission results from EC interactions with 
the stellar and disk photons rather than from a hot inner accretion flow or from a 
hot coronal plasma energized by magnetic flares above the disk, elements that are 
absent in their model. 





Figure 6.4: Optical depth for gamma-rays in the inner region of a LMMQ. In the upper 
plot both photon densities, from the disk and the corona, are included in the calculation. 
On the bottom, only the coronal photon density was considered showing that the opacity 
is completely due to the seed photons from the disk. The opacity goes down progressively 
and the last layers from r/55i?sch ~ are already optically thin to the gamma-rays of 
any energy. 
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Figure 6.5: Amplification factors on the SEDs: A = D 2+p , in the cases of the coronal 
and stellar contributions to the SED (top); A = D 2+P (l — cos (f>)( p+1 ^ 2 for the disk one 
(bottom), (p: viewing angle. 



Figures Ifi.fil and 16.71 show the SED obtained per steradian in the laboratory 
frame for a LMMQ with the parameters listed in Table 16.11 adding the contribu- 
tions from the three external photon fields. The EC coronal component (in the 
Klein-Nishina regime) is dominant in the COMPTEL 3-30 MeV band whereas the 
EC disk component (in the Thomson regime) takes over above 100 MeV. This is 
why the maximum gamma-ray luminosity is reached for close to 15°, reflecting 
the angle dependence of the disk amplification factor. This emission has a photon 
spectral index of 2.5 between 0.1 and 1 GeV. The IC stellar component is small 
at all viewing angles. These results show that even though the spectral index in 
the EGRET band matches that of the unidentified sources, the maximum predicted 
luminosity, L max ~ 4 x 10 29 (i?/100 MeV) -0 ' 5 erg s" 1 sr" 1 , is 5 orders of magnitude 
too faint to account for the typ ical halo source fluxes at distances of 5 to 10 kpc. 

The relative contributions of the three EC components change in the case of an 
extreme microblazar where the bulk Lorentz factor F reaches 10, electron energies 
extend up to 5 TeV, and the jet axis is close to the line of sight (0 = 1°). The 
generation of such a highly relativistic outflow, with T > 10, has been recently 
observed from Circinus X-l, a neutron star with a stellar- mass companion as it is 
discussed in Section l2.2.11 Particle energies as high as 10 TeV have been inferred for 
the large-scale jets of the low-mass microquasar XTE J1550-564 (see Section 12. 2. 4|) . 
0.1 pc away from the black hole (Corbel et al. 2002). 
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Figure 6.6: Spectral energy distribution of the EC emission from the jet of a microquasar 
with a K-M star companion, seen at angles of 5° (A), 15° (B), and 30° (C) from its axis, 
for T = 3, (/jet = 10~ 2 ; and a jet electron index p=2.3 with 7™ m = 2 and <-y™ ax = 10 4 . 



Fig. 1(181 shows the result of our calculations for a low-mass microblazar. The 
assumed parameters are indicated in the caption. The corona EC emission predom- 
inates in a short range of energies up to several hundred MeV for an F star and to 
several GeV for a K-M star, beyond which the harder, stellar EC component takes 
over. The disk EC emission is negligible because of the lesser amplification at small 
viewing angle. The ~ 35 times larger energy density that the jet encounters around 
an F star compared with the case of a K-M star results in a modest luminosity 
increment by a factor < 5 at the beginning of the EGRET band up to a maximum 
factor ~ 30 at the end of that range. Anyhow, the very low luminosities, as well 
as the generally too soft spectra, with photon indices ranging from 3 to 4, cannot 
account for the typical characteristics of the EGRET sources. Nor would these 
systems be detected above 100 GeV by the new-generation Cherenkov telescopes 
since the predicted luminosity falls orders of magnitude below their sensitivity in 
all cases, even when multi-TeV electrons are present in the jet. 



6.3 Discussion 

LMMQs can remain quiescent for years or decades before brightening by as much 
as a factor 10 7 in X-rays in a week. The accreting keV luminosity of ~ 10 35 erg 
s -i gr -i ac i pted in the previous section corresponds to such a flaring state, i.e. to a 
bright X-ray source with a vF v flux of 2.6 x 10~ 4 (-D/5kpc)~ 2 MeV cm -2 s _1 . The 
GeV emission obtained is mainly due to the EC of the seed photons from the disk 
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Figure 6.7: Spectral energy distribution of the EC emission from the jet of a microquasar 
with an F star companion, seen at angles of 5° (A), 15° (B), and 30° (C) from its axis, 
for T = 3, gjet = 10" 2 , and a jet electron index p=2.3 with 7™ 111 = 2 and 7™ ax = 10 4 . 

or the corona (depending on the conditions in the jet). The predicted luminosity 
should be considered as an upper limit for comparison with the unidentified sources 
and in that sense EC emission in low-mass microquasars largely fails to explain the 
variable EGRET sources at large scale heights. 

In contrast with high-mass systems where the external radiation energy density 
largely surpasses the magnetic one, SSC emission in low-mass systems is likely to 
dominate even for a modest field strength of ~ 10 G in the jet. For this value 
the magnetic energy density compares with that of a K-M star. Applying the 
model developed by Bosch-Ramon et al. (2005) to a cylindrical jet, we can get 
luminosities similar to those expected from the EGRET halo sources (see Fig. 16. 9|) . 
The luminosity per unit of frequency for the synchrotron radiation in the observer 
frame is: 



where R is the radius of the jet, s v is the source function of the synchrotron emission 
from an isotropic particle distribution and r v is the synchrotron optical depth of 
the jet 1 . 

The resulting SSC emission has a luminosity in the jet frame given by: 



The expressions for s v and t v can be found in Pacholczyk (1970). 




(6.8) 
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Figure 6.8: Spectral energy distribution of the EC emission from the jet of an extreme 
microblazar with an F or K-M star companion, seen at 1° from its axis, for V = 10, 
qjet = 1CT 2 ; and a jet electron index p=2.3 with 7™ m = 10 and 7™ ax = 10 7 . 



f2mn /"EOmax Hemax f/jV" g 

eL e = e / / EU B y n (6 )n(j e )-^—dyde dz, (6.9) 

" 2 min " e 0min ^7cmin ^ ^0 

where £ is the surface of a perpendicular jet slice, - p is the lab frame rate of IC 

cttae 

interactions per final photon energy (see 13 .7|) and the synchrotron radiation 
density is given by: 

£W « — ■ (6-io) 



Fig. l6.9l shows the luminosity of the SSC emission as calculated from the observer 
frame. The integration of (jfi.Dj) is first performed in the co-moving frame and then 
transformed to the lab: 

e'L' e , = D 2+p e'L £ ,. (6.11) 

The energy of the scattered photon in the jet's reference frame, e, is boosted to: 

e' — D e. (6.12) 



The synchrotron part of the spectrum is ~ 100 - 1000 times fainter than the 
near-IR to UV synchrotron emission recorded from XTE J1118+480 during the 
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Figure 6.9: Spectral energy distribution of the synchrotron and SSC emission from the 
jet of an extreme microblazar, with a 10 G magnetic field, seen at 1° from its axis, for 
T = 10, g-jet = 10 -2 , and a jet electron index p=2.3 with 7™ in = 10 and 7° iax = 10 5 . 

outburst 2 , the adopted electron index in Fig. 16.91 being softer than that of XTE 
Jl 118+480 to avoid a too hard spectrum at high energies. The latter softens be- 
yond 10 GeV because of the Klein- Nishina effect. One should keep in mind, however, 
that 77 absorption against the disk photons is not included in the calculation al- 
though it could efficiently limit the emerging gamma-ray flux (Bosch-Ramon et al. 
in preparation). Adiabatic losses in an expanding jet would also modify the result 
although not severely. 

SSC emission in a low-mass microquasar appears then an interesting possibility 
to be investigated in greater detail, bearing in mind that modeling the high-energy 
radiation from these fascinating objects is severely limited by the highly uncertain 
choice of the jet bulk motion and magnetic field. These first results of SSC emission 
in LMMQs, make of these objects a promising parent population for the halo group 
of unidentified EGRET sources. 



2 See Markoff et al. (2001) for a synchrotron model of this specific object. 
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The fact of being a scale-down version of quasars makes of microquasars a very 
interesting type of object, rich in an important variety of physical phenomena that 
were always difficult to study in detail in quasars considering their time scales and 
distances. For example, MQs have already provided an important insight into the 
connection between accretion disk instabilities and the formation of jets, since the 
characteristic times in the flow of matter onto a black hole are proportional to 
its mass, and therefore the accretion-ejection phenomena in quasars should last 
10 5 — 10 7 times longer than the analogous phenomena in MQs. In this way, varia- 
tions on scales of tens of minutes in MQs have revealed events that had been difficult 
to observe in quasars. 

Before the launch of the CGRO, only one AGN had been detected at high-energy 
gamma-rays. We are referring to the quasar 3C 273 that was observed in this wave- 
length by COS B (Swanenburg et al. 1978). After phases I (91/05 - 92/11) and II 
(92/11 - 93/09) of the CGRO, EGRET detected with a high degree of confidence, 
33 AGNs in high-energy gamma rays (E > 100 MeV) (von Montigny et al. 1995, 
and references therein). The strong connection between this high-energy emission 
and the presence of relativistic jets has been widely discussed, for example in the 
mentioned paper from von Montigny at al. (1995) and many other similar refer- 
ences. Following the line of the quasar - microquasar analogy, it looks reasonable 
at least to suspect the possibility that high gamma-ray emission could arise from 
MQs' jets as well. 

In fact, the association proposed by Paredes et al. (2000) of the microquasar LS 
5039 with the gamma-ray unidentified source 3EG J1824-1514, was a first observa- 
tional step towards the identification of the gamma-ray emitting galactic jet sources. 

This was the starting point for the research described in this thesis, also moti- 
vated by the significant number of unidentified gamma ray sources contained in the 
3 rd EGRET catalog. Many of these sources, as we saw in this text , are suspected to 
be of galactic origin because of their spatial distribution. In order to establish mi- 
croquasars as a parent population of some subset of unidentified gamma-ray sources 
we had first to look for a basic theoretical picture that could justify the gamma-ray 
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production in different scenarios involving diverse conditions given in micro quasars. 



We started with a very simple leptonic toy model, in the context of high-mass 
microquasars, where the emission is due to inverse Compton process between rela- 
tivistic leptons of the jet and the photon field of the high-mass stellar companion. 
We showed then that some variable unidentified EGRET sources in the galactic 
plane could be produced (with luminosities in the range ~ 1 o 35-37 erg s~ x ) by micro- 
quasars with precessing jets. When the jet points towards the observer, gamma-ray 
emission could be detectable yielding a variable source with weak or undetectable 
counterpart at longer wavelengths. The "microblazar" case was studied in some 
detail, reaching the conclusion that some weak or distant sources could periodically 
appear in the detectors like EGRET due to variable Doppler boosting magnification. 

The next step was to add some complexity to the leptonic model, by studying 
also the inverse Compton interactions with other external photon fields due to the 
accretion disk and the corona. We concluded that for HMMQs the predominant 
contribution comes from the scattering of the stellar photon field. This is explained 
by the higher luminosities obtained in this way for MeV gamma-ray energies. For 
higher energies, where the luminosities due to the upscattered disk and coronal 
photons might be important, the absorption by pair creation makes it difficult for 
these gamma-rays to escape in many cases. At this point, the model was applied 
to the recurrent flaring events that have been recently reported for the galactic 
black hole candidate Cygnus X-l at X-ray and soft gamma-ray energies, obtaining 
the observed luminosities and explaining the variability of these events through the 
precession of the jet. 

We then turned towards hadronic models for the gamma-ray emission in HM- 
MQs, where the jets are assumed to have a relativistic proton content. The jets of 
microquasars with massive, early-type stellar companions are exposed to the dense 
matter field of the stellar wind. We presented estimates of the gamma-ray emission 
expected from the jet-wind interaction (typical luminosities of ~ iq34-36 er g g -i^ 
spectra harder than in leptonic cases). The proposed mechanism could explain some 
of the unidentified gamma-ray sources detected by EGRET on the galactic plane 
as well as predict possible neutrino detections since those particles are produced in 
the successive decays that follow the p — p interaction. 

In the case of the subset of unidentified EGRET sources whose spatial distribu- 
tion forms a halo around the galactic center, we considered low-mass microquasars 
as their possible counterparts since these are old sources that, as we saw, can present 
large proper motions. We performed detailed calculations of the jet inverse Comp- 
ton emission in the seed photon fields from the star, the accretion disk, and the 
hot coronal region, in different configurations of parameters such as jet Lorentz fac- 
tors, powers, and angles with the line of sight. The conclusion was that unlike the 
HMMQs case, the external Compton emission largely fails to produce the required 
luminosities. Synchrotron-self-Compton emission appears as a promising alterna- 
tive. 
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In order to test these models in a statistically significant number of sources, new 
generation of gamma-ray detectors is needed. In fact these satellites are already 
planned for the quite near future, which is the case of the AGILE and GLAST 
missions. Their sensitivity is expected to be about 10 to 100 times better than the 
EGRET one. 

The variety of models proposed along the thesis was not only applied to the 
EGRET range of energies, but also to lower ones, like in the case of the soft gamma- 
ray emission from Cyg X-l, as well as to predict TeV gamma-ray production in the 
hadronic or the leptonic cases. The jets of microquasars could not only be the site 
of emission that would cover the electromagnetic spectra from radio to MeV, GeV 
or even TeV energy ranges, but as we proposed here they may also be a source of 
neutrinos. Furthermore, it has even been suggested that they may be important 
sites of particle acceleration in the interstellar medium 1 . 

The existence of TeV emission in MQs might be soon confirmed by new and 
powerful instruments like MAGIC, HESS and VERITAS. Neutrino observations 
with IceCube and ANTARES will be crucial with respect to the hadronic model 
predictions. 



In such a rich context of physical processes to be explained, observations to 
be interpreted, diverse theories to be tested, we certainly foresee new projects to 
straightforwardly continue and immediately complement the work developed during 
this thesis period. To begin with, a deep and detailed study of the SSC emission in 
low-mass microquasars systems is desirable, after realizing the promising prospect 
this process could have to explain the observed luminosities in some cases. Different 
combinations of parameters as well as scenarios with and without adiabatic losses 
will be studied. 

The generation of jets, either steady or transient, is also a pending issue that 
will be explored in the next stage of our research. This is a more difficult ground 
that involves the complex physics relating the accretion disk to energetic outflows 
of matter. The tight correlations observed in MQs between radio emission and 

1 Heinz & Sunyaev (2002) have discussed the possible contribution of X-ray binary jets to the 
production of galactic cosmic rays. They conclude that, whereas in terms of overall energetics 
such jets are still likely to inject less power into the ISM than supernovae, they may contribute 
a specific and detectable component to the cosmic ray spectrum. In particu lar, the shocks 
in the ISM associated with jets from X-ray binaries will be considerably more relativistic than 
those associated with the supernovae, and thus may be considerably more efficient for particle 
acceleration. 
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changes in the X-ray state hopefully will provide some new hints to make progress 
in this direction. 




The central aim of this thesis was to propose diverse models that could support the 
idea that microquasars are possible gamma-ray emitters. In the different models 
developed for this work, gamma-rays are produced by leptonic or hadronic processes 
always related with the interaction of the jets with diverse photon and matter fields. 
The next stage of our research will be connected with the generation of these jets 
as well as with the different lines of research opened along the thesis period. The 
discovery of this local scale-down version of quasars is allowing a better understand- 
ing of the physics related with the accretion- ejection relativistic phenomena in the 
Universe. 
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In this appendix we do not deal with the gamma-ray emission from MQs but 
we rather concentrate on a specific unidentified EGRET source, 3EG J0542+2610, 
that seems originated in a completely different kind of scenario where, nonetheless, 
similar physics to that described in Chapter take place. After a multi- wavelength 
study of the surroundings of this source we show that the only known object within 
the 95% confidence location contour of the source capable of generating the observed 
gamma-ray emission is the Be/X-ray transient A0535+26. We shall argue that the 
gamma-rays are produced during the accretion disk formation and loss phases in 
each orbit, through hadronic interactions between relativistic protons accelerated 
in an electrostatic gap in the pulsar magnetosphere and the matter in the disk. 



A.l Gamma-Ray Source 3EG J0542+2610 

The best estimated position of the gamma-ray source 3EG J0542+2610 is at (I, b) 
~ (182.02,-1.99). Its 95% confidence location contour overlaps with the shell-type 
supernova remnant (SNR) G180.0-1.7 (Romero et al. 1999). The gamma-ray flux 
for the combined EGRET viewing periods is (14.7 ± 3.2) 10 -8 ph cm -2 s -1 . Since 
the flux is highly variable on timescales of months, it should originate in a compact 
object and not in the extended SNR. In Fig. lA.ll we show the EGRET light curve. 
The source switches between periods of clear detections and periods when only up- 
per bounds to the flux can be determined. The variability analysis of Torres et 
al. (2001b) assigns to 3EG J0542+2610 a variability index I = 3.16, which means 
that the source variability level is 4.32a above the average (spurious) variability of 
all known gamma-ray pulsars. Tompkins' (1999) variability index for this source 
(r = 0.7) also indicates that the source is variable. 

If a pulsar origin is discarded due to the high variability and the steep spectral 
index (r = —2.67 ± 0.22), we are left with two main possibilities: 1) the source is 
a background, unnoticed gamma-ray blazar seen through the galactic plane, or 2) 
it is a galactic compact object with an energy budget high enough as to generate 
significant gamma-ray emission and, at the same time, does not have the stable 
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j= 50 
a 




Viewing Period 

Figure A.l: Flux evolution of 3EG J0542+2610 through single viewing periods. 



properties usually associated with isolated pulsars. 

In Fig. IA.2l lower panel, we show a 1.4-GHz VLA map made with data from the 
NVSS Sky Survey (Condon et al. 1998), where all point-like radio sources within 
the 68% confidence contour of 3EG J0542+2610 can be seen. The measured char- 
acteristics of these 29 sources are listed in Table IA.lt Most of them have no entry 
in any existing point source catalog. In those cases where we were able to find 
positional counterparts at other frequencies we have estimated the spectral indices, 
which are also shown in the table. Most of these sources are very weak, at the 
level of a few mJy. No strong (at Jy level), flat or nonthermal source is within the 
location error box of the gamma-ray source. The strongest radio source (No. 23 
in our table) has a rather steep spectrum and is a factor ~ 10 below the minimum 
flux density of firm gamma-ray blazar identifications given by Mattox et al. (1997). 
The nature of this source is not clear at present; it could be a background weak 
radio quasar. The fact that it is not seen at X-rays seems to argue against a galactic 
microquasar or any other kind of accreting source. 

At X-ray energies the most significant source within the EGRET error box is 
the X-ray transient A0535+26, which is discussed in the next section. This source 
does not present significant radio emission and consequently it cannot be seen in 
our maps. We have indicated its position with a star symbol in Fig. IA.21 middle 
panel. We also show in this figure the direction of the proper motion of the system, 
as determined by Lee Clark & Dolan (1999). The upper panel shows the entire 
radio field as determined from radio observations with Effelsberg 100-m single dish 
telescope at 1.408 GHz (data from Reich et al. 1997). The middle panel present an 
enhanced image obtained at 2.695 GHz with the same telescope (Fiirst et al. 1990), 
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Table A.l: Point radio sources within the inner location probability contours of the 
gamma-ray source 3EG J0542+2610 
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where the gamma-ray location probability contours have been superposed (Hart- 
man et al. 1999). We have processed these large-scale images using the background 
filtering techniques described by Combi et al. (1998). 
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Figure A. 2: Radio maps at different resolutions of the field containing the gamma-ray 
source 3EG 30542+2610. Upper panel: 1408-MHz image obtained with Effelsberg 100- 
m single dish telescope. The shell-type SNR G180.0-1.7 is clearly visible in the map. 
Contours are shown in steps of 25 mJy beam" 1 , starting from 20 mJy beam" 1 . Middle 
panel: 2695-MHz map obtained by the same telescope. Contours in steps of 30 mJy 
beam" 1 , starting from 30 mJy beam" 1 . The position of A0535+26 is marked by a star 
symbol. The arrow indicates the direction of the proper motion. EGRET confidence 
location contours are superposed to the radio image. Lower panel: VLA image of the 
inner region at 1.4 GHz. Contours in steps of 1 mJy beam" 1 , starting from 1 mJy beam" 1 . 
The diffuse background emission has been removed from the first two maps. 
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A0536+26 is a Be/X-ray transient where the compact object is a 104s pulsar in 
an eccentric orbit around the BOIII star HDE 245770 (Giovanelli & Sabau Graziati 
1992). Be stars are rapidly rotating objects which eject mass irregularly forming 
gaseous disks on their equatorial planes. If there is a compact companion in a close 
orbit, accretion from the star can result in strong X-ray emission. In the case of 
A0535+26, strong and recurrent X-ray outbursts are observed with a period of 111 
days, which has been identified with the orbital period (Giovanelli & Sabau Graziati 
1992). It is generally agreed that these outbursts occur when the accretion onto 
the neutron star increases at the periastron passage. The average ratio of the X-ray 
luminosity at the periastron to that of the apoastron is ~ 100 (Janot-Pacheco et 
al. 1987). In Table lA~2l we list the main characteristics of the A0535+26 system. 

During a major outburst of A0535+26 in 1994, the BATSE instrument of the 
Compton Gamma Ray Observatory detected a broad quasi-periodic oscillation 
(QPO) in the power spectra of the X-ray flux (Finger et al. 1996). The QPO 
component was detected during 33 days, with a central frequency that was well 
correlated with both the hard X-ray flux and neutron star spin-up rate inferred 
from pulse timing. Finger et al.'s (1996) observations are the first clear evidence 
that an accretion disk is formed during giant outbursts. 

Using the simultaneous variations of the spin and the QPO frequencies in the 
context of the beat frequency model (Alpar & Shaham 1985), Li (1997) has deter- 
mined the evolution of the ratio £ of the inner accretion disk radius to the Alfven 
radius. He found that during the initial rise of the outbursts £ quickly increased 
from ~ 0.6 to ~ 1, indicating a transition of the accretion process from spherical 
accretion before the outbursts to disk accretion during the high X-ray luminosity 
phase. Prior to the direct evidence for a transient accretion disk in A0535+26, 
Motch et al. (1991) had already suggested, on the basis of an analysis of the 
long-term X-ray, UV, and optical history of the system, that two different types of 
interactions exist between the Be star and the pulsar, accordingly with the dynam- 
ical state of the highly variable circumstellar envelope. "Normal" outbursts would 
occur for high equatorial wind velocities (~ 200 km s _1 ), whereas "giant" outburst 
would result from lower equatorial wind velocities (~ 20 — 80 km s _1 ) which allow 
the formation of a transient accretion disk. 

It is interesting to notice that, during the initial stage of the 1994 outburst, 
when QPOs were found in A0535+26, the gamma-ray source 3EG J0542+2610 was 
detected by EGRET with a flux of (39.1 ± 12.5) 10~ 8 ph cm -2 s _1 (viewing period 
321.1: February 8 - 15, 1994). But when A0535+26 was at the peak of its X-ray 
luminosity, on February 18, the gamma-ray source was not detected (viewing pe- 
riod 321.5). The gamma-ray emission seems to have been quenched precisely when 
the accretion disk was well-formed and maximally rotating. In the next section 
we present a model that can account for the gamma-ray production in A0535+26 
necessary to explain the EGRET source 3EG J0542+2610 and that is in agreement 
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Table A. 2: Physical parameters for A0535+26 (from Janot-Pacheco et al. 1987 and 
Giovanelli & Sabau Graziati 1992) 



Primary spectral type 
Primary mass 
Secondary mass 
Primary mass loss rate 
Distance 

L x peak (average) 
X-ray pulse period 
Orbital period 
Orbital eccentricity 

^max / j^min 



BOIII 

9 - 17 M Q 

< 2.7 M 

7.7 1(T 7 M G yr- 1 

2.6 ± 0.4 kpc 

7.5 ± 2.4 10 36 erg s" 1 

104 s 

111 ±0.5 d 
0.3 - 0.8 
~ 100 



with our present knowledge of the Be/X-ray transient source. 



A. 3 The Model 

Our purpose in this section is to show that there exist a plausible mechanism that 
could explain the gamma-ray emission of the source 3EG J0542+2610 as originated 
in the X-ray transient A0535+26. This mechanism should be capable of predict- 
ing the observed gamma-ray flux, the variability in the lightcurve, the fact that 
no gamma-ray emission was observed on February 18, 1994, when A0535+26 was 
at the peak of the X-ray outburst, but also that it was positively detected in the 
previous days when the X-ray flux was rising, and finally the fact that A0535+26 is 
not a non-thermal radio source. This latter restriction seems to suggest a hadronic 
origin for the gamma rays. Otherwise, relativistic electrons should also produce 
synchrotron radio emission. 

Cheng & Ruderman (1989, 1991) have studied the disturbances produced in 
the magnetosphere of an accreting pulsar when the Keplerian disk rotates more 
rapidly than the star. As in the case of equal angular velocities, when f2* < f^, 
inertial effects of electrons (— ) and ions (+) lead to a complete charge separation 
around the "null surface" f2* • B = 0. The equatorial plasma between the inner 
accretion disk radius ro and the Alfven radius ta co-rotates with the disk, whereas 
the rest of plasma co-rotates with the star. An electrostatic gap with no charge at 
all is then created around the "null surface" (see Fig. IA.3|) . In this gap E • B ^ 
and a strong potential drop is established (see Cheng & Ruderman 1991 for details). 

Cheng & Ruderman (1989) have shown that the potential drop along the mag- 
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Figure A. 3: Sketch of the magneto sphere model for A0535+26 when f2d > The gray 
region co-rotates with the star, whereas the hatched region co-rotates with the accretion 
disk. An empty plasma gap where E • B ^ exists around fi* • B = 0. Protons from the 
stellar surface are accelerated in this gap and collide with the disk, where they produce 
7T° -decay gamma-rays. Adapted from Cheng & Ruderman (1991). 



netic field lines through the gap is: 

B s R s n d (r ) 



AV m 



r c 

1/7 



~ 4 10l4r5/2 R " /7L ^ B ^ /7 V ' ^ 

where B s is the neutron star's surface dipole magnetic field, R is its radius, r is 
the inner accretion disk radius, M is the compact star mass, and L37 is the X-ray 
luminosity in units of 10 37 erg s _1 . The radius and magnetic field in the second 
expression are in units of 10 6 cm and 10 12 Gauss, respectively. The parameter 
(3 = 2ro/rA is twice the ratio of the inner accretion disk radius to the Alfven radius. 
For the particular case of A0535+26 we can adopt j3 ~ 1, according to the estimates 
by Li (1997). 

Protons entering into the gap are accelerated up to energies above E p ~ eAV mSuX 
whereas inverse Compton and curvature losses would limit the energy gain of elec- 
trons and positrons to lower values. The maximum current that can flow through 
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the gap can be determined from the requirement that the azimuthal magnetic field 
induced by the current cannot exceed that of the initial field B (Cheng & Ruderman 
1989): 



3^-2 



4ax ~ cB s R r Q 



1.5 1 24 /T 2 ( ^) ~ 2/ k /7 Li?B^ 7 esu s~\ (A.2) 



As it can be seen from Fig. IA.31 the proton current flux is directed from the po- 
lar cap of the star, where the accreting material impacts producing strong X-ray 
emission and abundant ions to be captured by the gap, towards the accretion disk. 
Electrons move in the opposite sense in order to cancel any net charge current flow, 
keeping in this way E • B ^ in the gap. 

The collision of the relativistic proton beam into the disk will produce hadronic 
interactions with copious 7r° production (Cheng et al. 1990). These tt° will quickly 
decay into gamma-rays that could escape only if the disk density is sufficiently low. 
Otherwise, they will be absorbed by the matter and then re-emitted as X-rays. The 
interaction of relativistic protons with a thin hydrogen layer has been studied by 
Cheng et al. (1990) in the context of their model for the Crab pulsar. The 7r°-decay 
gamma-rays will only escape insofar as the column density of the disk would be 
X < 100 g cm -2 . From Novikov & Thorne (1973), we have: 



7 IOV- 3 / 5 ' Z '\-^LT"B-^ g cm" 2 . (A.3) 



In a system like A0535+26, the column density of the disk will evolve with time 
following the variations in the X-ray luminosity. All other physical parameters in 
Eq. (jA.3|) remain constant along the orbital period. If at some point near the pe- 
riastron, when the X-ray luminosity is close to its maximum, the column density 
exceeds the value ~ 100 g cm -2 , the medium will be no longer transparent to 
gamma-rays and the gamma-emission will be quenched (we are assuming that the 
disk is seen from the opposite side to that where the current impacts). This could 
explain the fact that A0535+26 has not been detected by EGRET at the peak of 
its X-ray outbursts on February 18, 1994. The detection, instead, was clear a few 
days before, during the previous viewing period. 

The expected hadronic gamma-ray flux from A0535+26 on Earth will be: 

F(F\ 1 r fMM HF (A A) 

F{E ^-m^ L dtdE^o dE ^ (A - 4) 
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where D is the distance to the source, Af2 is the beaming solid angle and: 

d 2 N p (E n o) _ 2 J max S < m n o > da(E p , E n o) 
dt dE^o e m p dE n o 

is the differential production rate of neutral pions from p — p interactions. In this 
latter expression, da(E p , E n o)/dE n Q is the differential cross-section for the produc- 
tion of 7r°-mesons of energy E n o by a proton of energy E p in a p — p collision, and 
< m n o > is the mean multiplicity of 7r°. 

For calculation purposes we shall adopt the cross section given by Dermer 
(1986a, b) and a beaming factor AQ/Att « 0.3 as used by Cheng et al. (1991). 
We shall evaluate the total gamma-ray luminosity between 100 MeV and 20 GeV of 
A0535+2610 at an early epoch t of the disk formation, when the X-ray luminosity 
ratio is L^ ax / ' L l x ~ 10. If the absorption feature observed near 110 keV is a cy- 
clotron line, the polar magnetic field of the neutron star results 9.5 10 12 G (Finger 
et al. 1996). In our calculations we adopt this value along with /3 ~ 1, R§ — 1 
and M = 1.4 M & (Janot-Pacheco et al. 1987). We then obtain that the column 
density of the disk is ~ 42.5 g cm -2 at this stage and, consequently, gamma-rays 
are not absorbed in the disk material. The potential drop in the electrostatic gap 
results AVmax ~ 2.5 10 13 V whereas the proton current deposited into the disk is 
N p = J m ax/e ~ 4.8 10 32 s -1 . Using Eq. (fOjl for photons in EGRET's energy range 
we obtain: 

F(100 MeV < E 7 < 20 GeV) « 5.1 10~ 8 ph cm" 2 s~ x . (A.6) 

This flux is consistent with the lower EGRET detections. As the accretion onto 
the neutron star increases the gamma-ray flux also increases, reaching a maximum 
when S w 100 g cm -2 . At this point the flux has grown by a factor ~ 4 (notice the 
dependence on J max in addition to E) and photon absorption into the disk becomes 
important quenching the radiation. The gamma-ray source, then, is not detected 
when the peak of the X-ray luminosity occurs. The additional gamma-ray flux ex- 
pected from secondary electrons and positrons radiating in the the strong magnetic 
fields anchored in the disk can explain even higher EGRET fluxes (Cheng et al. 
1991). A more detailed analysis in this sense will be presented elsewhere, including 
spectral considerations. The observed spectrum depends on the proton injection 
spectrum at the gap. Strong shocks near the polar gap can accelerate the pro- 
tons up to a power law which will be preserved through the propagation across the 
gap and imprinted into the 7r° gamma-ray spectrum. Additionally, e + e~ cascades 
induced in the disk material can result in a relativistic Bremsstrahlung contribu- 
tion. Electromagnetic shower calculations are in progress in order to determine this 
contribution (Romero et al., in preparation). 
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A. 4 Discussion 

The model here outlined does not imply strict periodicity because of the highly 
chaotic nature of the Be stellar winds, which can significantly vary on short time- 
scales producing strong changes in the accretion rate. Notwithstanding, the general 
prediction that the peak of the gamma-ray emission should not be coincident with 
the maximum X-ray luminosity during a given outburst can be used to test the gen- 
eral scenario proposed here. At present, the poor time resolution of the gamma-ray 
lightcurve does not allow correlation studies. EGRET data for the best sampled 
X-ray outburst (in February 1994) consist only of two viewing periods, as it was 
mentioned, and in one of them the source was not detected. In the future, however, 
new instruments like GLAST could provide the tools for these kind of investigations. 

TeV emission should be produced in the accretion disk according to Eq. (jA.4|) . 
although degradation effects during the propagation in the strong magnetic and 
photon fields around the accretion disk could suppress much of it. The magnetic 
field at the inner accretion disk radius is (e.g. Cheng et al. 1991): 



B{r ) = 3 10 5 f3- 3 B u 5/7 R~ 9/7 L% 7 {M/M Q )- 3/7 G (A.7) 

which for the typical values of A0535+26 yields B ~ 10 4 G. TeV gamma-rays can 
be absorbed in the field through one-photon pair production above the threshold 
given by (e.g. Bednarek 1993): 

* = 2R = 77> ( A - 8 ) 
2m e c z B cv 15 

where 9 is the angle between the photon momentum and the magnetic field and 
B CT is the critical magnetic field given by B CT = m 2 c 3 /eh ~ 4.4 10 13 G. This pro- 
cess, then, suppress gamma-ray photons with energies higher than ~ 300 TeV in 
A0535+26. Since this value is well above the energy of the protons that impact on 
the disk, we find that absorption in the magnetic field should not occur. 

However, TeV photons with lower energies should be absorbed by two-photon 
pair production in the accretion disk X-ray photosphere. Calculations by Bednarek 
(1993) show that the optical depth quickly goes to values above 1 for disks with 
luminosities L37 ~ 1. The opacity effects can reach even GeV energies leading to 
a steepening in the spectrum respect to what is expected from a pure pion-decay 
mechanism. The fact that the observed spectrum in 3EG J0542+4610 has an in- 
dex T ~ —2.7 seems to support the idea that important absorption is occurring in 
the X-ray photosphere of this source. Additional spectral modifications should be 
produced by the emission of secondary pairs in the magnetic field close to the disk 
(Cheng et al. 1991, Romero et al., in preparation). 



ENGLISH A. 5 CONCLUSIONS 



In their original model, Cheng & Ruderman (1989) suggested that if the disk 
is sufficiently dense then the gamma-rays could be produced only in a moving 
low-density "window" . This window would collimate a pencil beam of gamma-rays 
aligned with the magnetic axis. Consequently, a pulsed emission could exist with the 
same period of the X-ray source (104 sec in the case of A0535+26). Small changes in 
the gamma-ray period might be produced by the radial motion of the "window" . It 
would be interesting to test whether this pulses are present in A0535+26. However, 
since the disk is a transient structure and only would excess the critical density 
during a few days per orbit, the number of photon counts in EGRET data are too 
low to allow a periodicity analysis as in the case of isolated and stable gamma-ray 
pulsars. Instruments with higher sensitivity like GLAST could sum up over several 
viewing periods in order to look for these features. 



A. 5 Conclusions 

We have shown that the Be/X-ray transient system A0535+26 can be also a tran- 
sient gamma-ray source under very reasonable assumptions. The existence of a 
QPO phenomenon detected by BATSE during the 1994 X-ray outburst provided 
direct evidence of the formation of a transient accretion disk near the periastron 
passage. In the beat frequency model for QPO the Keplerian orbital frequency of 
the material at the inner edge of the accretion disk should exceed the spin frequency 
of the neutron star. Cheng & Ruderman (1989) have shown that in such a circum- 
stance an electrostatic gap is open in the magnetosphere around the null surface 
determined by fi* • B = 0. This gap can accelerate protons up to energies of tens of 
TeV, producing a hadronic current that impacts into the accretion disk generating 
7T° gamma-rays. The transient character of the disk makes the high-energy gamma 
radiation highly variable, as was found in the observed gamma-ray flux evolution 
of 3EG J0542+2610. 

A specific prediction of the model is the suppression of gamma-ray emission 
when the column density of the disk exceeds a critical value, near the peak of the 
X-ray luminosity. Future GeV and TeV observations of this source with instruments 
of high temporal resolution, like GLAST or 5@5 (see Aharonian et al. 2001), could 
be used to test the proposed model and, if it is basically correct, to probe the 
evolution of the matter content on the accretion disk in this extraordinary X-ray 
binary. 
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The aim of this appendix is to provide some insight into the conceptual tools 
used in gamma-ray astrophysics and in particular in this thesis. We will summarize 
the mechanisms that produce gamma rays in the astronomical events that were 
discussed along this thesis, and how these gamma rays can be affected once they 
are created. We will closely follow the presentation made by K.S. Cheng & G.E. 
Romero in their recent book Cosmic Gamma-Ray Sources (Cheng & Romero 2004). 



B.l Basic Concepts 

It is convenient to start by introducing as a basic concept the number of particles 
incident per unit of surface area per unit of solid angle per unit of time arriving 
at a given, unspecified detector. We will call this basic quantity the intensity of 
particles, and we will denote it by I. In general, we will use a subscript to indicate 
the type of particles, e.g. 1 1 and I p denote intensities of gamma rays and protons, 
respectively. 

Once the intensity is introduced, we can define the particle flux as 



where the angle 6 is determined by the direction of motion of the particles with 
respect to the normal to the area, and the integration is performed over the solid 
angle. For isotropic radiation the flux is F = irl, and the number of particles per 
unit of volume is 



In most cases v = c is a good approximation because we deal with relativistic par- 
ticles. 




(B.l) 




(B.2) 
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Normally, we will have particles with different energies, so it is useful to introduce 
a particle energy distribution N(E) such that 

poo 

N= N(E)dE. (B.3) 



The number density of particles with energies greater than E is obtained just by 
integrating from E. In a similar way, the integrated flux density is 

poo 

F(> E)= F(E)dE. (B.4) 
Je 

The luminosity of a source located at a distance d that radiates isotropically is given 
by 

poo 

L(> E) = And 2 I F(E)dE, (B.5) 
Je 

where d is the distance to the source. 
The energy density of the particles is 



w 



/ EN(E)dE. (B.6) 

•J -Emm 



The energy flux is obtained from this expression just by multiplying by c/An, if we 
deal with relativistic particles or photons. 

Let us consider now that a flux of particles of type a, with velocity v a , interacts 
with some target formed by particles of type b within a volume dV. The number 
of particle interactions of a given type, dN iy occurring in a time dt in the volume 
dV will be proportional to the number of particles b in the volume dV and to the 
number of incident particles that traverse the cross section dA of that volume in 
the time dt: 

dN i = da i (n° b dV)(n a v a dt). (B.7) 

In this expression n° and n a are the densities of target and incident particles in a 
coordinate system with the target at rest. The differential cross section dai charac- 
terizes the number of reactions of type % occurring per unit of time in unit volume for 
a unit flux density of incident particles and unit density of the target. It is measured 
in units of area, the standard unit being the mb (i.e. 10~ 3 barn, 1 barn= 10 -24 cm 2 ). 

The total cross section for a given interaction a { is the sum over all possible 
momenta of the resulting particles after the interaction. Both <jj and dai are rel- 
ativistic invariants. The total cross section <7 to t is obtained by summing the cross 
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sections of all possible processes that occur upon the interaction of particles a and 
b, i.e. 



Ctot 



J>. (B.8) 



The relative probability of a given reaction channel is given simply by p = <Ji/a tot . 

In the case of gamma-ray emission, if the generation of the gamma rays is due to 
the interaction of particles of type % with a given intensity Ii(Ei, r) with a target of 
density n(r), we can write the intensity of the radiation from the resulting gamma- 
ray source as 

J 7 (£ 7 )= / / ntfaiE^EjIiiEirfdEidr, (B.9) 

Jl JE-f 

where / defines the direction along the line of sight (i.e. I = f/r). The emissivity of 
the gamma-ray source is defined as 

POO 

g 7 (£ 7 ,f) / ntfaiE^EjIiiEirfdEi, (B.10) 
Jej 

in such a way that 

J 7 (£ 7 ) = J^(E 7 ,r)dr. (B.ll) 



B.2 Gamma— Ray Production 
B.2.1 Thermal Mechanism 

Any body with a temperature different from zero emits thermal radiation. For a 
perfect absorber in thermodynamical equilibrium at temperature T (i.e. for a black 
body) the spectrum will be given by the Planck formula: 



Ibb(E ] 



2El 



ph 



(he) 



exp(E ph /kT) - 1 



(B.12) 



where h and k are the Planck and Boltzmann's constants, respectively. The corre- 
sponding number density of photons per unit energy is: 



riBB{E } 



ph; 



(E ph /m e c 2 ) 2 
m e c 2 (7r 2 A|?) 



1 



exp(E ph /kT) - 1 



;b.i3) 
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In this last expression, A c = (27im c c)~ 1 h is the Compton wavelength of the electron. 
The maximum of the intensity occurs at 

£ph, m ax~4.7xl(T lo (-0 MeV, (B.14) 

and the average energy of the photons is 

(E ph ) = 2.7kT « 2.3 x 1(T 10 (J^J MeV. (B.15) 

In order to have photons with average energies of 1 GeV, temperatures of ~ 10 13 
K are necessary. These temperatures cannot be found in steady astrophysical ob- 
jects, but only in explosive events and in the Big Bang. In addition, the photon 
density in a source with T ~ 10 13 K would be ~ 3 x 10 34 cm~ 3 . Since the mean 
free path of a photon in this radiation field is A 7 ~ (ncr 77 ) _1 << 1 cm, the source 
would be self-absorbed by photon-photon pair creation. 

Typical astrophysical gamma-ray sources in the continuum are non-thermal 
sources where the gamma rays are produced by the interaction of relativistic par- 
ticles with radiation or matter fields. We describe next the interactions that are 
present in our models. 



B.2.2 Non-Thermal Mechanism I: Particle— Field Interac- 
tions 

SYNCHROTRON RADIATION 

A relativistic particle moving in a magnetic field will emit photons within an angle 
6 ~ mc 2 / E of its direction of motion. In a magnetic field B an electron moves 
along a helical path with an angular frequency Ub given by 

m e c hi 

The radiation spectrum of the electron is given by (e.g. Ginzburg & Syrovatskii 
1964): 

P{E) = %-zB ± - / K B/a { V )dn, (B.17) 
m e c cj c Je/e c 



where E = hv is the energy of the radiation, B ± = B sin 6, 6 is the pitch angle, and 
K 5 / 3 is a modified Bessel function of the second kind. The characteristic energy of 
the photons is given by 
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3/i eB 1 / £ x 2 



^c = t j • (B.18) 

47r m e c \m e r 

The maximum of P(E) occurs at E max = 1.9 x lO~ n B±(E/GeV) 2 GeV. We see 
then that only for extremely energetic particles and strong magnetic fields we can 
get gamma-ray photons from synchrotron radiation. 

The total energy rate loss by synchrotron radiation of an electron moving in a 
field B can be obtained by integrating equation (jB.17|) . The result is 



(IF., \ 2 / e 2 



dt J 3 V vrifC 2 

' syn x 



BiY, (B.19) 



where 7 = E e /m e c 2 is the Lorentz factor of the particle. 

Introducing the Thomson cross section <tt = 87re 4 /3m 2 c 4 ~ 0.665 x 10~ 24 cm -2 
and averaging over an isotropic pitch angle distribution, the expression for the 
energy losses can be set in the following convenient form: 

<U " \(J T cw m ^ 2 = 0.66 x 10 3 5 2 7 2 eV/s, (B.20) 



dt Ayn 3 

where w m3ig = B 2 /8ir is the magnetic energy density and B is measured in Gauss. 
If we have a homogeneous and isotropic power-law electron distribution given 

by 

N e (E e )dE e = K e E' p dE e (B.21) 
in a random magnetic field, the resulting spectrum is (Ginzburg & Syrovatskii 1964): 

,3 / 3e x(p-D/2 



^ 7 ^ ^m e c 2 \Airmlc 



B^ 2 K e LE~^ 2 . (B.22) 



In this expression L is the characteristic size of the emitting region and a(p) is given 
by 



2(p-D/2 ygr (M) r (^g) r (zf) 
80F(p + i)r 



«(P) = ; r-) 12 ,\^/ p 2 +7 { V 4 - (B.23) 



a(p) is 0.147, 0.103, 0.0852, and 0.0742 for p = 1.5, 2, 2.5, and 3, respectively. The 
emission is a power law with index a = (p — l)/2. 



In the case of a homogeneous magnetic field, the degree of polarization is 
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D.W = ^j, (B-24) 

which typically yields values in the range 69-75 % for p between 2 and 3. If the 
magnetic field has random component B T , then the degree of polarization will be 

n(p) = n ^^, (B.25) 

where _B r is the random component. Additional details can be found in Ginzburg 
& Syrovatskii (1964, 1965), in Blumenthal & Gould (1970) and in Longair (1992, 
1997). 



INVERSE COMPTON (IC) INTERACTIONS 

The scattering of relativistic electrons on soft photons can produce gamma rays. 
According to equation (|B.9|) the intensity of the radiation from this process when 
the soft photon field has a density n^E^r) is 

p roc /*oo 

lf(E,)= / I e (E e ,r)a(E e ,E^E ph )n ph (E ph ,f)dE e dE ph dr. (B.26) 

Jl JEry JO 

We can introduce a parameter £ = E e E ph /(m e c 2 ) 2 , such that for £ << 1 the 
scattering is classical. In such a limit the cross section can be approximated by the 
Thomson cross section <tt and the average energy of the emerging photons will be 

(E y ) = ^(E ph ) 7 \ (B.27) 

where (-Eph) is the average energy of the target photons. The energy losses for an 
electron in a photon field of energy density u> p h when £ << 1 can be approximated 
as (e.g. Ginzburg & Syrovatskii 1964): 

dE \ 

M - ca T w ph 7 2 w 2 x KT 14 u; ph 7 2 eV/s. (B.28) 



dt 



ic 



Comparing with equation (|B.20J) . we see that at these energies the ratio of syn- 
chrotron to IC cooling times is simply ti C /t syn w w mag /w p h. 

If the incident electron spectrum is a power law and the photon field can be 
approximated by a monoenergetic distribution, then we get from equation (|B.26j) : 

I™(E 7 ) = \n pY MT{rn e c 2 f-v Q (E ph )\ " V K e E-^' 2 . (B.29) 
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Here, L is the typical source dimension, and (E ph ) and n ph are average values for 
the photon energy and the photon density in the source (Ginzburg & Syrovatskii 
1964). If the photon field is thermal radiation, then we get from equation (jB.26|) : 



I?(&r) = (rg/47rW) LK e (kT)W*F(p)E: 



(P+l)/2 



(B.30) 



where T is the temperature and 



F(p) 



2 p+ V + 4p + ll)r [\{p + 5)] C [\ip + 5 
(p + 3) 2 (p + l)(p + 5) 



(B.3i; 



where ( is the Riemann function. For p = 1.5, 2.0 and 2.5, F(p) = 3.91, 5.25, and 
7.57 respectively (Blumenthal & Gould 1970). 

If £ >> 1, then the electron gives most of its energy to the photon: E 1 ~ E e . 
The cross section at these energies decreases drastically and can be represented by 
the well-known Klein-Nishina formula (e.g. Heitler 1954): 



3 m e c 2 
&KN{E e , E ph ) = g°T^ — yln 



2jE, 



m P c z 



(B.32) 



Here, 7 = E e /m e c 2 is the electron Lorentz factor, as usual. The electron energy 
losses are now given by 



dE e 
~dt 



KN 



IO 



3 (m e c 2 \\ f2E e (E ph ) 1 

.2 \ 2 



io- 14 w ph 



A useful general expression for IC scattering by an electron moving in a mo- 
noenergetic, isotropic photon field has been obtained by Jones (1968), applying the 
so- called "head-on" approximation: the seed photons are treated as coming from 
the direction opposite to the electron velocity. The corresponding cross section can 
be written as (e.g. Blumenthal & Gould 1970): 



^ic(£,e P h,7) = A 3<TT 2 f( x ) 
4e ph 7 2 



(B.34) 



where 



/(*) 



2x In x + x + 1 — 2x + 



(4e ph 7x) 2 (l - x) 
2(l + 4e ph7 x) 



P(l/4 7 2 ,l,a;), 



(B.35) 
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being e p h = E p ^/m e c 2 the target photon energy and x a function of the energy 
e 7 = E 1 /m e c 2 of the scattered photons given by 

4e ph 7 2 (l - e 7 /7j 

The function P in equation (jB.35|) is 1 for I/47 2 < x < 1 and otherwise. It con- 
strains the cross section to the physical case, where the energy of the scattered pho- 
tons cannot be lower than that of the seed photons or higher than the energy of elec- 
tron. The maximum energy of the scattered photons is ~ (4e p h7max)/(l + 4e P h7max), 
with 7 m ax the maximum Lorentz factor for the electrons. 



B.2.3 Non-Thermal Mechanism II: Particle— Matter Inter- 
actions 

HADRONIC GAMMA RAY EMISSION 
7T°-decays from proton-proton interactions 

The dominant 7r-producing channels in hadronic interactions are (E t hi ~ m^-c 2 ): 

p + p -> p + p + air + b(n + + tt~) (B.37) 
p + p -» p + n + ir + + avr° + b(ir + + vr~) (B.38) 
p + p -» n + n + 2-K + + a7r° + 6(tt + + tt~), (B.39) 

where a and h are positive integers. Neutral pions decay into gamma rays with a 
proper lifetime of only 9 x 10~ 17 s. The gamma-ray emissivity generated through 
such decays at a source with a proton spectrum I p (E p ) = (c/4ii)N p (E p ) given by a 
power law: 

I P (E P ) = KEf (B.40) 



is 



*r(^) = 2 f % m dK, (B.41) 



where 



2 4 

mzc 



ET(EJ = E, + - r V, (B.42) 



4£ 7 



and 
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da*(E w , E PJ 

v 



<h(E K ) = 4vr / I P (E P ) nK v> dE p . (B.43) 

?min CLU/jr 



Here, da n (E n , E p )/dE 7T is the differential cross section for the production of 7r°- 
mesons of energy E n by a proton of energy E p in a p — p collision. 

The gamma-ray emissivity can be approximated at high energies by 

2Z r 

g 7 (£ 7 ,f) « 4^0-^(7=)— p^/ P (S 7 )77, (B.44) 

where Z p _^a are the so-called spectrum weighted moments (Gaisser 1990, Drury et 
al. 1994) and is ~ 0.275, 0.26, 0.245 for Y ~ 1.6, 1.8, 2.0 respectively, rj ~ 1.5 is 
a parameter that takes into account the contribution of nuclei other than protons 
into the gamma-ray production, and the cross section for inelastic p — p interactions 
can be approximated at E p « 10E^, with E^ > 1 GeV, by 



a pp (E p ) fa 30 x 



0.95 + 0.06 log 



GeV 



mb. (B.45) 



The spectral gamma-ray intensity (photons per unit of time per unit of energy- 
band) is 

J 7 (£ 7 ) = / g 7 (E 7 ,f)rf 3 r, (B.46) 
Jv 

where V is the interaction volume. The photon flux observed at the Earth from a 
source at a distance d then results 

= (B-47) 

Since the injection proton spectrum was a power law, we can also expect a power 
law spectrum for the observed gamma rays: 

F 7 (£ 7 ) oc E~ r . (B.48) 

For energies lower than E p ~ 10 GeV (i.e. which corresponds to E 1 ~ 1 GeV) 
this simple approach is not possible and we must use more complex parameteriza- 
tions of the cross section (see Blattnig et al. 2000). 

As indicated by the reaction chains ()B.37j) - (|B.39|) . electron-positron pairs from 
the decay of the charged pions are also created in p — p interactions. These pairs 
can in turn generate gamma rays through relativistic Bremsstrahlung, IC, and/or 



156 



B RELEVANT RADIATIVE PROCESSES 



synchrotron processes depending on the environment. The total number of pi- 
ons generated in a p — p interaction depends on the energy: N W (E P ) = N E^, 
with one third each of ir°, tt + , and 7r~. The total energy that goes to pions is 
£jtot _ (E n ) = kE p , with k = k Ep. Typical values for the parameters are (e.g. 
Ginzburg & Syrovatskii 1964): iVo = 3.3, 5 = 1/4, fc = 1/3, and a — 0. Hence, 
(E n ) = kE/N 7T (E p ) pa O.lEp^ 4 . If we have an incident proton spectrum with a 
power law intensity I p (E p ) = K P E~ T interacting with a medium of particle density 
n, the pion emissivity can be approximated by 

q^E n )dE n = a pp nN E 5 p I{Ep)dE p = a pp nK; T *dE n , (B.49) 

where 

r - = < B50 > 

1 + a — o 

and 

rv-i 



' (R51) 



The decays ir^ ^ ^ + v and ^ — > e ± + v + v will lead to the injection of 
leptons and neutrinos in the source. The electron emissivity can be approximated by 
(Ginzburg & Syrovatskii 1964, see also Dermer 1986c for a more general expression): 



q e {E e )dE e = o vv nK e {T)K^E- v dE e , (B.52) 

where 

K.<?)=(<) r ~ 1 , 2(F + 5) v (B.53) 

1 ; \mj r(r + 2)(r + 3) 1 ; 



Here, m^c 2 pa 140 MeV and m^c 2 pa 106 MeV. The mean energies of the produced 
secondaries in the lab frame are (E e ±) = E T ±/A, (E v ^) = E w ±/A, and (E 7 ) = E w o/2. 
Since we have equal luminosities for the different pions, then L 7 ~ L v>p ~ L n /2. 



ELECTRON POSITRON ANNIHILATION 

Electron-positron annihilation can be an important source of gamma rays through 
the reaction e + e + — > 7 + 7. When the two particles are at rest the energy of the 
resulting photons is simply E 1 = m e c 2 = 0.511 MeV. Line radiation at this energy 
is usually referred as annihilation radiation. If one of the leptons moves at high 
velocity when it collides with the other at rest, then one of the photons will have a 
high energy whereas the other will have an energy ~ 0.511 MeV. The cross section 
for pair annihilation of an electron with energy r ym e c 2 with a positron at rest is 
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(B.55) 
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Expressed in terms of the center-of-mass frame (cm) this becomes 
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(B.57) 



where all quantities are referred to the center-of-mass system of the colliding parti- 
cles. 

Interpolating between the non-relativistic and ultra-relativistic limits of the cor- 
responding pair annihilation rates, a simple expression with an accuracy within 14 
per cent can be found (Coppi & Blandford 1990): 

R e± « |^[lna; + x- 1/2 ], (B.58) 

where x = 7 e +7 e -. Then, the annihilation luminosity for a plasma with particle 
density given by N e ±(E e ±)dE e ± is: 

L (ann) = J ( ^ + + ^_ } ^ jV e+ (£ e+ ) jV e - (£ e -) dE e+ dE e -dV. (B.59) 

Electron-positron annihilation can also occur with the emission of a single pho- 
ton, but in this case the electron must be bound to an atom. The one-photon 
annihilation cross section when the electron is bound to an atom with charge eZ is 



lph 4vrZ 5 a 4 r 2 



2 1 2 7 4 7 + 2 
/ + 3 + 3 ^7 



ln[(l + /3) 7 ] 



(B.60) 



£7(7+ I) 2 

AnZW ° for 7 »l (B.61) 

7 

— ^ for /3 « 1, (B.62) 
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where as before a is the fine structure constant and the energy of the positron is 

7m e c 2 . 

Three or more photons can be produced in the annihilation of free electron- 
positron pairs, but the cross section is down by a factor ~ a % ~ 2 ~ (1/137)*" 2 from 
the free two photon annihilation case, where i is the number of resulting photons. 



B.3 Gamma-ray absorption processes 

Gamma rays, once created at the source can be absorbed by photon and matter 
fields either in the same source or in the medium between the source and the de- 
tector. We worked along the thesis with one of the main mechanisms that results 
in the absorption of gamma rays: pair creation in photon-photon interaction. 

If gamma rays with an initial intensity I®(E y ) are injected into a medium of 
particle density n, the intensity after traversing a distance x will be 

J 7 (£ 7 ) = /?(^ 7 )e- T , (B.63) 

where 

r = anx (B.64) 

is the optical depth and a is the cross section of the relevant interaction for the 
photons. The quantity A 7 = (cm) -1 is the mean free path of the photon in the 
medium. The probability for a photon to interact after traveling a distance L is 

1 - e- L /\ 



B.3.1 Phot on- Phot on pair creation 

A gamma-ray photon of energy E Jt i can produce an electron-positron pair in a 
collision with a photon of energy E lt2 if E 1 ^E 1 ^ 2 > ijn e c 2 ) 2 . The pair creation cross 
section is (Dirac 1930): 



(7 77 (£' 7) i, E Jt 2) 



P 2 ) 



2(3((3 2 - 2) + (3(3 4 ) In 



1 + P 
l-(3 



(B.65) 



where r is as before the classical radius of the electron and 

,2\2 1 1/2 



(m e c 



(B.66) 



The outgoing electron (positron) has an energy m e c 2 / >/l[3 2 in the center-of-mass 
system. It is also useful to express the cross section in terms of the total energy 
square in the center of mass frame (s = (p 7j i + p^^) 2 )'- 
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where p 7 is the four-momentum of photon and 
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4 (m P c 2 ) 2 
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(B.68) 



The threshold of pair creation is s=4?7igC 4 . 

If a gamma ray of energy E 1 should traverse a region of size R with a photon 
field of number density N(E ph ,r)dE ph , the optical depth is given by: 



'^7ph(-E'7) 



oo pR 



JO 



N(E ph , r)(x 77 (£ 7 , E ph )dr dE ph . 



(B.69) 



Gould & Schreder (1967) present estimates of the absorption by a blackbody pho- 
ton gas and a power-law photon spectrum. 



In the case of a gamma-ray source with an intrinsic luminosity L 7 , the intrinsic 
77-attenuation will be determined by the compactness parameter I, which is defined 
as the ratio of the intrinsic luminosity to the source radius R, and the mean photon 
energy (E y ): 

r 77 « a 77 n 7 i? = 

For (E^) = 1 MeV, we have r ~ 1.7 x 10 -31 /, where / is measured in erg s^ 1 cm -1 . 
We see, then, that compact and luminous gamma ray sources can be self-absorbed 
through pair production. If the radiation is beamed, the inferred (e.g. through 
variability observations) luminosity might be significantly larger than the intrinsic 
value (£ app = D n L int , with 3 < n < 4, where D is the Doppler factor 1 ). In such a 
case the opacity constraints are not so strong. 



(B.70) 



1 The Doppler factor is defined as D = [r(l — /Jcos^)] -1 ), where T is the bulk Lorentz factor 
of the radiating plasma, is its velocity in units of c and (f> is the angle between the beamed 
radiation and the line of sight. 
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Many explanations have been proposed to understand the existence of these ap- 
parently superluminal, v > c, motions in radio sources. The most reliable is that 
apparent superluminal velocities are attributable to bulk relativistic motion along 
the line of sight to the continuum source. We will expand this idea following the 
model presented in Blandford et al. (1977). We shall use Fig. IC.ll to a better 
understanding of the situation. 



B 




Figure C.l: The geometry assumed in the model. The observer at A sees a radio source 
moving from B to B' at a speed apparently exceeding the speed of light, c. This effect can 
happen if v — > c and the angle to the line of sight, 9, is small, but non-zero. 

Let us consider that the source is at point B at time t%. Its emission is detected 
by an observer at point A at time t[ after the light has crossed the distance AB. 
Later, at t 2 = t\ + St, the source has moved a distance vSt. Observations are made 
again, recorded at t' 2 , taking into account the light-travel time between the observer 
and the new position of the source. The distance between A and B is d + v St cos 6 
while the distance between A and B' can be approximated by d if the angle to the 
line of sight, 9, is small. 

The angular separation between B and B' is given by, 
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Figure C.2: The apparent transverse velocity /3t = vt/c of a source moving at an angle 
9 to the observer's line of sight as a function of the Lorentz factor T. It can be seen how, 
for r>l, the apparent transverse velocity can exceed the speed of light. 



. , v St sind 

A0 = ^— (CD 



and the times of the observations recorded by the observer are, 

, d + v St cos 9 

A = h + - 



c 

t' 2 =t 2 + - (C.2) 

c 

Then, the measured interval between the observations is: 

, , v St cos 6 

At =t 2 - t\ = t 2 - h 

c 

= St(l -13 cos 6) (C.3) 

where f3 = v/c. Therefore, the transverse velocity inferred by the observer is 

t>T d A(f> vsin9 j3 sin 9 , . 

^ T = ~ = c At = c(l-Pcos0) = l-/?cos^ ^ ' ' 

We can notice that /3t — ► oo if 6* is small and — > c. 

Making equal to zero the differentiation of eq IC.4I with respect to 9, we obtain 
at what value of 9, /?t is maximized: # max = cos -1 (3. Inserting this result into 
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eqElwe get that /?^ ax = (3Y, where V is the Lorentz factor, T = (1 - (3 2 )- 1 / 2 . 
This is another way to realize that /3p can be arbitrarily high, since if (3 — > 1, (3 « T. 

We show in Fig. 10.21 the measured transverse velocity /3j as a function of the 
angle to the line of sight, 9, for different values of the Lorentz factor, T. 

The situation of an approaching source have been analyzed up to here. The case 
of a receding source yields the following result: 



receding 



d A0 v sin 9 (3 sin 9 



(C.5) 



receding 



c 



c At c(l + /9 cos 6») 1 + (3 cos 6> 
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